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

Timestamp:

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

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

@@ -14 +14 @@
   !$OMP THREADPRIVATE(flag_emit_distrib)
   
-  ! flag to choose nucleation nucleation method
-  LOGICAL,SAVE :: flag_new_nucl   ! T=new routine from A. Maattanen (LATMOS), F=older routine from H. Vehkamäki (FMI)
-  !$OMP THREADPRIVATE(flag_new_nucl)
-  
-  ! Use relative humidity from 2D model stratospheric H2O because LMDz is too dry in the stratosphere
-  ! (no CH4 oxidation)
-  LOGICAL,SAVE :: flag_H2O2d_nucleation
-  !$OMP THREADPRIVATE(flag_H2O2d_nucleation)
+  ! flag to choose nucleation method
+  INTEGER,SAVE :: flag_nucl   ! 1: routine from H. Vehkamäki (FMI), 2: routine from A. Maattanen (LATMOS)
+  !$OMP THREADPRIVATE(flag_nucl)
   
   ! OH reduction from SO2. OH is reduced when its reaction with SO2 competes sufficiently with its reaction
@@ -33 +28 @@
   !$OMP THREADPRIVATE(flag_H2SO4_photolysis)
   
-  ! flag for minimum lifetime (=1.5 pdt phys)
-  LOGICAL,SAVE :: flag_min_rreduce
-  !$OMP THREADPRIVATE(flag_min_rreduce)
-  
-  ! flag to read new climato (O3, H2O & H2SO4_LIFET)
-  LOGICAL,SAVE :: flag_newclim_file
-  !$OMP THREADPRIVATE(flag_newclim_file)
-  
-  ! flag to choose new H2SO4 density and weight percent from Tabazadeh et al. (1994).
-  LOGICAL,SAVE :: flag_new_strat_compo
-  !$OMP THREADPRIVATE(flag_new_strat_compo)
+  ! Parameterization method to compute H2SO4/H2O aerosol composition
+  INTEGER,SAVE :: flag_strat_compo  !  1: S. Bekki et al., 2: Tabazadeh et al. 1997
+  !$OMP THREADPRIVATE(flag_strat_compo)
   
   ! Verbose mode to get more print info
@@ -184 +171 @@
     flag_emit = 0                   ! Background (default)
     flag_emit_distrib = 0           ! Gaussian (default)
-    flag_new_nucl = .TRUE.          ! Define nucleation routine (default: A. Maattanen - LATMOS)
+    flag_nucl = 2                   ! Define nucleation routine (default: A. Maattanen - LATMOS)
     flag_verbose_strataer = .FALSE. ! verbose mode
-    flag_newclim_file = .TRUE.      ! Define input climato file (default: all climato)
-    flag_H2O2d_nucleation = .FALSE. ! Use H2O 2D climato (default: No)
     flag_OH_reduced = .FALSE.       ! OH reduce (default: No)
     flag_H2SO4_photolysis = .FALSE. ! H2SO4 photolysis (default: No)
-    flag_min_rreduce = .TRUE.       ! Minimum lifetime=1.5 pdt phys (default: Yes)
-    flag_new_strat_compo =.FALSE.   ! H2SO4/H2O weight percent & density routine (default: S. Bekki)
+    flag_strat_compo = 2            ! H2SO4/H2O composition routine (default: Tabazadeh et al. 1997)
     ok_qemiss = .FALSE.             ! H2O emission flag
     
@@ -219 +203 @@
     CALL getin_p('flag_emit_distrib',flag_emit_distrib)
     CALL getin_p('flag_verbose_strataer',flag_verbose_strataer)
-    CALL getin_p('flag_new_nucl',flag_new_nucl)
-    CALL getin_p('flag_newclim_file',flag_newclim_file)
-    CALL getin_p('flag_H2O2d_nucleation',flag_H2O2d_nucleation)
+    CALL getin_p('flag_nucl',flag_nucl)
     CALL getin_p('flag_OH_reduced',flag_OH_reduced)
     CALL getin_p('flag_H2SO4_photolysis',flag_H2SO4_photolysis)
-    CALL getin_p('flag_min_rreduce',flag_min_rreduce)
-    CALL getin_p('flag_new_strat_compo',flag_new_strat_compo)
+    CALL getin_p('flag_strat_compo',flag_strat_compo)
     CALL getin_p('ok_qemiss',ok_qemiss)
-    
-    !============= Test flag coherence =============
-    IF (.NOT. flag_newclim_file) THEN
-       IF (flag_H2SO4_photolysis .OR. flag_OH_reduced .OR. flag_H2O2d_nucleation) THEN
-          WRITE(lunout,*) 'ERROR : flag_newclim_file=',flag_newclim_file, &
-               ' whereas flag_H2SO4_photolysis=',flag_H2SO4_photolysis,', flag_OH_reduced=',flag_OH_reduced, &
-               ' and flag_H2O2d_nucleation=',flag_H2O2d_nucleation
-          CALL abort_physic('strataer_local_var_mod','Incompatible options in physiq_def file !',1)
-       ENDIF
-       IF(flag_min_rreduce) THEN
-          WRITE(lunout,*) 'Warning : flag_min_rreduce will be ignored with old climato file !'
-       ENDIF
-    ENDIF
     
     !============= Print params =============
     IF (is_master) THEN
        WRITE(lunout,*) 'flag_emit = ',flag_emit
-       WRITE(lunout,*) 'IN STRATAER : flag_new_nucl = ',flag_new_nucl
-       WRITE(lunout,*) 'IN STRATAER : flag_newclim_file = ',flag_newclim_file
+       WRITE(lunout,*) 'IN STRATAER : flag_nucl = ',flag_nucl
        WRITE(lunout,*) 'IN STRATAER : flag_emit_distrib = ',flag_emit_distrib
        WRITE(lunout,*) 'IN STRATAER : flag_verbose_strataer = ',flag_verbose_strataer
        IF (flag_emit == 1 .OR. flag_emit == 4) THEN
-          WRITE(lunout,*) 'IN STRATAER : flag_H2O2d_nucleation = ',flag_H2O2d_nucleation
           WRITE(lunout,*) 'IN STRATAER : flag_OH_reduced = ',flag_OH_reduced
           WRITE(lunout,*) 'IN STRATAER : flag_H2SO4_photolysis = ',flag_H2SO4_photolysis
-          WRITE(lunout,*) 'IN STRATAER : flag_min_rreduce = ',flag_min_rreduce
-          WRITE(lunout,*) 'IN STRATAER : flag_new_strat_compo = ',flag_new_strat_compo
+          WRITE(lunout,*) 'IN STRATAER : flag_strat_compo = ',flag_strat_compo
           WRITE(lunout,*) 'IN STRATAER : ok_qemiss = ',ok_qemiss
        ENDIF