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

Timestamp:

Oct 30, 2016, 4:35:25 PM (8 years ago)

Author:

oboucher

Message:

Adding a module for stratospheric aerosols with a bin scheme.
The module gets activated with -strataer true compiling option.
May not quite work yet, more testing needed, but should not affect
the rest of LMDz as everything is under a CPP_StratAer key.

File:

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

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -745 +745 @@
     REAL conv_q(klon,klev) ! convergence de l'humidite (kg/kg/s)
     REAL conv_t(klon,klev) ! convergence de la temperature(K/s)
-
-    !
-    !  REAL zxsnow(klon)
+    !
     REAL zxsnow_dummy(klon)
     REAL zsav_tsol(klon)
@@ -763 +761 @@
     real zqsat(klon,klev)
     !
-    INTEGER i, k, iq, ig, j, nsrf, ll, l, iiq
+    INTEGER i, k, iq, j, nsrf, ll, l
     !
     REAL t_coup
@@ -884 +882 @@
     !IM 141004 END
     !IM 190504 BEG
-    INTEGER ij
     !  INTEGER imp1jmp1
     !  PARAMETER(imp1jmp1=(iim+1)*jjmp1)
@@ -893 +890 @@
     LOGICAL ok_msk
     REAL msk(klon)
-    !IM 
-    REAL airetot, pi
     !ym A voir plus tard
     !ym      REAL zm_wo(jjmp1, klev)
@@ -932 +927 @@
     !$OMP THREADPRIVATE(ok_sync)
     real date0
-    integer idayref
 
     ! essai writephys
@@ -953 +947 @@
     DATA      ip_ebil/0/
     !$OMP THREADPRIVATE(ip_ebil)
-    INTEGER   if_ebil ! level for energy conserv. dignostics
-    SAVE      if_ebil
-    !$OMP THREADPRIVATE(if_ebil)
     REAL q2m(klon,nbsrf)  ! humidite a 2m
 
@@ -3483 +3474 @@
           ELSE
 #ifdef CPP_RRTM
+#ifndef CPP_StratAer
+          !--prescribed strat aerosols 
+          !--only in the case of non-interactive strat aerosols
             IF (flag_aerosol_strat.EQ.1) THEN
              CALL readaerosolstrato1_rrtm(debut)
@@ -3492 +3486 @@
              CALL abort_physic(modname,abort_message,1)
             ENDIF
+#endif
 #else
              abort_message='You should compile with -rrtm if running ' &
@@ -3499 +3494 @@
           ENDIF
        ENDIF
+!
+#ifdef CPP_RRTM
+#ifdef CPP_StratAer
+       !--interactive strat aerosols
+       CALL calcaerosolstrato_rrtm(pplay,t_seri,paprs,debut)
+#endif
+#endif
        !--fin STRAT AEROSOL
        !