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

Timestamp:

Dec 9, 2025, 3:08:05 PM (10 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/traccoag_mod.f90 (modified) (2 diffs)

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

@@ -121 +121 @@
     WHERE (stratomask.GT.0.5) is_strato=.TRUE. 
 
-    IF(flag_new_strat_compo) THEN
-       IF(debutphy) WRITE(lunout,*) 'traccoag: COMPO/DENSITY (Tabazadeh 97) + H2O kelvin effect', flag_new_strat_compo
+    IF(flag_strat_compo) THEN
+       IF(debutphy) WRITE(lunout,*) 'traccoag: COMPO/DENSITY (Tabazadeh 97) + H2O kelvin effect', flag_strat_compo
        ! STRACOMP (H2O, P, t_seri, R -> R2SO4 + Kelvin effect) : Taba97, Socol, etc...
        CALL stracomp_kelvin(sh,t_seri,pplay)
     ELSE
-       IF(debutphy) WRITE(lunout,*) 'traccoag: COMPO from Bekki 2D model', flag_new_strat_compo
+       IF(debutphy) WRITE(lunout,*) 'traccoag: COMPO from Bekki 2D model', flag_strat_compo
        ! STRACOMP (H2O, P, t_seri -> aerosol composition (R2SO4)) 
        ! H2SO4 mass fraction in aerosol (%)
@@ -383 +383 @@
                    & *pplay(i,j)/t_seri(i,j)/RD                           ! [air mass concentration in kg air /m3A]
               
-              IF(flag_new_strat_compo) THEN
+              IF(flag_strat_compo) THEN
                  !     SAD_sulfate: SAD WET sulfate aerosols (cm2/cm3)
                  SAD_sulfate(i,j) = SAD_sulfate(i,j) + nd_mode(i,j,it) &     ! [DRY part/m3A (in bin it)]