program test_nc_gwforcing
      use nc_gwforcing
      implicit none
      character(len=120) :: 
     |  file_bin,   ! binary file to read data
     |  file_txt,   ! text file to read dimensions and coord arrays
     |  file_nc     ! new netcdf file to create
      integer :: ntime

      file_txt = 'LonLatZT.txt.com'
      file_bin = 'forcefxy.com.bin'
      file_nc  = 'forcefxy.nc'
!
! Read dimensions and coordinate arrays from text file:
! This defines spatial coord arrays in the nc_gwforcing module, and
!   allocates and reads the time(ntime) coordinate in the module.
!
      call read_text(file_txt,ntime)
!
! Allocate and initialize data:
      call init_data(ntime)
!
! Read data from fortran binary file:
      call read_binary(file_bin,ntime)
!
! Create netcdf data set:
      call nc_write(file_nc,ntime)
!
! Read the new netcdf 
      call init_data(ntime) ! reallocate and init fu,fv
      call nc_read(file_nc)

      end program test_nc_gwforcing