Changeset 2644 for LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

Timestamp:

Oct 2, 2016, 6:55:08 PM (8 years ago)

Author:

oboucher

Message:

I introduced flag_bc_internal_mixture for BC/sulphate internal mixture.
Only works for iflag_rrtm=y, NSW=6 and flag_aerosol=6 or aerosol_couple.
It has no impact at the moment as LUT for aerosol properties will have to
be changed. But the plumbing work is already done....

File:

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

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -19 +19 @@
        ok_ade, ok_aie, ok_cdnc, aerosol_couple, &
        flag_aerosol, flag_aerosol_strat, new_aod, &
-       bl95_b0, bl95_b1,&
+       flag_bc_internal_mixture, bl95_b0, bl95_b1,&
        read_climoz, &
        alp_offset)
@@ -64 +64 @@
     ! ok_cdnc, ok cloud droplet number concentration
     ! flag_aerosol_strat : flag pour les aerosols stratos
+    ! flag_bc_internal_mixture : use BC internal mixture if true
     ! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc 
     !
@@ -77 +78 @@
     INTEGER              :: flag_aerosol
     INTEGER              :: flag_aerosol_strat
+    LOGICAL              :: flag_bc_internal_mixture
     LOGICAL              :: new_aod
     REAL                 :: bl95_b0, bl95_b1
@@ -95 +97 @@
     INTEGER, SAVE       :: flag_aerosol_omp
     INTEGER, SAVE       :: flag_aerosol_strat_omp
+    LOGICAL, SAVE       :: flag_bc_internal_mixture_omp
     LOGICAL, SAVE       :: new_aod_omp
     REAL,SAVE           :: bl95_b0_omp, bl95_b1_omp
@@ -398 +401 @@
     flag_aerosol_strat_omp = 0
     CALL getin('flag_aerosol_strat',flag_aerosol_strat_omp)
+
+    !
+    !Config Key  = flag_bc_internal_mixture
+    !Config Desc = state of mixture for BC aerosols
+    ! - n = external mixture
+    ! - y = internal mixture
+    !Config Def  = n
+    !Config Help = Used in physiq.F / aeropt
+    !
+    flag_bc_internal_mixture_omp = .false.
+    CALL getin('flag_bc_internal_mixture',flag_bc_internal_mixture_omp)
 
     ! Temporary variable for testing purpose!
@@ -2136 +2150 @@
     flag_aerosol=flag_aerosol_omp
     flag_aerosol_strat=flag_aerosol_strat_omp
+    flag_bc_internal_mixture=flag_bc_internal_mixture_omp
     new_aod=new_aod_omp
     aer_type = aer_type_omp
@@ -2311 +2326 @@
     IF (ok_aie .AND. .NOT. ok_cdnc) THEN
        CALL abort_physic('conf_phys', 'ok_cdnc must be set to y if ok_aie is activated',1)
+    ENDIF
+
+    ! BC internal mixture is only possible with RRTM & NSW=6 & flag_aerosol=6 or aerosol_couple
+    IF (flag_bc_internal_mixture .AND. NSW.NE.6) THEN 
+       CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with NSW=6',1)
+    ENDIF
+    IF (flag_bc_internal_mixture .AND. iflag_rrtm.NE.1) THEN 
+       CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with RRTM',1)
+    ENDIF
+    IF (flag_bc_internal_mixture .AND. flag_aerosol.NE.6) THEN 
+       CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with flag_aerosol=6',1)
     ENDIF