Changeset 2536 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Jun 3, 2016, 3:04:40 PM (8 years ago)

Author:

oboucher

Message:

Introducing stratomask diagnostic on where the stratosphere is
flag_aerosol_strat = 2 for CMIP6 strat aerosol forcing

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -3079 +3079 @@
              CALL readaerosolstrato(debut)
            ELSE
-             abort_message='flag_aerosol_strat must equal 1 for rrtm=0' 
+             abort_message='flag_aerosol_strat must equal 1 for rrtm=0'
              call abort_physic(modname,abort_message,1)
            ENDIF
           ELSE
 #ifdef CPP_RRTM
-             CALL readaerosolstrato_rrtm(debut)
+            IF (flag_aerosol_strat.EQ.1) THEN
+             CALL readaerosolstrato1_rrtm(debut)
+            ELSEIF (flag_aerosol_strat.EQ.2) THEN
+             CALL stratosphere_mask(t_seri, pplay, latitude_deg)
+             CALL readaerosolstrato2_rrtm(debut)
+            ELSE
+             abort_message='flag_aerosol_strat must equal 1 or 2 for rrtm=1'
+             call abort_physic(modname,abort_message,1)
+            ENDIF
 #else
              abort_message='You should compile with -rrtm if running ' &
@@ -3093 +3101 @@
        ENDIF
        !--fin STRAT AEROSOL
-
 
        !   On prend la somme des fractions nuageuses et des contenus en eau