*=*=*=*= WRITEFI.html =*=*=*=*
subroutine WRITEFI

subroutine WRITEFI


      subroutine WRITEFI(ngrid,nsoil,
     s  unit,ldrs,time,nbetat,co2ice,tsurf,tsoil,emis,q2)
      implicit none
c=================================================================
c
c  Ecriture des fichiers de redemarrage sous format DRS
c
c  L.Fairhead 04/93
c
c  parametres:
c     unit:                      unite logique d'ecriture
c     time:                      pas de temps de la simulation
c     vcov, ucov, h, pext, phis: champs a ecrire
c     nbetat:                    Nbre d'etats qui seront ecrits
c
c=================================================================
#include "dimensions.h"
#include "dimphys.h"
#include "drsdef.h"
#include "description.h"
#include "drscoorfi.h"

      integer ngrid,nsoil
      logical ldrs
      real co2ice(ngrid),tsurf(ngrid),tsoil(ngrid,nsoil),emis(ngrid)
      real q2(ngrid, nlayermx+1)
      real time
      integer nbetat
      integer unit
      integer ierr, setvdim, setname, putdat, putvdim, cluvdb
      integer it1, it2
      real zz

      IF(ldrs) THEN
c  Ecriture/extension de la coordonnee temps
      ierr = cluvdb()
      ierr = setname(descript,'temps','temps simule','jours',' ')
      ierr = putvdim(unit,1,time,it1,it2)
      write(6,*)' WRITEDEM: it1, it2, ',it1 ,it2

C  ecriture des champs

c co2ice
      ierr = cluvdb()
      ierr = setvdim(1,
     .  ' ','hor_coor',' ',' ',hor_coor(1),hor_coor(ngrid))
      ierr = setvdim(2,
     .  ' ','temps',' ',' ',time,time)
      ierr = setname(descript,'co2ice','CO2 ice cover ',
     .  'Pa',' ')
      ierr = putdat(unit,co2ice)

c tsurf
      ierr = cluvdb()
      ierr = setvdim(1,
     .  ' ','hor_coor',' ',' ',hor_coor(1),hor_coor(ngrid))
      ierr = setvdim(2,
     .  ' ','temps',' ',' ',time,time)
      ierr = setname(descript,'tsurf','Surface temperature  ',
     .  'K',' ')
      ierr = putdat(unit,tsurf)

c tsoil
      ierr = cluvdb()
      ierr = setvdim(1,
     .  ' ','hor_coor',' ',' ',hor_coor(1),hor_coor(ngrid))
      ierr = setvdim(2,
     .  ' ','surf_coor',' ',' ',
     .  surf_coor(1),surf_coor(nsoil))
      ierr = setvdim(3,
     .  ' ','temps',' ',' ',time,time)
      ierr = setname(descript,'tsoil','Soil temperature ',
     .  'K',' ')
      ierr = putdat(unit,tsoil)

c emis
      ierr = cluvdb()
      ierr = setvdim(1,
     .  ' ','hor_coor',' ',' ',hor_coor(1),hor_coor(ngrid))
      ierr = setvdim(2,
     .  ' ','temps',' ',' ',time,time)
      ierr = setname(descript,'emis','surface emissivity  ',
     .  ' ',' ')
      ierr = putdat(unit,emis)

c q2
      ierr = cluvdb()
      ierr = setvdim(1,
     .  ' ','hor_coor',' ',' ',hor_coor(1),hor_coor(ngrid))
      ierr = setvdim(2,
     .  ' ','vert2_coor',' ',' ',
     .  vert2_coor(1),vert2_coor(nlayermx+1))
      ierr = setvdim(3,
     .  ' ','temps',' ',' ',time,time)

      ierr = setname(descript,'q2','pbl wind variance',
     .  ' ',' ')
      ierr = putdat(unit,q2)

      ELSE
c        WRITE(unit) time
c!!! WARNING test 1+1=2
         zz=0.
c        WRITE(unit+1,*) 'zz'
c        WRITE(unit+1,'(3e25.15)') zz
c        WRITE(unit+1,*) 'co2ice'
c        WRITE(unit+1,'(3e25.15)') co2ice
c        WRITE(unit+1,*) 'tsurf'
c        WRITE(unit+1,'(3e25.15)') tsurf
c        WRITE(unit+1,*) 'tsoil'
c        WRITE(unit+1,'(3e25.15)') tsoil
c        WRITE(unit+1,*) 'emis'
c        WRITE(unit+1,'(3e25.15)') emis
c        WRITE(unit+1,*) 'q2'
c        WRITE(unit+1,'(3e25.15)') q2
         WRITE(unit) zz
         WRITE(unit) co2ice,tsurf,tsoil,emis,q2
      ENDIF

      return
      end