Changeset 29 for LMDZ.3.3/trunk/libf


Ignore:
Timestamp:
Jan 24, 2000, 1:07:41 PM (25 years ago)
Author:
lmdz
Message:

Le tableau ndex, qui permet de "flagger" les valeurs a sortir dans IOIPSL, est
dimmensionne maniere correcte, sinon sur Nec les champs 3D physique ne sont
pas correctement rempli. (probleme non visible avec les champs dynamiques mais
pour faire propre ...) LF

Location:
LMDZ.3.3/trunk/libf
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/trunk/libf/bibio/writedynav.F

    r2 r29  
    5858C   Variables locales
    5959C
    60       integer ndex(1),iq, ii, ll
     60      integer ndex2d(iip1*jjp1),ndex3d(iip1*jjp1*llm),iq, ii, ll
    6161      character*3 str
    6262      real us(ip1jmp1*llm), vs(ip1jmp1*llm)
    … …  
    6868C
    6969      str='q  '
    70       ndex(1) = 0
     70      ndex3d = 0
     71      ndex2d = 0
    7172      ok_sync = .TRUE.
    7273
    … …  
    8182      call gr_u_scal(llm, unat, us)
    8283      call histwrite(histid, 'u', time, us,
    83      .               iip1*jjp1*llm, ndex)
     84     .               iip1*jjp1*llm, ndex3d)
    8485C
    8586C  Vents V scalaire
    … …  
    8788      call gr_v_scal(llm, vnat, vs)
    8889      call histwrite(histid, 'v', time, vs,
    89      .               iip1*jjp1*llm, ndex)
     90     .               iip1*jjp1*llm, ndex3d)
    9091C
    9192C  Temperature moyennee
    … …  
    9596      enddo
    9697      call histwrite(histid, 'temp', time, tm,
    97      .                iip1*jjp1*llm, ndex)
     98     .                iip1*jjp1*llm, ndex3d)
    9899C
    99100C  Geopotentiel
    100101C
    101102      call histwrite(histid, 'phi', time, phi,
    102      .                iip1*jjp1*llm, ndex)
     103     .                iip1*jjp1*llm, ndex3d)
    103104C
    104105C  Traceurs
    … …  
    112113          ENDIF
    113114          call histwrite(histid, str, time, q(:,:,iq),
    114      .                   iip1*jjp1*llm, ndex)
     115     .                   iip1*jjp1*llm, ndex3d)
    115116        enddo
    116117      endif
    … …  
    118119C  Masse
    119120C
    120        call histwrite(histid, 'masse', time, masse, iip1*jjp1, ndex)
     121       call histwrite(histid, 'masse', time, masse, iip1*jjp1, ndex2d)
    121122C
    122123C  Pression au sol
    123124C
    124        call histwrite(histid, 'ps', time, ps, iip1*jjp1, ndex)
     125       call histwrite(histid, 'ps', time, ps, iip1*jjp1, ndex2d)
    125126C
    126127C  Geopotentiel au sol
    127128C
    128        call histwrite(histid, 'phis', time, phis, iip1*jjp1, ndex)
     129       call histwrite(histid, 'phis', time, phis, iip1*jjp1, ndex2d)
    129130C
    130131C  Fin
    131132C
    132       if (ok_sync) call histsync
     133      if (ok_sync) call histsync(histid)
    133134      return
    134135      end

  • LMDZ.3.3/trunk/libf/bibio/writehist.F

    r2 r29  
    5959C   Variables locales
    6060C
    61       integer ndex(1),iq, ii, ll
     61      integer iq, ii, ll
     62      integer ndexu(ip1jmp1,llm),ndexv(ip1jm,llm),ndex2d(ip1jmp1)
    6263      character*3 str
    6364      logical ok_sync
    … …  
    6667C
    6768      str='q  '
    68       ndex(1) = 0
     69      ndexu = 0
     70      ndexv = 0
     71      ndex2d = 0
    6972      ok_sync =.TRUE.
    7073C
    … …  
    7477C
    7578      call histwrite(histid, 'ucov', time, ucov,
    76      .               iip1*jjp1*llm, ndex)
     79     .               iip1*jjp1*llm, ndexu)
    7780
    7881C
    … …  
    8083C
    8184      call histwrite(histvid, 'vcov', time, vcov,
    82      .               iip1*jjm*llm, ndex)
     85     .               iip1*jjm*llm, ndexv)
    8386
    8487C
    … …  
    8689C
    8790      call histwrite(histid, 'teta', time, teta,
    88      .                iip1*jjp1*llm, ndex)
     91     .                iip1*jjp1*llm, ndexu)
    8992C
    9093C  Geopotentiel
    9194C
    9295      call histwrite(histid, 'phi', time, phi,
    93      .                iip1*jjp1*llm, ndex)
     96     .                iip1*jjp1*llm, ndexu)
    9497C
    9598C  Traceurs
    … …  
    103106          ENDIF
    104107          call histwrite(histid, str, time, q(:,:,iq),
    105      .                   iip1*jjp1*llm, ndex)
     108     .                   iip1*jjp1*llm, ndexu)
    106109        enddo
    107110      endif
    … …  
    109112C  Masse
    110113C
    111       call histwrite(histid, 'masse', time, masse, iip1*jjp1, ndex)
     114      call histwrite(histid, 'masse', time, masse, iip1*jjp1, ndex2d)
    112115C
    113116C  Pression au sol
    114117C
    115       call histwrite(histid, 'ps', time, ps, iip1*jjp1, ndex)
     118      call histwrite(histid, 'ps', time, ps, iip1*jjp1, ndex2d)
    116119C
    117120C  Geopotentiel au sol
    118121C
    119       call histwrite(histid, 'phis', time, phis, iip1*jjp1, ndex)
     122      call histwrite(histid, 'phis', time, phis, iip1*jjp1, ndex2d)
    120123C
    121124C  Fin
    122125C
    123       if (ok_sync) call histsync
     126      if (ok_sync) then
     127        call histsync(histid)
     128        call histsync(histvid)
     129      endif
    124130      return
    125131      end

  • LMDZ.3.3/trunk/libf/phylmd/physiq.F

    r23 r29  
    506506      REAL tabcntr0( length       )
    507507c
    508       INTEGER ndex(1)
     508      INTEGER ndex2d(iim*(jjm+1)),ndex3d(iim*(jjm+1)*klev)
    509509      REAL zx_tmp_fi2d(klon)
    510510      REAL zx_tmp_2d(iim,jjm+1), zx_tmp_3d(iim,jjm+1,klev)
    … …  
    810810         CALL histend(nid_day)
    811811c
    812          ndex(1) = 0
    813 c
    814          i = NINT(zout/zsto)
    815          CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
    816          CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex)
    817 c
    818          i = NINT(zout/zsto)
    819          CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
    820          CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex)
     812         ndex2d = 0
     813         ndex3d = 0
     814c
     815Cess         i = NINT(zout/zsto)
     816Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
     817Cess        CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
     818c
     819Cess         i = NINT(zout/zsto)
     820Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
     821Cess        CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
    821822c
    822823      ENDIF ! fin de test sur ok_journe
    … …  
    11571158         CALL histend(nid_mth)
    11581159c
    1159          ndex(1) = 0
    1160 c
    1161          i = NINT(zout/zsto)
    1162          CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
    1163          CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex)
     1160         ndex2d = 0
     1161         ndex3d = 0
     1162c
     1163Cess         i = NINT(zout/zsto)
     1164Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
     1165Cess        CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
    11641166C
    1165          i = NINT(zout/zsto)
    1166          CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
    1167          CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex)
     1167Cess         i = NINT(zout/zsto)
     1168Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
     1169Cess        CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
    11681170c
    11691171      ENDIF ! fin de test sur ok_mensuel
    … …  
    12381240         CALL histend(nid_ins)
    12391241c
    1240          ndex(1) = 0
    1241 c
    1242          i = NINT(zout/zsto)
    1243          CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
    1244          CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex)
    1245 c
    1246          i = NINT(zout/zsto)
    1247          CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
    1248          CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex)
     1242         ndex2d = 0
     1243         ndex3d = 0
     1244c
     1245Cess         i = NINT(zout/zsto)
     1246Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
     1247Cess        CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
     1248c
     1249Cess         i = NINT(zout/zsto)
     1250Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
     1251Cess        CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
    12491252c
    12501253      ENDIF
    … …  
    17971800C la physique s'applique
    17981801C
     1802       write(*,*) 'Phytrac= '
    17991803      call phytrac (rnpb,
    18001804     I                   debut,
    … …  
    18081812     O                   tr_seri)
    18091813
     1814       write(*,*) 'OFFLINE= ', offline
    18101815      IF (offline) THEN
     1816       write(*,*) 'OFFLINE= ', offline
    18111817      call phystoke (
    18121818     I                   nlon,nlev,pdtphys,
    … …  
    18451851      IF (ok_journe) THEN
    18461852c
    1847       ndex(1) = 0
     1853      ndex2d = 0
     1854      ndex3d = 0
    18481855c
    18491856c Champs 2D:
    18501857c
     1858         i = NINT(zout/zsto)
     1859         CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
     1860         CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
     1861c
     1862         i = NINT(zout/zsto)
     1863         CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
     1864         CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
     1865C
    18511866      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zxtsol,zx_tmp_2d)
    1852       CALL histwrite(nid_day,"tsol",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1867      CALL histwrite(nid_day,"tsol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18531868c
    18541869      DO i = 1, klon
    … …  
    18561871      ENDDO
    18571872      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    1858       CALL histwrite(nid_day,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1873      CALL histwrite(nid_day,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18591874c
    18601875      DO i = 1, klon
    … …  
    18621877      ENDDO
    18631878      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    1864       CALL histwrite(nid_day,"rain",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1879      CALL histwrite(nid_day,"rain",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18651880c
    18661881      CALL gr_fi_ecrit(1, klon,iim,jjm+1, snow_fall,zx_tmp_2d)
    1867       CALL histwrite(nid_day,"snow",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1882      CALL histwrite(nid_day,"snow",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18681883c
    18691884      CALL gr_fi_ecrit(1, klon,iim,jjm+1, evap,zx_tmp_2d)
    1870       CALL histwrite(nid_day,"evap",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1885      CALL histwrite(nid_day,"evap",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18711886c
    18721887      CALL gr_fi_ecrit(1, klon,iim,jjm+1, topsw,zx_tmp_2d)
    1873       CALL histwrite(nid_day,"tops",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1888      CALL histwrite(nid_day,"tops",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18741889c
    18751890      CALL gr_fi_ecrit(1, klon,iim,jjm+1, toplw,zx_tmp_2d)
    1876       CALL histwrite(nid_day,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1891      CALL histwrite(nid_day,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18771892c
    18781893      CALL gr_fi_ecrit(1, klon,iim,jjm+1, solsw,zx_tmp_2d)
    1879       CALL histwrite(nid_day,"sols",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1894      CALL histwrite(nid_day,"sols",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18801895c
    18811896      CALL gr_fi_ecrit(1, klon,iim,jjm+1, sollw,zx_tmp_2d)
    1882       CALL histwrite(nid_day,"soll",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1897      CALL histwrite(nid_day,"soll",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18831898c
    18841899      CALL gr_fi_ecrit(1, klon,iim,jjm+1, bils,zx_tmp_2d)
    1885       CALL histwrite(nid_day,"bils",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1900      CALL histwrite(nid_day,"bils",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18861901c
    18871902      CALL gr_fi_ecrit(1, klon,iim,jjm+1, sens,zx_tmp_2d)
    1888       CALL histwrite(nid_day,"sens",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1903      CALL histwrite(nid_day,"sens",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18891904c
    18901905      CALL gr_fi_ecrit(1, klon,iim,jjm+1, fder,zx_tmp_2d)
    1891       CALL histwrite(nid_day,"fder",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1906      CALL histwrite(nid_day,"fder",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18921907c
    18931908      CALL gr_fi_ecrit(1, klon,iim,jjm+1, ruis,zx_tmp_2d)
    1894       CALL histwrite(nid_day,"ruis",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1909      CALL histwrite(nid_day,"ruis",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    18951910c
    18961911      DO i = 1, klon
    … …  
    18981913      ENDDO
    18991914      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    1900       CALL histwrite(nid_day,"frtu",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1915      CALL histwrite(nid_day,"frtu",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19011916c
    19021917      DO i = 1, klon
    … …  
    19041919      ENDDO
    19051920      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    1906       CALL histwrite(nid_day,"frtv",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1921      CALL histwrite(nid_day,"frtv",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19071922c
    19081923      DO i = 1, klon
    … …  
    19101925      ENDDO
    19111926      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    1912       CALL histwrite(nid_day,"sicf",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1927      CALL histwrite(nid_day,"sicf",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19131928c
    19141929      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldl,zx_tmp_2d)
    1915       CALL histwrite(nid_day,"cldl",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1930      CALL histwrite(nid_day,"cldl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19161931c
    19171932      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldm,zx_tmp_2d)
    1918       CALL histwrite(nid_day,"cldm",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1933      CALL histwrite(nid_day,"cldm",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19191934c
    19201935      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldh,zx_tmp_2d)
    1921       CALL histwrite(nid_day,"cldh",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1936      CALL histwrite(nid_day,"cldh",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19221937c
    19231938      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldt,zx_tmp_2d)
    1924       CALL histwrite(nid_day,"cldt",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1939      CALL histwrite(nid_day,"cldt",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19251940c
    19261941      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldq,zx_tmp_2d)
    1927       CALL histwrite(nid_day,"cldq",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1942      CALL histwrite(nid_day,"cldq",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19281943c
    19291944c Champs 3D:
    … …  
    19311946      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t_seri, zx_tmp_3d)
    19321947      CALL histwrite(nid_day,"temp",itap,zx_tmp_3d,
    1933      .                                   iim*(jjm+1)*klev,ndex)
     1948     .                                   iim*(jjm+1)*klev,ndex3d)
    19341949c
    19351950      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, qx(1,1,ivap), zx_tmp_3d)
    19361951      CALL histwrite(nid_day,"ovap",itap,zx_tmp_3d,
    1937      .                                   iim*(jjm+1)*klev,ndex)
     1952     .                                   iim*(jjm+1)*klev,ndex3d)
    19381953c
    19391954      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, zphi, zx_tmp_3d)
    19401955      CALL histwrite(nid_day,"geop",itap,zx_tmp_3d,
    1941      .                                   iim*(jjm+1)*klev,ndex)
     1956     .                                   iim*(jjm+1)*klev,ndex3d)
    19421957c
    19431958      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, u_seri, zx_tmp_3d)
    19441959      CALL histwrite(nid_day,"vitu",itap,zx_tmp_3d,
    1945      .                                   iim*(jjm+1)*klev,ndex)
     1960     .                                   iim*(jjm+1)*klev,ndex3d)
    19461961c
    19471962      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, v_seri, zx_tmp_3d)
    19481963      CALL histwrite(nid_day,"vitv",itap,zx_tmp_3d,
    1949      .                                   iim*(jjm+1)*klev,ndex)
     1964     .                                   iim*(jjm+1)*klev,ndex3d)
    19501965c
    19511966      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, omega, zx_tmp_3d)
    19521967      CALL histwrite(nid_day,"vitw",itap,zx_tmp_3d,
    1953      .                                   iim*(jjm+1)*klev,ndex)
     1968     .                                   iim*(jjm+1)*klev,ndex3d)
    19541969c
    19551970      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, pplay, zx_tmp_3d)
    19561971      CALL histwrite(nid_day,"pres",itap,zx_tmp_3d,
    1957      .                                   iim*(jjm+1)*klev,ndex)
    1958 c
    1959       if (ok_sync) call histsync
     1972     .                                   iim*(jjm+1)*klev,ndex3d)
     1973c
     1974      if (ok_sync) then
     1975        call histsync(nid_day)
     1976      endif
    19601977      ENDIF
    19611978C
    19621979      IF (ok_mensuel) THEN
    19631980c
    1964       ndex(1) = 0
     1981      ndex2d = 0
     1982      ndex3d = 0
    19651983c
    19661984c Champs 2D:
    19671985c
     1986         i = NINT(zout/zsto)
     1987         CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
     1988         CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
     1989C
     1990         i = NINT(zout/zsto)
     1991         CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
     1992         CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
     1993
    19681994      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zxtsol,zx_tmp_2d)
    1969       CALL histwrite(nid_mth,"tsol",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     1995      CALL histwrite(nid_mth,"tsol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19701996c
    19711997      DO i = 1, klon
    … …  
    19731999      ENDDO
    19742000      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    1975       CALL histwrite(nid_mth,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2001      CALL histwrite(nid_mth,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19762002c
    19772003      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zxqsol,zx_tmp_2d)
    1978       CALL histwrite(nid_mth,"qsol",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2004      CALL histwrite(nid_mth,"qsol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19792005c
    19802006      DO i = 1, klon
    … …  
    19822008      ENDDO
    19832009      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    1984       CALL histwrite(nid_mth,"rain",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2010      CALL histwrite(nid_mth,"rain",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19852011c
    19862012      DO i = 1, klon
    … …  
    19882014      ENDDO
    19892015      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    1990       CALL histwrite(nid_mth,"plul",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2016      CALL histwrite(nid_mth,"plul",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19912017c
    19922018      DO i = 1, klon
    … …  
    19942020      ENDDO
    19952021      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    1996       CALL histwrite(nid_mth,"pluc",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2022      CALL histwrite(nid_mth,"pluc",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    19972023c
    19982024      CALL gr_fi_ecrit(1, klon,iim,jjm+1, snow_fall,zx_tmp_2d)
    1999       CALL histwrite(nid_mth,"snow",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2025      CALL histwrite(nid_mth,"snow",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20002026c
    20012027      CALL gr_fi_ecrit(1, klon,iim,jjm+1, agesno,zx_tmp_2d)
    2002       CALL histwrite(nid_mth,"ages",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2028      CALL histwrite(nid_mth,"ages",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20032029c
    20042030      CALL gr_fi_ecrit(1, klon,iim,jjm+1, evap,zx_tmp_2d)
    2005       CALL histwrite(nid_mth,"evap",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2031      CALL histwrite(nid_mth,"evap",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20062032c
    20072033      CALL gr_fi_ecrit(1, klon,iim,jjm+1, topsw,zx_tmp_2d)
    2008       CALL histwrite(nid_mth,"tops",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2034      CALL histwrite(nid_mth,"tops",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20092035c
    20102036      CALL gr_fi_ecrit(1, klon,iim,jjm+1, toplw,zx_tmp_2d)
    2011       CALL histwrite(nid_mth,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2037      CALL histwrite(nid_mth,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20122038c
    20132039      CALL gr_fi_ecrit(1, klon,iim,jjm+1, solsw,zx_tmp_2d)
    2014       CALL histwrite(nid_mth,"sols",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2040      CALL histwrite(nid_mth,"sols",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20152041c
    20162042      CALL gr_fi_ecrit(1, klon,iim,jjm+1, sollw,zx_tmp_2d)
    2017       CALL histwrite(nid_mth,"soll",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2043      CALL histwrite(nid_mth,"soll",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20182044c
    20192045      CALL gr_fi_ecrit(1, klon,iim,jjm+1, topsw0,zx_tmp_2d)
    2020       CALL histwrite(nid_mth,"tops0",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2046      CALL histwrite(nid_mth,"tops0",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20212047c
    20222048      CALL gr_fi_ecrit(1, klon,iim,jjm+1, toplw0,zx_tmp_2d)
    2023       CALL histwrite(nid_mth,"topl0",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2049      CALL histwrite(nid_mth,"topl0",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20242050c
    20252051      CALL gr_fi_ecrit(1, klon,iim,jjm+1, solsw0,zx_tmp_2d)
    2026       CALL histwrite(nid_mth,"sols0",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2052      CALL histwrite(nid_mth,"sols0",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20272053c
    20282054      CALL gr_fi_ecrit(1, klon,iim,jjm+1, sollw0,zx_tmp_2d)
    2029       CALL histwrite(nid_mth,"soll0",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2055      CALL histwrite(nid_mth,"soll0",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20302056c
    20312057      CALL gr_fi_ecrit(1, klon,iim,jjm+1, bils,zx_tmp_2d)
    2032       CALL histwrite(nid_mth,"bils",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2058      CALL histwrite(nid_mth,"bils",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20332059c
    20342060      CALL gr_fi_ecrit(1, klon,iim,jjm+1, sens,zx_tmp_2d)
    2035       CALL histwrite(nid_mth,"sens",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2061      CALL histwrite(nid_mth,"sens",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20362062c
    20372063      CALL gr_fi_ecrit(1, klon,iim,jjm+1, fder,zx_tmp_2d)
    2038       CALL histwrite(nid_mth,"fder",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2064      CALL histwrite(nid_mth,"fder",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20392065c
    20402066      CALL gr_fi_ecrit(1, klon,iim,jjm+1, ruis,zx_tmp_2d)
    2041       CALL histwrite(nid_mth,"ruis",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2067      CALL histwrite(nid_mth,"ruis",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20422068c
    20432069      DO i = 1, klon
    … …  
    20452071      ENDDO
    20462072      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    2047       CALL histwrite(nid_mth,"frtu",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2073      CALL histwrite(nid_mth,"frtu",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20482074c
    20492075      DO i = 1, klon
    … …  
    20512077      ENDDO
    20522078      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    2053       CALL histwrite(nid_mth,"frtv",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2079      CALL histwrite(nid_mth,"frtv",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20542080c
    20552081      DO i = 1, klon
    … …  
    20572083      ENDDO
    20582084      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    2059       CALL histwrite(nid_mth,"sicf",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2085      CALL histwrite(nid_mth,"sicf",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20602086c
    20612087      CALL gr_fi_ecrit(1, klon,iim,jjm+1, albsol,zx_tmp_2d)
    2062       CALL histwrite(nid_mth,"albs",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2088      CALL histwrite(nid_mth,"albs",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20632089c
    20642090      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cdragm,zx_tmp_2d)
    2065       CALL histwrite(nid_mth,"cdrm",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2091      CALL histwrite(nid_mth,"cdrm",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20662092c
    20672093      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cdragh,zx_tmp_2d)
    2068       CALL histwrite(nid_mth,"cdrh",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2094      CALL histwrite(nid_mth,"cdrh",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20692095c
    20702096      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldl,zx_tmp_2d)
    2071       CALL histwrite(nid_mth,"cldl",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2097      CALL histwrite(nid_mth,"cldl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20722098c
    20732099      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldm,zx_tmp_2d)
    2074       CALL histwrite(nid_mth,"cldm",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2100      CALL histwrite(nid_mth,"cldm",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20752101c
    20762102      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldh,zx_tmp_2d)
    2077       CALL histwrite(nid_mth,"cldh",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2103      CALL histwrite(nid_mth,"cldh",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20782104c
    20792105      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldt,zx_tmp_2d)
    2080       CALL histwrite(nid_mth,"cldt",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2106      CALL histwrite(nid_mth,"cldt",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20812107c
    20822108      CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldq,zx_tmp_2d)
    2083       CALL histwrite(nid_mth,"cldq",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2109      CALL histwrite(nid_mth,"cldq",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20842110c
    20852111      CALL gr_fi_ecrit(1, klon,iim,jjm+1, ue,zx_tmp_2d)
    2086       CALL histwrite(nid_mth,"ue",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2112      CALL histwrite(nid_mth,"ue",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20872113c
    20882114      CALL gr_fi_ecrit(1, klon,iim,jjm+1, ve,zx_tmp_2d)
    2089       CALL histwrite(nid_mth,"ve",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2115      CALL histwrite(nid_mth,"ve",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20902116c
    20912117      CALL gr_fi_ecrit(1, klon,iim,jjm+1, uq,zx_tmp_2d)
    2092       CALL histwrite(nid_mth,"uq",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2118      CALL histwrite(nid_mth,"uq",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20932119c
    20942120      CALL gr_fi_ecrit(1, klon,iim,jjm+1, vq,zx_tmp_2d)
    2095       CALL histwrite(nid_mth,"vq",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2121      CALL histwrite(nid_mth,"vq",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    20962122c
    20972123c Champs 3D:
    … …  
    20992125      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t_seri, zx_tmp_3d)
    21002126      CALL histwrite(nid_mth,"temp",itap,zx_tmp_3d,
    2101      .                                   iim*(jjm+1)*klev,ndex)
     2127     .                                   iim*(jjm+1)*klev,ndex3d)
    21022128c
    21032129      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, qx(1,1,ivap), zx_tmp_3d)
    21042130      CALL histwrite(nid_mth,"ovap",itap,zx_tmp_3d,
    2105      .                                   iim*(jjm+1)*klev,ndex)
     2131     .                                   iim*(jjm+1)*klev,ndex3d)
    21062132c
    21072133      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, zphi, zx_tmp_3d)
    21082134      CALL histwrite(nid_mth,"geop",itap,zx_tmp_3d,
    2109      .                                   iim*(jjm+1)*klev,ndex)
     2135     .                                   iim*(jjm+1)*klev,ndex3d)
    21102136c
    21112137      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, u_seri, zx_tmp_3d)
    21122138      CALL histwrite(nid_mth,"vitu",itap,zx_tmp_3d,
    2113      .                                   iim*(jjm+1)*klev,ndex)
     2139     .                                   iim*(jjm+1)*klev,ndex3d)
    21142140c
    21152141      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, v_seri, zx_tmp_3d)
    21162142      CALL histwrite(nid_mth,"vitv",itap,zx_tmp_3d,
    2117      .                                   iim*(jjm+1)*klev,ndex)
     2143     .                                   iim*(jjm+1)*klev,ndex3d)
    21182144c
    21192145      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, omega, zx_tmp_3d)
    21202146      CALL histwrite(nid_mth,"vitw",itap,zx_tmp_3d,
    2121      .                                   iim*(jjm+1)*klev,ndex)
     2147     .                                   iim*(jjm+1)*klev,ndex3d)
    21222148c
    21232149      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, pplay, zx_tmp_3d)
    21242150      CALL histwrite(nid_mth,"pres",itap,zx_tmp_3d,
    2125      .                                   iim*(jjm+1)*klev,ndex)
     2151     .                                   iim*(jjm+1)*klev,ndex3d)
    21262152c
    21272153      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, cldfra, zx_tmp_3d)
    21282154      CALL histwrite(nid_mth,"rneb",itap,zx_tmp_3d,
    2129      .                                   iim*(jjm+1)*klev,ndex)
     2155     .                                   iim*(jjm+1)*klev,ndex3d)
    21302156c
    21312157      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, zx_rh, zx_tmp_3d)
    21322158      CALL histwrite(nid_mth,"rhum",itap,zx_tmp_3d,
    2133      .                                   iim*(jjm+1)*klev,ndex)
     2159     .                                   iim*(jjm+1)*klev,ndex3d)
    21342160c
    21352161      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, cldliq, zx_tmp_3d)
    21362162      CALL histwrite(nid_mth,"oliq",itap,zx_tmp_3d,
    2137      .                                   iim*(jjm+1)*klev,ndex)
     2163     .                                   iim*(jjm+1)*klev,ndex3d)
    21382164c
    21392165      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_dyn, zx_tmp_3d)
    21402166      CALL histwrite(nid_mth,"dtdyn",itap,zx_tmp_3d,
    2141      .                                   iim*(jjm+1)*klev,ndex)
     2167     .                                   iim*(jjm+1)*klev,ndex3d)
    21422168c
    21432169      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_qx_dyn(1,1,ivap),
    21442170     .                 zx_tmp_3d)
    21452171      CALL histwrite(nid_mth,"dqdyn",itap,zx_tmp_3d,
    2146      .                                   iim*(jjm+1)*klev,ndex)
     2172     .                                   iim*(jjm+1)*klev,ndex3d)
    21472173c
    21482174      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_con, zx_tmp_3d)
    21492175      CALL histwrite(nid_mth,"dtcon",itap,zx_tmp_3d,
    2150      .                                   iim*(jjm+1)*klev,ndex)
     2176     .                                   iim*(jjm+1)*klev,ndex3d)
    21512177c
    21522178      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_con, zx_tmp_3d)
    21532179      CALL histwrite(nid_mth,"dqcon",itap,zx_tmp_3d,
    2154      .                                   iim*(jjm+1)*klev,ndex)
     2180     .                                   iim*(jjm+1)*klev,ndex3d)
    21552181c
    21562182      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_lsc, zx_tmp_3d)
    21572183      CALL histwrite(nid_mth,"dtlsc",itap,zx_tmp_3d,
    2158      .                                   iim*(jjm+1)*klev,ndex)
     2184     .                                   iim*(jjm+1)*klev,ndex3d)
    21592185c
    21602186      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_lsc, zx_tmp_3d)
    21612187      CALL histwrite(nid_mth,"dqlsc",itap,zx_tmp_3d,
    2162      .                                   iim*(jjm+1)*klev,ndex)
     2188     .                                   iim*(jjm+1)*klev,ndex3d)
    21632189c
    21642190      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_vdf, zx_tmp_3d)
    21652191      CALL histwrite(nid_mth,"dtvdf",itap,zx_tmp_3d,
    2166      .                                   iim*(jjm+1)*klev,ndex)
     2192     .                                   iim*(jjm+1)*klev,ndex3d)
    21672193c
    21682194      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_vdf, zx_tmp_3d)
    21692195      CALL histwrite(nid_mth,"dqvdf",itap,zx_tmp_3d,
    2170      .                                   iim*(jjm+1)*klev,ndex)
     2196     .                                   iim*(jjm+1)*klev,ndex3d)
    21712197c
    21722198      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_eva, zx_tmp_3d)
    21732199      CALL histwrite(nid_mth,"dteva",itap,zx_tmp_3d,
    2174      .                                   iim*(jjm+1)*klev,ndex)
     2200     .                                   iim*(jjm+1)*klev,ndex3d)
    21752201c
    21762202      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_eva, zx_tmp_3d)
    21772203      CALL histwrite(nid_mth,"dqeva",itap,zx_tmp_3d,
    2178      .                                   iim*(jjm+1)*klev,ndex)
     2204     .                                   iim*(jjm+1)*klev,ndex3d)
    21792205c
    21802206C      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_ajs, zx_tmp_3d)
    21812207C      CALL histwrite(nid_mth,"dtajs",itap,zx_tmp_3d,
    2182 C    .                                   iim*(jjm+1)*klev,ndex)
     2208C    .                                   iim*(jjm+1)*klev,ndex3d)
    21832209c
    21842210C      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_ajs, zx_tmp_3d)
    21852211C      CALL histwrite(nid_mth,"dqajs",itap,zx_tmp_3d,
    2186 C     .                                   iim*(jjm+1)*klev,ndex)
     2212C     .                                   iim*(jjm+1)*klev,ndex3d)
    21872213c
    21882214      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, heat, zx_tmp_3d)
    21892215      CALL histwrite(nid_mth,"dtswr",itap,zx_tmp_3d,
    2190      .                                   iim*(jjm+1)*klev,ndex)
     2216     .                                   iim*(jjm+1)*klev,ndex3d)
    21912217c
    21922218      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, heat0, zx_tmp_3d)
    21932219      CALL histwrite(nid_mth,"dtsw0",itap,zx_tmp_3d,
    2194      .                                   iim*(jjm+1)*klev,ndex)
     2220     .                                   iim*(jjm+1)*klev,ndex3d)
    21952221c
    21962222      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, cool, zx_tmp_3d)
    21972223      CALL histwrite(nid_mth,"dtlwr",itap,zx_tmp_3d,
    2198      .                                   iim*(jjm+1)*klev,ndex)
     2224     .                                   iim*(jjm+1)*klev,ndex3d)
    21992225c
    22002226      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, cool0, zx_tmp_3d)
    22012227      CALL histwrite(nid_mth,"dtlw0",itap,zx_tmp_3d,
    2202      .                                   iim*(jjm+1)*klev,ndex)
     2228     .                                   iim*(jjm+1)*klev,ndex3d)
    22032229c
    22042230      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_u_vdf, zx_tmp_3d)
    22052231      CALL histwrite(nid_mth,"duvdf",itap,zx_tmp_3d,
    2206      .                                   iim*(jjm+1)*klev,ndex)
     2232     .                                   iim*(jjm+1)*klev,ndex3d)
    22072233c
    22082234      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_v_vdf, zx_tmp_3d)
    22092235      CALL histwrite(nid_mth,"dvvdf",itap,zx_tmp_3d,
    2210      .                                   iim*(jjm+1)*klev,ndex)
     2236     .                                   iim*(jjm+1)*klev,ndex3d)
    22112237c
    22122238      IF (ok_orodr) THEN
    22132239      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_u_oro, zx_tmp_3d)
    22142240      CALL histwrite(nid_mth,"duoro",itap,zx_tmp_3d,
    2215      .                                   iim*(jjm+1)*klev,ndex)
     2241     .                                   iim*(jjm+1)*klev,ndex3d)
    22162242c
    22172243      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_v_oro, zx_tmp_3d)
    22182244      CALL histwrite(nid_mth,"dvoro",itap,zx_tmp_3d,
    2219      .                                   iim*(jjm+1)*klev,ndex)
     2245     .                                   iim*(jjm+1)*klev,ndex3d)
    22202246c
    22212247      ENDIF
    … …  
    22242250      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_u_lif, zx_tmp_3d)
    22252251      CALL histwrite(nid_mth,"dulif",itap,zx_tmp_3d,
    2226      .                                   iim*(jjm+1)*klev,ndex)
     2252     .                                   iim*(jjm+1)*klev,ndex3d)
    22272253c
    22282254      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_v_lif, zx_tmp_3d)
    22292255      CALL histwrite(nid_mth,"dvlif",itap,zx_tmp_3d,
    2230      .                                   iim*(jjm+1)*klev,ndex)
     2256     .                                   iim*(jjm+1)*klev,ndex3d)
    22312257      ENDIF
    22322258C
    22332259      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wo, zx_tmp_3d)
    22342260      CALL histwrite(nid_mth,"ozone",itap,zx_tmp_3d,
    2235      .                                   iim*(jjm+1)*klev,ndex)
     2261     .                                   iim*(jjm+1)*klev,ndex3d)
    22362262c
    22372263      IF (nqmax.GE.3) THEN
    … …  
    22412267         WRITE(str2,'(i2.2)') iq
    22422268         CALL histwrite(nid_mth,"trac"//str2,itap,zx_tmp_3d,
    2243      .                                   iim*(jjm+1)*klev,ndex)
     2269     .                                   iim*(jjm+1)*klev,ndex3d)
    22442270      ELSE
    22452271         PRINT*, "Trop de traceurs"
    … …  
    22492275      ENDIF
    22502276c
    2251       if (ok_sync) call histsync
     2277      if (ok_sync) then
     2278        call histsync(nid_mth)
     2279      endif
    22522280      ENDIF
    22532281c
    22542282      IF (ok_instan) THEN
    22552283c
    2256       ndex(1) = 0
     2284      ndex2d = 0
     2285      ndex3d = 0
    22572286c
    22582287c Champs 2D:
    22592288c
     2289         i = NINT(zout/zsto)
     2290         CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
     2291         CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
     2292c
     2293         i = NINT(zout/zsto)
     2294         CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
     2295         CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
     2296
    22602297      DO i = 1, klon
    22612298         zx_tmp_fi2d(i) = paprs(i,1)
    22622299      ENDDO
    22632300      CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d)
    2264       CALL histwrite(nid_ins,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2301      CALL histwrite(nid_ins,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    22652302c
    22662303      CALL gr_fi_ecrit(1, klon,iim,jjm+1, toplw,zx_tmp_2d)
    2267       CALL histwrite(nid_ins,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex)
     2304      CALL histwrite(nid_ins,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d)
    22682305c
    22692306c Champs 3D:
    … …  
    22712308      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t_seri, zx_tmp_3d)
    22722309      CALL histwrite(nid_ins,"temp",itap,zx_tmp_3d,
    2273      .                                   iim*(jjm+1)*klev,ndex)
     2310     .                                   iim*(jjm+1)*klev,ndex3d)
    22742311c
    22752312      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, u_seri, zx_tmp_3d)
    22762313      CALL histwrite(nid_ins,"vitu",itap,zx_tmp_3d,
    2277      .                                   iim*(jjm+1)*klev,ndex)
     2314     .                                   iim*(jjm+1)*klev,ndex3d)
    22782315c
    22792316      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, v_seri, zx_tmp_3d)
    22802317      CALL histwrite(nid_ins,"vitv",itap,zx_tmp_3d,
    2281      .                                   iim*(jjm+1)*klev,ndex)
     2318     .                                   iim*(jjm+1)*klev,ndex3d)
    22822319c
    22832320      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, zphi, zx_tmp_3d)
    22842321      CALL histwrite(nid_ins,"geop",itap,zx_tmp_3d,
    2285      .                                   iim*(jjm+1)*klev,ndex)
     2322     .                                   iim*(jjm+1)*klev,ndex3d)
    22862323c
    22872324      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, pplay, zx_tmp_3d)
    22882325      CALL histwrite(nid_ins,"pres",itap,zx_tmp_3d,
    2289      .                                   iim*(jjm+1)*klev,ndex)
    2290 c
    2291       if (ok_sync) call histsync
     2326     .                                   iim*(jjm+1)*klev,ndex3d)
     2327c
     2328      if (ok_sync) then
     2329        call histsync(nid_ins)
     2330      endif
    22922331      ENDIF
    22932332c
Note: See TracChangeset for help on using the changeset viewer.