Changeset 3412 for LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

Timestamp:

Nov 8, 2018, 10:39:16 PM (6 years ago)

Author:

oboucher

Message:

Externalizing a controlling parameter for activating aerosol radiative feedback
Need to pass the info all the way to recmwf_aero. Default is True as before.

File:

• LMDZ6/trunk/libf/phylmd/conf_phys_m.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

@@ -18 +18 @@
        iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
        ok_ade, ok_aie, ok_alw, ok_cdnc, aerosol_couple, chemistry_couple, &
-       flag_aerosol, flag_aerosol_strat, new_aod, &
+       flag_aerosol, flag_aerosol_strat, flag_aer_feedback, new_aod, &
        flag_bc_internal_mixture, bl95_b0, bl95_b1,&
        read_climoz, &
@@ -74 +74 @@
     INTEGER              :: flag_aerosol
     INTEGER              :: flag_aerosol_strat
+    LOGICAL              :: flag_aer_feedback
     LOGICAL              :: flag_bc_internal_mixture
     LOGICAL              :: new_aod
@@ -95 +96 @@
     INTEGER, SAVE       :: flag_aerosol_omp
     INTEGER, SAVE       :: flag_aerosol_strat_omp
+    LOGICAL, SAVE       :: flag_aer_feedback_omp
     LOGICAL, SAVE       :: flag_bc_internal_mixture_omp
     LOGICAL, SAVE       :: new_aod_omp
@@ -1022 +1024 @@
       CALL getin('ok_suntime_rrtm',ok_suntime_rrtm_omp)
     ENDIF
-    !
+    
     !Config Key  = flag_aerosol_strat
     !Config Desc = use stratospheric aerosols 0, 1, 2
@@ -1032 +1034 @@
     !Config Help = Used in physiq.F
     !
-    !
     flag_aerosol_strat_omp = 0
     CALL getin('flag_aerosol_strat',flag_aerosol_strat_omp)
 
-    !
+    !Config Key  = flag_aer_feedback
+    !Config Desc = (des)activate aerosol radiative feedback
+    ! - F = no aerosol radiative feedback
+    ! - T = aerosol radiative feedback
+    !Config Def  = T
+    !Config Help = Used in physiq.F
+    !
+    flag_aer_feedback_omp = .TRUE. 
+    IF (iflag_rrtm_omp==1) THEN
+       CALL getin('flag_aer_feedback',flag_aer_feedback_omp)
+    ENDIF
+
     !Config Key  = iflag_cld_th 
     !Config Desc =  
@@ -2297 +2309 @@
     flag_aerosol=flag_aerosol_omp
     flag_aerosol_strat=flag_aerosol_strat_omp
+    flag_aer_feedback=flag_aer_feedback_omp
     flag_bc_internal_mixture=flag_bc_internal_mixture_omp
     new_aod=new_aod_omp
@@ -2635 +2648 @@
     WRITE(lunout,*) ' flag_aerosol = ', flag_aerosol
     WRITE(lunout,*) ' flag_aerosol_strat= ', flag_aerosol_strat
+    WRITE(lunout,*) ' flag_aer_feedback= ', flag_aer_feedback
     WRITE(lunout,*) ' new_aod = ', new_aod
     WRITE(lunout,*) ' aer_type = ',aer_type