Changeset 653 for LMDZ4/trunk/libf/phylmd

Timestamp:

Jun 13, 2005, 2:58:03 PM (19 years ago)

Author:

Laurent Fairhead

Message:

On avait perdu des modifs de Anne et Marie-Angele
LF

Location:

LMDZ4/trunk/libf/phylmd

Files:

• ini_histmth.h (modified) (3 diffs)
• write_histmth.h (modified) (6 diffs)

LMDZ4/trunk/libf/phylmd/ini_histmth.h

@@ -1392 +1392 @@
      .                "ave(X)", zstophy,zout)
 c
+         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
+     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zstophy,zout)
+c
       ENDIF !iflag_con .GE. 3
 c34EK
@@ -1578 +1582 @@
      .                "ave(X)", zstophy,zout)
 c
-         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
-     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
-     .                "ave(X)", zstophy,zout)
-c
          CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
@@ -1740 +1740 @@
 C
          if (nqmax.GE.3) THEN
-         DO iq=1,nqmax-2
-         IF (iq.LE.99) THEN
-         WRITE(str2,'(i2.2)') iq
-         CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
-     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
-     .                "ave(X)", zstophy,zout)
-         ELSE
-         PRINT*, "Trop de traceurs"
-         CALL abort
-         ENDIF
+         DO iq=3,nqmax
+         iiq=niadv(iq)
+         CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
+     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zstophy,zout)
          ENDDO
          ENDIF
-c
+
       ENDIF !lev_histmth.GE.4
 c

LMDZ4/trunk/libf/phylmd/write_histmth.h

@@ -719 +719 @@
      .                                   iim*jjmp1*klev,ndex3d)
 c
+      IF (iflag_con.GE.3) THEN
+c
       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
       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
       CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
@@ -735 +733 @@
      .                                   iim*jjmp1*klev,ndex3d)
 c
+      ENDIF !iflag_con.GE.3
 c
       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
@@ -881 +880 @@
 C
       IF (nqmax.GE.3) THEN
-      DO iq=1,nqmax-2
-      IF (iq.LE.99) THEN
-         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
-         WRITE(str2,'(i2.2)') iq
-         CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d,
-     .                                   iim*jjmp1*klev,ndex3d)
-      ELSE
-         PRINT*, "Trop de traceurs"
-         CALL abort
-      ENDIF
+      DO iq=3,nqmax
+         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
+         CALL histwrite(nid_mth,tnom(iq),itau_w,zx_tmp_3d,
+     .                                   iim*jjmp1*klev,ndex3d)
       ENDDO
       ENDIF
@@ -1798 +1791 @@
      .                                   iim*jjmp1*klev,ndex3d)
 c
+      IF (iflag_con.GE.3) THEN
+c
       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
       CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
@@ -1814 +1809 @@
      .                                   iim*jjmp1*klev,ndex3d)
 c
+      ENDIF !iflag_con.GE.3
 c
       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
@@ -1945 +1941 @@
 C
       IF (nqmax.GE.3) THEN
-      DO iq=1,nqmax-2
-      IF (iq.LE.99) THEN
-         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
-         WRITE(str2,'(i2.2)') iq
-         CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d,
-     .                                   iim*jjmp1*klev,ndex3d)
-      ELSE
-         PRINT*, "Trop de traceurs"
-         CALL abort
-      ENDIF
+      DO iq=3,nqmax
+         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
+         CALL histwrite(nid_mth,tnom(iq),itau_w,zx_tmp_3d,
+     .                                   iim*jjmp1*klev,ndex3d)
       ENDDO
       ENDIF