Changeset 1685 for trunk/LMDZ.MARS/libf/phymars/nucleaCO2.F

Timestamp:

Apr 3, 2017, 11:08:39 AM (8 years ago)

Author:

jaudouard

Message:

Further modifications on CO2 clouds scheme. Water ice clouds can now serve as CCN for CO2 clouds

File:

• trunk/LMDZ.MARS/libf/phymars/nucleaCO2.F (modified) (6 diffs)

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

@@ -2 +2 @@
 *                                                     *
       subroutine nucleaCO2(pco2,temp,sat,n_ccn,nucrate,
-     &           n_ccn_h2oice,rad_h2oice,nucrate_h2oice)
+     &           n_ccn_h2oice,rad_h2oice,nucrate_h2oice,
+     &           vo2co2)
       USE comcstfi_h
 
@@ -24 +25 @@
 
 c     Inputs
-      DOUBLE PRECISION pco2,sat
+      DOUBLE PRECISION pco2,sat,vo2co2
       DOUBLE PRECISION n_ccn(nbinco2_cld), n_ccn_h2oice(nbinco2_cld)
       REAL temp
@@ -49 +50 @@
 c     double precision xratio
       
-      double precision mtetalocal ! local mteta in double precision
+      double precision mtetalocal,mtetalocalh ! local mteta in double precision
 
       double precision fshapeco2simple,zefshapeco2
@@ -63 +64 @@
 
       mtetalocal = dble(mtetaco2)  !! use mtetalocal for better performance
+      mtetalocalh=dble(mteta)
 
 cccccccccccccccccccccccccccccccccccccccccccccccccc
@@ -99 +101 @@
 
         nco2   = pco2 / kbz / temp
-        rstar  = 2. * sigco2 * vo1co2 / (kbz*temp*dlog(sat))
-        gstar  = 4. * pi * (rstar * rstar * rstar) / (3.*vo1co2)
+        rstar  = 2. * sigco2 * vo2co2 / (kbz*temp*dlog(sat))
+        gstar  = 4. * pi * (rstar * rstar * rstar) / (3.*vo2co2)
         
        fshapeco2simple = (2.+mtetalocal)*(1.-mtetalocal)*(1.-mtetalocal)
@@ -152 +154 @@
 
           if (rad_h2oice(i).gt.3000.*rstar) then
-            zefshapeco2 = fshapeco2simple
-          else
-            zefshapeco2 = fshapeco2(mtetalocal,rad_h2oice(i)/rstar) ! same m for dust/h2o ice
+            zefshapeco2 = (2.+mtetalocalh)*(1.-mtetalocalh)*
+     &                (1.-mtetalocalh) / 4.
+          else
+            zefshapeco2 = fshapeco2(mtetalocalh,rad_h2oice(i)/rstar) ! same m for dust/h2o ice
           endif