program make_zero_oflux

! Purpose
! -------
! Creates a flux restart file for the CSIRO coupled model, containing zero
! fluxes
!
! Usage
! -----
! make_zero_oflux <filename>
!
! History
! -------
! 2003 Sep 6	Steve Phipps	Original version

  implicit none

  ! Define constants
  integer, parameter :: outunit=11, lon=64, lat=28

  ! Declare variables
  character(len=80) :: outfile
  integer :: itt
  real, dimension(lon, lat, 2) :: data

  ! Set timestep number and fluxes to zero
  itt = 0
  data = 0.0

  ! Get name of output file
  call getarg(1, outfile)

  ! Open output file
  open(unit=outunit, file=outfile, status='new', form='unformatted', &
       action='write')

  ! Write data to file
  write (outunit) 'surface fluxes'
  write (outunit) itt
  write (outunit) data
  write (outunit) data
  write (outunit) data
  write (outunit) data
  write (outunit) 'end surface fluxes'

  ! Close output file
  close (unit=outunit)

end program make_zero_oflux