c.....This program serves as the driver routine to the
c.....subroutines that define the Desai and Wahr ocean tide model
c.....For provided lat lon, the m2 ocean tide height
c.....should be defined by:
c..... tide1 = coamplitude = 45.0142 cm
c..... tide2 = Greenwich phase = 123.8352 degrees
      program driver
      implicit none
      double precision lat, lon, time
      double precision tide1, tide2
      integer itype
      logical isdata

      open(50,file='tidepot_dw.ampm2',form='formatted',
     &access='sequential',status='old')
      open(60,file='cuocean_dw3.2.100',
     &form='formatted',access='sequential',status='old')
c.....MJD for Jan 1, 1992 = 48622.0d0
c.....NOTE: time not needed if ITYPE <> 1
      time=0.0d0
      itype=4
      lat=10.5d0
      lon=180.5d0
      call dwotide(lat,lon,time,itype,tide1,tide2,isdata)
      print*,lat,lon,tide1,tide2,isdata

      stop
      end

