Changeset 684 for LMDZ4/trunk/libf/phylmd/write_histday.h

Timestamp:

Apr 4, 2006, 5:00:40 PM (19 years ago)

Author:

lmdzadmin

Message:

Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def +
Lecture freq. ecriture en nombre de jours
IM

File:

• LMDZ4/trunk/libf/phylmd/write_histday.h (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/write_histday.h

@@ -135 +135 @@
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
       CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
+c
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
+      CALL histwrite(nid_day,"radsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
 c
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
@@ -274 +277 @@
         zx_tmp_fi2d(i) = fluxo(i)
        ELSE
-        zx_tmp_fi2d(i) = 1.E+20
+        zx_tmp_fi2d(i) = 0.
        ENDIF
       ENDDO
@@ -286 +289 @@
         zx_tmp_fi2d(i) = fluxg(i)
        ELSE
-        zx_tmp_fi2d(i) = 1.E+20
+        zx_tmp_fi2d(i) = 0.
        ENDIF
       ENDDO
@@ -386 +389 @@
       CALL histwrite(nid_day, "SWdnSFCclr",itau_w,zx_tmp_2d,
      .                               iim*jjmp1,ndex2d)
+c
+cIM AMMA-MIP
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
+      CALL histwrite(nid_day,"LWdnSFC",itau_w,zx_tmp_2d,iim*jjmp1,
+     .               ndex2d)
+c
+      zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
+      CALL histwrite(nid_day,"LWupSFC",itau_w,zx_tmp_2d,iim*jjmp1,
+     .               ndex2d)
+c
+cIM   IF (iflag_con.GE.3) THEN
+cIM ema_pct= 0.     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ema_pct,zx_tmp_2d)
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1, meanptop(:,1),zx_tmp_2d)
+      CALL histwrite(nid_day,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,
+     .               ndex2d)
+cIM   ENDIF
 c
       zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
@@ -578 +598 @@
      .               iim*jjmp1*klev,ndex3d)
 c
+c
+cIM: 101003 : K/30min ==> K/s
+      zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
+      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
+      CALL histwrite(nid_day,"dtcon",itau_w,zx_tmp_3d,
+     .                                   iim*jjmp1*klev,ndex3d)
+c
+      zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
+      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
+      CALL histwrite(nid_day,"dqcon",itau_w,zx_tmp_3d,
+     .                                   iim*jjmp1*klev,ndex3d)
+c
+      DO i=1, klon
+       zx_tmp_2d(i)=MIN(100.,rh2m(i)*100.)
+      ENDDO
+c
+      CALL gr_fi_ecrit(1,klon,iim,jjmp1, zx_tmp_2d,zx_tmp_2d)
+      CALL histwrite(nid_day,"rh2m",itau_w,zx_tmp_2d,
+     .               iim*jjmp1,ndex2d)
+c
+      CALL gr_fi_ecrit(1,klon,iim,jjmp1, qsat2m,zx_tmp_2d)
+      CALL histwrite(nid_day,"qsat2m",itau_w,zx_tmp_2d,
+     .               iim*jjmp1,ndex2d)
+c
+      CALL gr_fi_ecrit(1,klon,iim,jjmp1, tpot,zx_tmp_2d)
+      CALL histwrite(nid_day,"tpot",itau_w,zx_tmp_2d,
+     .               iim*jjmp1,ndex2d)
+c
+      CALL gr_fi_ecrit(1,klon,iim,jjmp1, tpote,zx_tmp_2d)
+      CALL histwrite(nid_day,"tpote",itau_w,zx_tmp_2d,
+     .               iim*jjmp1,ndex2d)
+c
       ENDIF !lev_histday.GE.5
 c=================================================================