Changeset 455

Timestamp:

Dec 6, 2011, 4:40:55 PM (13 years ago)

Author:

tnavarro

Message:

tiny change : nuice_sed initialisation is now done in inifis. Also changed initracer and improvedclouds.

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/improvedclouds.F (modified) (5 diffs)
• libf/phymars/inifis.F (modified) (1 diff)
• libf/phymars/initracer.F (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -1299 +1299 @@
 
 
-
-
-
+== 06/12/11 == TN
+
+>> tiny change : nuice_sed initialisation is now done in inifis. Also changed initracer and improvedclouds.

trunk/LMDZ.MARS/libf/phymars/improvedclouds.F

@@ -126 +126 @@
       REAL Mccn_col(ngridmx)         ! total column ccn mass
       REAL Nccn_col(ngridmx)         ! total column ccn mass
+      REAL rate_out(ngridmx,nlayermx) ! nucleation rate
       INTEGER count
       
@@ -191 +192 @@
         
         write(*,*) 'Variance of ice & CCN distribs :', sigma_ice
+        write(*,*) 'nuice for sedimentation:', nuice_sed
         write(*,*) 'Volume of a water molecule:', vo1
 
@@ -316 +318 @@
         dN = 0.
         dM = 0.
+        rate_out(ig,l) = 0
         do i = 1, nbin_cld
           n_aer(i) = n_aer(i) / ( 1. + rate(i)*ptimestep )
@@ -321 +324 @@
           dN       = dN + n_aer(i) * rate(i) * ptimestep
           dM       = dM + m_aer(i) * rate(i) * ptimestep
+          rate_out(ig,l)=rate_out(ig,l)+rate(i)
         enddo
         
@@ -539 +543 @@
          call WRITEDIAGFI(ngrid,"growthrate","growth rate","m^2/s",1,
      &                    gr_out)
+         call WRITEDIAGFI(ngrid,"nuclearate","nucleation rate","",1,
+     &                    rate_out)
          call WRITEDIAGFI(ngrid,"dM","ccn variation","kg",1,
      &                    dM_out)

trunk/LMDZ.MARS/libf/phymars/inifis.F

@@ -419 +419 @@
          endif
          
+! water ice clouds effective variance distribution for sedimentaion       
+        write(*,*) "effective variance for water ice clouds ?"
+        nuice_sed=0.45 
+        call getin("nuice_sed",nuice_sed)
+        write(*,*) "water_param nueff Sedimentation:", nuice_sed
+         
 ! ccn factor if no scavenging         
         write(*,*) "water param CCN reduc. factor ?", ccn_factor

trunk/LMDZ.MARS/libf/phymars/initracer.F

@@ -496 +496 @@
       nuice_ref=0.1   ! Effective variance nueff of the
                       ! water-ice size distribution
-      nuice_sed=0.45   ! Sedimentation effective variance
+      !!!nuice_sed=0.45   ! Sedimentation effective variance
                       ! of the water-ice size distribution