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

Timestamp:

Dec 9, 2025, 3:08:05 PM (8 hours ago)

Author:

lebasn

Message:

StratAer?: remove old params and rename others for coherence (flag_nucl and flag_strat_compo). Use integer instead of boolean to choose compo and nucleation methods.

File:

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

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

@@ -28 +28 @@
   USE infotrac_phy
   USE phys_local_var_mod, ONLY: DENSO4, DENSO4B, f_r_wet, f_r_wetB
-  USE strataer_local_var_mod, ONLY: flag_new_strat_compo
+  USE strataer_local_var_mod, ONLY: flag_strat_compo
 
   IMPLICIT NONE
@@ -129 +129 @@
   IF (is_strato(ilon,ilev)) THEN
   !compute actual wet particle radius & volume for every grid box
-  IF(flag_new_strat_compo) THEN
+  IF(flag_strat_compo) 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_new_strat_compo) THEN
+  IF(flag_strat_compo) 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_new_strat_compo), it was correct
+!       so in old STRATAER (.not.flag_strat_compo), it was correct
   ENDIF