c.....This program serves as the driver routine to the
c.....subroutines that define the Desai and Wahr ocean tide model
c.....At the given time (Jan 1, 1992) and lat and lon
c.....the total elastic ocean (ocean +load) tide from the DW model
c.....should be tide1 = 16.2102 cm
      program driver
      implicit none
      double precision lat, lon, time
      double precision tide1, tide2
      integer itype
      logical isdata

      open(50,file='tidepot_dw.amprec',form='formatted',
     &access='sequential',status='old')
      open(60,file='cuelastic_dw3.2.100',
     &form='formatted',
     &access='sequential',status='old')
c.....MJD for Jan 1, 1992 = 48622.0d0
      time=48622.0d0
      itype=1
      lat=10.5d0
      lon=180.5d0
      call dwotide(lat,lon,time,itype,tide1,tide2,isdata)
      print*,lat,lon,tide1

      stop
      end
