if (ok_hf) then

c   Comprendre comment marche el i=nint(zout/zsto)
c
      print*,'ACRITURE HF !!! ACRITURE HF !!! ACRITURE HF !!! '
      ndex2d = 0
      ndex3d = 0
c
c
c     i = NINT(zout/zsto)
c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
c     CALL histwrite(nid_hf,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d)
c
c     i = NINT(zout/zsto)
c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
c     CALL histwrite(nid_hf,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d)
C
      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
      CALL histwrite(nid_hf,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
c
      DO i = 1, klon
         zx_tmp_fi2d(i) = paprs(i,1)
      ENDDO
      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
      CALL histwrite(nid_hf,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
c
      DO i = 1, klon
         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
      ENDDO
      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
      CALL histwrite(nid_hf,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d)

      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
      CALL histwrite(nid_hf,"u850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)

      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
      CALL histwrite(nid_hf,"v850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)

      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
      CALL histwrite(nid_hf,"u500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)

      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
      CALL histwrite(nid_hf,"v500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)

      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
      CALL histwrite(nid_hf,"u200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)

      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
      CALL histwrite(nid_hf,"v200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)

      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
      CALL histwrite(nid_hf,"phi500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)

      if (ok_sync) then
        call histsync(nid_hf)
      endif

      endif