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

Timestamp:

Apr 4, 2006, 5:00:40 PM (18 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_histmth.h (modified) (11 diffs)

LMDZ4/trunk/libf/phylmd/write_histmth.h

@@ -207 +207 @@
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
       CALL histwrite(nid_mth,"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_mth,"radsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
 c
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
@@ -355 +358 @@
      $      zx_tmp_2d,iim*jjmp1,ndex2d)
 C
+        zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
+        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
+        CALL histwrite(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
+     $      zx_tmp_2d,iim*jjmp1,ndex2d)
+C
        IF(1.EQ.0) THEN
 C
@@ -517 +525 @@
         zx_tmp_fi2d(i) = fluxo(i)
        ELSE
-        zx_tmp_fi2d(i) = 1.E+20
+        zx_tmp_fi2d(i) = 0.
        ENDIF
       ENDDO
@@ -529 +537 @@
         zx_tmp_fi2d(i) = fluxg(i)
        ELSE
-        zx_tmp_fi2d(i) = 1.E+20
+        zx_tmp_fi2d(i) = 0.
        ENDIF
       ENDDO
@@ -535 +543 @@
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
       CALL histwrite(nid_mth,"fluxg",itau_w,zx_tmp_2d,
+     .     iim*jjmp1,ndex2d)
+c
+      DO i=1, klon
+       IF (pctsrf(i,is_oce).GT.epsfra.OR.
+     .     pctsrf(i,is_sic).GT.epsfra) THEN
+        zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ 
+     .                   ftsol(i, is_sic) * pctsrf(i,is_sic))/
+     .                   (pctsrf(i,is_oce)+pctsrf(i,is_sic))
+       ELSE
+        zx_tmp_fi2d(i) = 273.15
+       ENDIF
+      ENDDO
+c
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
+      CALL histwrite(nid_mth,"t_oce_sic",itau_w,zx_tmp_2d,
      .     iim*jjmp1,ndex2d)
 c
@@ -723 +746 @@
       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
       CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
+     .                                   iim*jjmp1*klev,ndex3d)
+c
+      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
+      CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
      .                                   iim*jjmp1*klev,ndex3d)
 c
@@ -1280 +1307 @@
       CALL histwrite(nid_mth,"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_mth,"radsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
+c
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
       CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
@@ -1426 +1456 @@
      $      zx_tmp_2d,iim*jjmp1,ndex2d)
 C
+        zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
+        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
+        CALL histwrite(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
+     $      zx_tmp_2d,iim*jjmp1,ndex2d)
+C
        IF(1.EQ.0) THEN
 C
@@ -1588 +1623 @@
         zx_tmp_fi2d(i) = fluxo(i)
        ELSE
-        zx_tmp_fi2d(i) = 1.E+20
+        zx_tmp_fi2d(i) = 0.
        ENDIF
       ENDDO
@@ -1600 +1635 @@
         zx_tmp_fi2d(i) = fluxg(i)
        ELSE
-        zx_tmp_fi2d(i) = 1.E+20
+        zx_tmp_fi2d(i) = 0.
        ENDIF
       ENDDO
@@ -1606 +1641 @@
       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
       CALL histwrite(nid_mth,"fluxg",itau_w,zx_tmp_2d,
+     .     iim*jjmp1,ndex2d)
+c
+      DO i=1, klon
+       IF (pctsrf(i,is_oce).GT.epsfra.OR.
+     .     pctsrf(i,is_sic).GT.epsfra) THEN
+        zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ 
+     .                   ftsol(i, is_sic) * pctsrf(i,is_sic))/
+     .                   (pctsrf(i,is_oce)+pctsrf(i,is_sic))
+       ELSE
+        zx_tmp_fi2d(i) = 273.15
+       ENDIF
+      ENDDO
+c
+      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
+      CALL histwrite(nid_mth,"t_oce_sic",itau_w,zx_tmp_2d,
      .     iim*jjmp1,ndex2d)
 c