IF (ok_hf) THEN
c
        zout = dtime * REAL(NINT(86400./dtime*ecrit_hf))
        zsto = zout
        zsto1 = dtime
        PRINT*, 'La frequence de sortie instant. est de ', ecrit_hf
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

         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)

cccIM      CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat,
         CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
     .                 nhori, nid_hf)

         CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb",
     .                 klev, presnivs, nvert)
c

      IF(lev_histhf.GE.1) THEN
c
c        CALL histdef(nid_hf, "phis", "Surface geop. height", "-",
c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
c    .                "once", zsto,zout)
c
c        CALL histdef(nid_hf, "aire", "Grid area", "-",
c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
c    .                "once", zsto,zout)
c
         CALL histdef(nid_hf, "aireTER","Grid area CONT","-",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "once", zsto,zout)
c
         CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","-",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "once", zsto,zout)
c
         CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "-",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "once", zsto,zout)
c
c Champs 2D:
c
         CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
     .                "inst(X)", zsto,zout)
c
         CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
         CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
     .                "inst(X)", zsto,zout)
c
         print*,'ATTENTION METTRE AVE(X) POUR LES PRECIPS'

         CALL histdef(nid_hf, "rain", "Precipitation", "kg/m^2s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
     .                "ave(X)", zsto1,zout)
c
cIM ENSEMBLES BEG
c
         CALL histdef(nid_hf, "tsol", "Surface Temperature", "K",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
         CALL histdef(nid_hf, "slp", "Sea Level Pressure", "Pa",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
c
         CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
         CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
         CALL histdef(nid_hf, "wind10m","10-m wind speed","m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
         DO k=1, nlevENS 
          IF(clev(k).EQ."500") THEN
          CALL histdef(nid_hf, "phi"//clev(k),
     .                "Geopotential"//clev(k)//"mb", "m2/s2",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
          ENDIF !clev(k).EQ."500"
         ENDDO !k=1, nlevENS
c
      ENDIF !lev_histhf.GE.1
c
      IF(lev_histhf.GE.2) THEN
c
         CALL histdef(nid_hf, "cldt", "Total cloudiness", "-",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto1,zout)
c
c
         CALL histdef(nid_hf, "SWnetOR",
     .                "Sfce net SW radiation OR", "W/m2",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto1,zout)
c
c
         CALL histdef(nid_hf, "SWdownOR",
     .                "Sfce incident SW radiation OR", "W/m^2",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto1,zout)
c
         CALL histdef(nid_hf, "LWdownOR",
     .                "Sfce incident LW radiation OR", "W/m^2",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zsto1,zout)
      ENDIF !lev_histhf.GE.2
c
cIM ENSEMBLES END
c
      IF(lev_histhf.GE.3) THEN
c
         DO k=1, nlevENS
c
          CALL histdef(nid_hf, "t"//clev(k),
     .                 "Temperature"//clev(k)//"mb","K",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
          IF(clev(k).NE."500") THEN
          CALL histdef(nid_hf, "phi"//clev(k), 
     .                "Geopotential"//clev(k)//"mb", "m2/s2",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
          ENDIF
c
          CALL histdef(nid_hf, "q"//clev(k),
     .                 "Specific humidity"//clev(k)//"mb","kg/kg",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
         IF(1.EQ.0) THEN
         CALL histdef(nid_hf, "rh"//clev(k),
     .                 "Relative humidity"//clev(k)//"mb", "%",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
         ENDIF
c
          CALL histdef(nid_hf, "u"//clev(k),
     .                 "Zonal wind"//clev(k)//"mb","m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
          CALL histdef(nid_hf, "v"//clev(k),
     .                 "Meridional wind"//clev(k)//"mb","m/s",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
         ENDDO !nlevENS
         IF(1.EQ.0) THEN
         CALL histdef(nid_hf, "cdrm", " Momentum drag coef.", "-",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
c
         CALL histdef(nid_hf, "cdrh", "Heat drag coef.", "-",
     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     .                "inst(X)", zsto,zout)
         ENDIF !(1.EQ.0) THEN
c
      ENDIF !lev_histhf.GE.3
c
      IF(lev_histhf.GE.4) THEN
#define histhf3d
#ifdef histhf3d
#include "ini_histhf3d.h"
#endif
      ENDIF !lev_histhf.GE.4
c
c#define histhf3d
c#ifdef histhf3d
c#include "ini_histhf3d.h"
c#endif
c
         CALL histend(nid_hf)
c
      endif ! ok_hf