Changeset 903 for LMDZ4

Timestamp:

Feb 7, 2008, 3:32:33 PM (17 years ago)

Author:

Laurent Fairhead

Message:

On avait perdu les sorties aerosols de Johannes dans un merge
LF

Location:

LMDZ4/trunk/libf/phylmd

Files:

• ini_histmth.h (modified) (2 diffs)
• write_histmth.h (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/ini_histmth.h

@@ -300 +300 @@
      .                "ave(X)", zstophy,zout)
 c
+c Effets des aerosols
+c
+      IF (ok_ade.OR.ok_aie) THEN
+         CALL histdef(nid_mth, "topsad", "ADE at TOA", "W/m2",
+     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
+     .                "ave(X)", zstophy,zout)
+c
+         CALL histdef(nid_mth, "solsad", "ADE at sfc", "W/m2",
+     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
+     .                "ave(X)", zstophy,zout)
+c
+         CALL histdef(nid_mth, "topsai", "AIE at TOA", "W/m2",
+     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
+     .                "ave(X)", zstophy,zout)
+c
+         CALL histdef(nid_mth, "solsai", "AIE at sfc", "W/m2",
+     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
+     .                "ave(X)", zstophy,zout)
+      endif
+c
 c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
 c    .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
@@ -1303 +1323 @@
      .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
      .                "ave(X)", zstophy,zout)
+c
+c Effets des aerosols
+c
+         IF (ok_ade.OR.ok_aie) THEN
+             CALL histdef(nid_mth, "topsad", "ADE at TOA", "W/m2",
+     .                    iim,jj_nb,nhori, 1,1,1, -99, 32,
+     .                    "ave(X)", zsto,zout)
+c
+             CALL histdef(nid_mth, "solsad", "ADE at sfc", "W/m2",
+     .                    iim,jj_nb,nhori, 1,1,1, -99, 32,
+     .                    "ave(X)", zsto,zout)
+c
+             CALL histdef(nid_mth, "topsai", "AIE at TOA", "W/m2",
+     .                    iim,jj_nb,nhori, 1,1,1, -99, 32,
+     .                    "ave(X)", zsto,zout)
+c
+             CALL histdef(nid_mth, "solsai", "AIE at sfc", "W/m2",
+     .                    iim,jj_nb,nhori, 1,1,1, -99, 32,
+     .                    "ave(X)", zsto,zout)
+         ENDIF
 c
 c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",

LMDZ4/trunk/libf/phylmd/write_histmth.h

@@ -1273 +1273 @@
 cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
       CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
+c effets des aerosols
+c
+      IF (ok_ade.OR.ok_aie) THEN
+      zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
+      CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_2d)
+c
+      zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
+      CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_2d)
+c
+      zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
+      CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_2d)
+c
+      zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
+      CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_2d)
+      endif
+
 c
 cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
@@ -1284 +1300 @@
 cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
       CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
+c
+c Effets des aerosols
+c
+      IF (ok_ade.OR.ok_aie) THEN
+        zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
+        CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d)
+c
+        zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
+        CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_fi2d)
+c
+        zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
+        CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d)
+c
+        zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
+        CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d)
+      ENDIF
 c
 c