Changeset 517 for LMDZ.3.3/branches/rel-LF/libf/phylmd/write_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/write_histmth.h (modified) (2 diffs)

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

@@ -262 +262 @@
      $               iim*jjmp1,ndex2d)
 c
+c
+c effets des aerosols
+c
+c     IF (ok_ade.OR.ok_aie) THEN
+      zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
+c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topswad,zx_tmp_2d)
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
+      CALL histwrite(nid_mth,"topsad",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
+c
+      zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
+c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solswad,zx_tmp_2d)
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
+      CALL histwrite(nid_mth,"solsad",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
+c
+      zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
+c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topswai,zx_tmp_2d)
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
+      CALL histwrite(nid_mth,"topsai",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
+c
+      zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
+c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solswai,zx_tmp_2d)
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
+      CALL histwrite(nid_mth,"solsai",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
+c     endif
+c
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
       CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
@@ -712 +737 @@
       ENDIF
 C
+c
+c effets des aerosols
+c
+c     IF (ok_ade.OR.ok_aie) THEN
+      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, re, zx_tmp_3d)
+      CALL histwrite(nid_mth,"re",itau_w,zx_tmp_3d,
+     .                                   iim*jjmp1*klev,ndex3d)
+c
+      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fl, zx_tmp_3d)
+      CALL histwrite(nid_mth,"redenom",itau_w,zx_tmp_3d,
+     .                                   iim*jjmp1*klev,ndex3d)
+c
+      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtau, zx_tmp_3d)
+      CALL histwrite(nid_mth,"tau",itau_w,zx_tmp_3d,
+     .                                   iim*jjmp1*klev,ndex3d)
+c
+      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtaupi, zx_tmp_3d)
+      CALL histwrite(nid_mth,"taupi",itau_w,zx_tmp_3d,
+     .                                   iim*jjmp1*klev,ndex3d)
+c     endif
+c
       IF (nqmax.GE.3) THEN
       DO iq=1,nqmax-2