Index: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histmth.h
===================================================================
--- LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histmth.h	(revision 894)
+++ LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histmth.h	(revision 896)
@@ -291,4 +291,25 @@
      .                iim,jj_nb,nhori, 1,1,1, -99, 32, 
      .                "ave(X)", zstophy,zout)
+c
+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,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)
+c     endif
 c
 c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
Index: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histmth.h
===================================================================
--- LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histmth.h	(revision 894)
+++ LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histmth.h	(revision 896)
@@ -244,4 +244,22 @@
 cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
       CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
+c
+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_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)