Changeset 411 for trunk/LMDZ.MARS/libf/phymars/suaer.F90

Timestamp:

Nov 22, 2011, 5:34:42 PM (13 years ago)

Author:

tnavarro

Message:

changed scavenging in improvedclouds.F, updated ice radius in callsedim.F and commented outputs in suaer.F90 & aeropacity.F

File:

• trunk/LMDZ.MARS/libf/phymars/suaer.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/suaer.F90

@@ -376 +376 @@
   END DO
 ! 2.5 Output display
-  WRITE(*,*) 'Les donnees spectrales :'
-  WRITE(*,*) 'Solaire (SW) ---->'
-  WRITE(*,*) 'Aerosol number: ', iaer
-  WRITE(*,*) 'Rayon aerosol: ', radiustab(iaer,idomain,isize)
-  WRITE(*,*) '<Qext>/Qext(longrefvis) ; omega ; g'
-  DO isun=1,nsun
-    WRITE(*,*) QVISsQREF(isun,iaer,isize),&
-         omegavis(isun,iaer,isize),&
-         gvis(isun,iaer,isize)
-  ENDDO
-  WRITE(*,*) 'QREFvis(',iaer,isize,') = ',QREFvis(iaer,isize)
-  WRITE(*,*) 'omegaREFvis(',iaer,isize,') = ',&
-                                      omegaREFvis(iaer,isize)
+!  WRITE(*,*) 'Les donnees spectrales :'
+!  WRITE(*,*) 'Solaire (SW) ---->'
+!  WRITE(*,*) 'Aerosol number: ', iaer
+!  WRITE(*,*) 'Rayon aerosol: ', radiustab(iaer,idomain,isize)
+!  WRITE(*,*) '<Qext>/Qext(longrefvis) ; omega ; g'
+!  DO isun=1,nsun
+!    WRITE(*,*) QVISsQREF(isun,iaer,isize),&
+!         omegavis(isun,iaer,isize),&
+!         gvis(isun,iaer,isize)
+!  ENDDO
+!  WRITE(*,*) 'QREFvis(',iaer,isize,') = ',QREFvis(iaer,isize)
+!  WRITE(*,*) 'omegaREFvis(',iaer,isize,') = ',&
+!                                      omegaREFvis(iaer,isize)
 ! ------------------------------------------------
 ENDDO
@@ -427 +427 @@
        nir-1,longir,epav,omegav,gav,&
        QREFir(iaer,isize),omegaREFir(iaer,isize) )
-  WRITE(*,*) 'QREFir(',iaer,isize,') = ',QREFir(iaer,isize)
-  WRITE(*,*) 'omegaREFir(',iaer,isize,') = ',&
-                                      omegaREFir(iaer,isize)
+!  WRITE(*,*) 'QREFir(',iaer,isize,') = ',QREFir(iaer,isize)
+!  WRITE(*,*) 'omegaREFir(',iaer,isize,') = ',&
+!                                      omegaREFir(iaer,isize)
 
 ! 2.5 Computing  <QIR>/Qext(longrefvis)
 
   DO iir=1,nir-1
-    WRITE(*,*) 'QIRsQREFir Channel ',iir,': ',epav(iir)
+!    WRITE(*,*) 'QIRsQREFir Channel ',iir,': ',epav(iir)
     epav(iir)=  epav(iir) * QREFir(iaer,isize) / &
                             QREFvis(iaer,isize)
   ENDDO
-  WRITE(*,*) 'Aerosol number', iaer
-  WRITE(*,*) 'Particle size: ',radiustab(iaer,idomain,isize)
-  WRITE(*,*) 'Rapport Solaire/IR:',&
-             QREFvis(iaer,isize) / QREFir(iaer,isize)
+!  WRITE(*,*) 'Aerosol number', iaer
+!  WRITE(*,*) 'Particle size: ',radiustab(iaer,idomain,isize)
+!  WRITE(*,*) 'Rapport Solaire/IR:',&
+!             QREFvis(iaer,isize) / QREFir(iaer,isize)
 
 ! 2.6 Variable assignements
@@ -473 +473 @@
 ! 2.7 Output display
 
-  WRITE(*,*) 'AEROSOL PROPERTIES: Number ',iaer
-  WRITE(*,*) 'Thermal IR (LW) ---->'
-  WRITE(*,*) 'Particle size: ',radiustab(iaer,idomain,isize)
-  WRITE(*,*) '<Qext>/Qext(longrefvis) ; omega ; g'
-  DO iir=1,nir
-    WRITE(*,*) QIRsQREF(iir,iaer,isize),omegaIR(iir,iaer,isize),&
-         gIR(iir,iaer,isize)
-  ENDDO
-  WRITE(*,*) 'CO2: <Qabs>/Qext(longrefvis) = ',&
-       QIRsQREF(1,iaer,isize)*(1-omegaIR(1,iaer,isize))
-  WRITE(*,*) ''
+!  WRITE(*,*) 'AEROSOL PROPERTIES: Number ',iaer
+!  WRITE(*,*) 'Thermal IR (LW) ---->'
+!  WRITE(*,*) 'Particle size: ',radiustab(iaer,idomain,isize)
+!  WRITE(*,*) '<Qext>/Qext(longrefvis) ; omega ; g'
+!  DO iir=1,nir
+!    WRITE(*,*) QIRsQREF(iir,iaer,isize),omegaIR(iir,iaer,isize),&
+!         gIR(iir,iaer,isize)
+!  ENDDO
+!  WRITE(*,*) 'CO2: <Qabs>/Qext(longrefvis) = ',&
+!       QIRsQREF(1,iaer,isize)*(1-omegaIR(1,iaer,isize))
+!  WRITE(*,*) ''
 
 ENDDO ! isize (particle size) -------------------------------------