Changeset 3127 for trunk/LMDZ.MARS/libf/phymars/improvedclouds_mod.F

Timestamp:

Nov 17, 2023, 11:24:09 AM (12 months ago)

Author:

csegonne

Message:

MARS PCM:

• Bug fixes for improvedcloud_mod.F :
• Initialization of dMice when low ccn number
• Initialization of zdq at first sub time step (spenttime=0)
• Move 'coef' to fix the formula for zimicro

CS

File:

• trunk/LMDZ.MARS/libf/phymars/improvedclouds_mod.F (modified) (3 diffs)

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

@@ -516 +516 @@
           endif
          
+          ELSE
+          ! Initialization of dMice when it's not computed 
+            dMice=0 ! no condensation/sublimation to account for
           ENDIF !of if Nccn>1
 
@@ -570 +573 @@
           IF (cloud_adapt_ts) then
 c           Estimation of how much is actually condensing/subliming
-                zdq=(dMice/spenttime)*(ptimestep-spenttime)
+                IF (spenttime.ne.0) then
+                   zdq=(dMice/spenttime)*(ptimestep-spenttime)
+                ELSE
+                !Initialization for spenttime=0
+                   zdq=zpotcond(ig,l)*((ptimestep-spenttime)/ 
+     &          ptimestep)
+                ENDIF
 c           Threshold with powerlaw (sanity check)
                 zdq=min(abs(zdq),
@@ -689 +698 @@
        alpha=1.81846504e+11
        beta=1.54550140e+00
-       zimicro=ceiling(min(max(coef*alpha*abs(potcond)**beta,5.),7000.))
+       zimicro=ceiling(coef*min(max(alpha*abs(potcond)**beta,5.),7000.))
 
       END SUBROUTINE adapt_imicro