c     if (ok_hf) then

c   Comprendre comment marche el i=nint(zout/zsto)
c
c     print*,'ACRITURE HF !!! ACRITURE HF !!! ACRITURE HF !!! '
      ndex2d = 0
      ndex3d = 0
c
      zsto = dtime
      zout = dtime * ecrit_hf
      itau_w = itau_phy + itap
c
c     IF(lev_histhf.GE.1) THEN
c
c Champs 3D:
c
      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
      CALL histwrite(nid_hf3d,"temp",itau_w,zx_tmp_3d,
     .                                   iim*jjmp1*klev,ndex3d)
c
      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
      CALL histwrite(nid_hf3d,"ovap",itau_w,zx_tmp_3d,
     .                                   iim*jjmp1*klev,ndex3d)
c
      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
      CALL histwrite(nid_hf3d,"vitu",itau_w,zx_tmp_3d,
     .                                   iim*jjmp1*klev,ndex3d)
c
      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
      CALL histwrite(nid_hf3d,"vitv",itau_w,zx_tmp_3d,
     .                                   iim*jjmp1*klev,ndex3d)
c
c     ENDIF !lev_histhf.GE.1
c
      if (ok_sync) then
        call histsync(nid_hf3d)
      endif

c     endif