Changeset 5440 for LMDZ6

Timestamp:

Dec 20, 2024, 8:54:12 PM (7 hours ago)

Author:

evignon

Message:

small bug fixes in lscp

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_lscp.f90 (modified) (1 diff)
• lmdz_lscp_condensation.f90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp.f90

@@ -984 +984 @@
             ENDIF
 
-            zfrac_lessi = 1. - EXP(zalpha_tr*zprec_cond(i)/zneb(i))
+            zfrac_lessi = 1. - EXP(zalpha_tr*zprec_cond(i)/MAX(rneb(i,k),seuil_neb))
             frac_nucl(i,k)= 1.-zneb(i)*zfrac_lessi 
 
             ! Nucleation with a factor of -1 instead of -0.5
-            zfrac_lessi = 1. - EXP(-zprec_cond(i)/zneb(i))
+            zfrac_lessi = 1. - EXP(-zprec_cond(i)/MAX(rneb(i,k),seuil_neb))
 
         ENDIF

LMDZ6/trunk/libf/phylmd/lmdz_lscp_condensation.f90

@@ -292 +292 @@
   !--formed elsewhere)
   IF (keepgoing(i)) THEN
+
+    !--Initialisation
+    issrfra(i)  = 0.
+    qissr(i)    = 0.
 
     !--If the temperature is higher than the threshold below which
@@ -363 +367 @@
       dqi_mix(i)  = 0.
       dqvc_mix(i) = 0.
-
-      issrfra(i)  = 0.
-      qissr(i)    = 0.
-      subfra(i)   = 0.
-      qsub(i)     = 0.
 
       !--Initialisation of the cell properties