Changeset 517 for LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h

Timestamp:

Apr 16, 2004, 5:43:38 PM (20 years ago)

Author:

lmdzadmin

Message:

Inclusion des modifications de O. Boucher et de J. Quaas pour le calcul des
premiers effets directs et indirects dus aux aerosols
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h (modified) (2 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h

@@ -277 +277 @@
      .                "ave(X)", zsto,zout)
 c
+c Effets des aerosols
+c
+c     IF (ok_ade.OR.ok_aie) THEN
+         CALL histdef(nid_mth, "topsad", "ADE at TOA", "W/m2",
+     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth, "solsad", "ADE at sfc", "W/m2",
+     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth, "topsai", "AIE at TOA", "W/m2",
+     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth, "solsai", "AIE at sfc", "W/m2",
+     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
+     .                "ave(X)", zsto,zout)
+c     endif
+c
+
+c
 c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
@@ -693 +715 @@
      .                "ave(X)", zsto,zout)
          ENDIF
-C
+c
+c Effets des aerosols
+c
+c     IF (ok_ade.OR.ok_aie) THEN
+         CALL histdef(nid_mth, "re", "CDR", "um",
+     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth, "redenom", "CDR denominator", "-",
+     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth, "tau", "cloud opt thickness", "-",
+     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth, "taupi", "cloud opt thickn. (pi)", "-",
+     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c     endif
+c
+         CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
+     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c
          if (nqmax.GE.3) THEN
          DO iq=1,nqmax-2