Changeset 653 for LMDZ4/trunk


Ignore:
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:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk/libf/phylmd/ini_histmth.h

    r644 r653  
    13921392     .                "ave(X)", zstophy,zout)
    13931393c
     1394         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
     1395     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1396     .                "ave(X)", zstophy,zout)
     1397c
    13941398      ENDIF !iflag_con .GE. 3
    13951399c34EK
    … …  
    15781582     .                "ave(X)", zstophy,zout)
    15791583c
    1580          CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
    1581      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    1582      .                "ave(X)", zstophy,zout)
    1583 c
    15841584         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
    15851585     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    … …  
    17401740C
    17411741         if (nqmax.GE.3) THEN
    1742          DO iq=1,nqmax-2
    1743          IF (iq.LE.99) THEN
    1744          WRITE(str2,'(i2.2)') iq
    1745          CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
    1746      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    1747      .                "ave(X)", zstophy,zout)
    1748          ELSE
    1749          PRINT*, "Trop de traceurs"
    1750          CALL abort
    1751          ENDIF
     1742         DO iq=3,nqmax
     1743         iiq=niadv(iq)
     1744         CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
     1745     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1746     .                "ave(X)", zstophy,zout)
    17521747         ENDDO
    17531748         ENDIF
    1754 c
     1749
    17551750      ENDIF !lev_histmth.GE.4
    17561751c

  • LMDZ4/trunk/libf/phylmd/write_histmth.h

    r644 r653  
    719719     .                                   iim*jjmp1*klev,ndex3d)
    720720c
     721      IF (iflag_con.GE.3) THEN
     722c
    721723      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
    722724      CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
    723725     .                                   iim*jjmp1*klev,ndex3d)
    724726c
    725       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
    726       CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
    727      .                                   iim*jjmp1*klev,ndex3d)
    728 c
    729727      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
    730728      CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
    … …  
    735733     .                                   iim*jjmp1*klev,ndex3d)
    736734c
     735      ENDIF !iflag_con.GE.3
    737736c
    738737      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
    … …  
    881880C
    882881      IF (nqmax.GE.3) THEN
    883       DO iq=1,nqmax-2
    884       IF (iq.LE.99) THEN
    885          CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
    886          WRITE(str2,'(i2.2)') iq
    887          CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d,
    888      .                                   iim*jjmp1*klev,ndex3d)
    889       ELSE
    890          PRINT*, "Trop de traceurs"
    891          CALL abort
    892       ENDIF
     882      DO iq=3,nqmax
     883         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
     884         CALL histwrite(nid_mth,tnom(iq),itau_w,zx_tmp_3d,
     885     .                                   iim*jjmp1*klev,ndex3d)
    893886      ENDDO
    894887      ENDIF
    … …  
    17981791     .                                   iim*jjmp1*klev,ndex3d)
    17991792c
     1793      IF (iflag_con.GE.3) THEN
     1794c
    18001795      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
    18011796      CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
    … …  
    18141809     .                                   iim*jjmp1*klev,ndex3d)
    18151810c
     1811      ENDIF !iflag_con.GE.3
    18161812c
    18171813      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
    … …  
    19451941C
    19461942      IF (nqmax.GE.3) THEN
    1947       DO iq=1,nqmax-2
    1948       IF (iq.LE.99) THEN
    1949          CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
    1950          WRITE(str2,'(i2.2)') iq
    1951          CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d,
    1952      .                                   iim*jjmp1*klev,ndex3d)
    1953       ELSE
    1954          PRINT*, "Trop de traceurs"
    1955          CALL abort
    1956       ENDIF
     1943      DO iq=3,nqmax
     1944         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
     1945         CALL histwrite(nid_mth,tnom(iq),itau_w,zx_tmp_3d,
     1946     .                                   iim*jjmp1*klev,ndex3d)
    19571947      ENDDO
    19581948      ENDIF
Note: See TracChangeset for help on using the changeset viewer.