program make_uniform_mlod

! Description:  Creates an "mlod" file containing uniform mixed-layer ocean
!               depths of 100m
!
! Written by:   Steve Phipps  8 November 2002

  implicit none

  ! Define parameters
  integer, parameter :: lon=64, lat2=56

  ! Declare local variables
  integer :: i
  real, dimension(lon) :: arlon
  real, dimension(lat2) :: arlat
  real, dimension(lon, lat2) :: hmo2
  character(len=80) :: h2

  ! Create some longitude and latitude values, even though these are not used
  do i = 1, lon
     arlon(i) = 5.625 * real(i-1)
  end do
  do i = 1, lat2
     arlat(i) = (2.0 * 87.561 / 55.0) * real(i-1) - 87.561
  end do

  ! Specify mixed-layer ocean depth of 100m at all gridpoints
  hmo2 = 100.0

  ! Specify a header
  h2 = "Uniform mixed-layer ocean depth of 100m"

  ! Write data to output file
  open (unit=10, file="mlod", form="formatted", status="new", action="write")
  write (10, *) lat2
  write (10, '(6e12.5)') arlat
  write (10, *) lon
  write (10, '(6e12.5)') arlon
  write (10, *) h2
  write (10, '(6e12.5)') hmo2
  close (10)

end program make_uniform_mlod