c     IF (ok_mensuelNMC) THEN
      IF (ok_mensuel) THEN
c
       IF(region.EQ.'3d') THEN
         zsto = dtime
         zout = dtime * ecrit_mth
       ELSEIF(region.EQ.'1d') THEN
         zsto = dtime
c        zout = dtime
         zout = dtime *4.
       ENDIF
c
         idayref = day_ref
         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
c
         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
         DO i = 1, iim
            zx_lon(i,1) = rlon(i+1)
            zx_lon(i,jjmp1) = rlon(i+1)
         ENDDO
         DO ll=1,klev
            znivsig(ll)=float(ll)
         ENDDO
         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
         CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
c071103
c    .                 1,iim,1,jjmp1, 
     .                 imin_debut,nbpti,
     .                 jmin_debut,nbptj,
c071103
     .                 itau_phy, zjulian, dtime, 
     .                 nhori, nid_nmc)
         write(*,*)'Mensuel NMC ', itau_phy, zjulian
         CALL histvert(nid_nmc, "presnivs", "Vertical levels", "mb",
     .                 klev, presnivs, nvert)
c        call histvert(nid_nmc, 'sig_s', 'Niveaux sigma','-',
c    .              klev, znivsig, nvert)
c
c Champs 2D:

c Champs interpolles sur des niveaux de pression du NMC
         CALL histdef(nid_nmc, "u1000",
     . "Zonal wind 1000mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v1000",
     .                "Meridional wind 1000mb", "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "u925", "Zonal wind 925mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v925", "Meridional wind 925mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "u850", "Zonal wind 850mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v850", "Meridional wind 850mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u700", "Zonal wind 700mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v700", "Meridional wind 700mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u600", "Zonal wind 600mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v600", "Meridional wind 600mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u500", "Zonal wind 500mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v500", "Meridional wind 500mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u400", "Zonal wind 400mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v400", "Meridional wind 400mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u300", "Zonal wind 300mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v300", "Meridional wind 300mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u250", "Zonal wind 250mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v250", "Meridional wind 250mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u200", "Zonal wind 200mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v200", "Meridional wind 200mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u150", "Zonal wind 150mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v150", "Meridional wind 150mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u100", "Zonal wind 100mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v100", "Meridional wind 100mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u70", "Zonal wind 70mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v70", "Meridional wind 70mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u50", "Zonal wind 50mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v50", "Meridional wind 50mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u30", "Zonal wind 30mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v30", "Meridional wind 30mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u20", "Zonal wind 20mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v20", "Meridional wind 20mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histdef(nid_nmc, "u10", "Zonal wind 10mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)

         CALL histdef(nid_nmc, "v10", "Meridional wind 10mb", 
     .                "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto,zout)
c
         CALL histend(nid_nmc)
c
c        ndex2d = 0
c
      ENDIF ! fin de test sur ok_mensuel