Changeset 5926 for LMDZ6/trunk/libf/phylmd/StratAer/coagulate.f90

Timestamp:

Dec 11, 2025, 4:17:28 PM (5 weeks ago)

Author:

lebasn

Message:

StratAer?: bugfix on last commit + SAI mode

File:

• LMDZ6/trunk/libf/phylmd/StratAer/coagulate.f90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/StratAer/coagulate.f90

@@ -129 +129 @@
   IF (is_strato(ilon,ilev)) THEN
   !compute actual wet particle radius & volume for every grid box
-  IF(flag_strat_compo) THEN
+  IF(flag_strat_compo.EQ.2) THEN
      DO i=1, nbtr_bin
         radiuswet(i)=f_r_wetB(ilon,ilev,i)*mdw(i)/2.
@@ -174 +174 @@
 !--pre-compute the thermal velocity of a particle thvelpar(i) from equation 20
   thvelpar=0.0
-  IF(flag_strat_compo) THEN
+  IF(flag_strat_compo.EQ.2) THEN
      DO i=1, nbtr_bin
         m_par(i)=4./3.*RPI*radiuswet(i)**3.*DENSO4B(ilon,ilev,i)*1000.
@@ -275 +275 @@
 !                     =(Vdry(k)*tr_t(ilon,ilev,k)+pdtcoag*num_dry)/( (1.+pdtcoag*denom)*Vdry(k) ) 
 !          with num_dry=...beta(i,j)*Vdry(i)*....
-!       so in old STRATAER (.not.flag_strat_compo), it was correct
+!       so in old STRATAER (flag_strat_compo.EQ.1), it was correct
   ENDIF