Changeset 5924 for LMDZ6/trunk/libf/phylmd/StratAer/so2_to_h2so4.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/so2_to_h2so4.f90 (modified) (3 diffs)

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

@@ -10 +10 @@
   ! lifetime (sec) et O3_clim (VMR)
   USE phys_local_var_mod, ONLY : SO2_lifetime,H2SO4_lifetime,O3_clim,budg_3D_so2_to_h2so4,budg_so2_to_h2so4,SO2_chlm
-  USE strataer_local_var_mod, ONLY : flag_OH_reduced, flag_H2SO4_photolysis, flag_min_rreduce
+  USE strataer_local_var_mod, ONLY : flag_OH_reduced, flag_H2SO4_photolysis
   
   IMPLICIT NONE
@@ -100 +100 @@
               ! Check lifetime rreduce < timestep*1.5 (such as SO2 loss > 0.5*SO2) with exp(-1/1.5)=0.52
               ! Check lifetime rreduce < timestep*3 (such as SO2 loss > 0.28*SO2) with exp(-1/3)=0.72
-              IF(flag_min_rreduce) THEN
-                 IF (rreduce .LT. (3.*pdtphys)) rreduce = 3.*pdtphys
-              ENDIF
+              IF (rreduce .LT. (3.*pdtphys)) rreduce = 3.*pdtphys
               budg_3D_so2_to_h2so4(ilon,ilev)=tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/rreduce))
               tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) - budg_3D_so2_to_h2so4(ilon,ilev)
@@ -132 +130 @@
                  ! Check lifetime rreduce < timestep*1.5 (such as H2SO4 loss > 0.5*H2SO4) with exp(-1/1.5)=0.52
                  ! Check lifetime rreduce < timestep*3 (such as H2SO4 loss > 0.28*H2SO4) with exp(-1/3)=0.72
-                 IF(flag_min_rreduce) THEN
-                    IF (rreduce .LT. (3.*pdtphys)) rreduce = 3.*pdtphys
-                 ENDIF
+                 IF (rreduce .LT. (3.*pdtphys)) rreduce = 3.*pdtphys
                  dummyso4toso2 = (mSO2mol/mH2SO4mol)*tr_seri(ilon,ilev,id_H2SO4_strat)*(1.0-exp(-pdtphys/rreduce))
                  budg_3D_so2_to_h2so4(ilon,ilev) = budg_3D_so2_to_h2so4(ilon,ilev) + dummyso4toso2