program info_DW95 

      implicit
      
      dimension time(2),time1(2),rlat(2),rlon(2),u(3,2),v(3,2)

c...Read sample input data file to get time and position. 
c    rlat - north latitude (in degrees, -90 to  90)
c    rlon - east longitude (in degrees,  0  to 360)
c    time - desired time,  Modified Julian Date in days. e.g.,
c     Jan 1, 1992 is 48622.

      rlat(1) = -65.44040
      rlon(1) = 277.80185 
      time(1) = 49354.41960

      rlat(2) = -63.10150
      rlon(2) = 291.97490
      time(2) = 49354.50178

      do i = 1,2
c...Convert time to seconds after Jan. 1,1992,0 hours 
c...which desired by the model(MJD = 48622)
      time1(i) = (time(i)-48622.d0)*86400.d0

c...Call DW95.0 model to compute the tide height
c       in_sol - logical denoting whether tide data exists
c                at desired location
c       tpd    - computed ocean tide height(ocean+load+long period) in cm
c       tide   - predicted ocean tide in cm

      call tptide_cur(rlat(i),rlon(i),time1(i),1,tpd,tjuk,in_sol)

      tide = tpd

      if (in_sol) then
        write(6,*) rlat(i),rlon(i),time(i),tide,tpd
      else
        write(6,*) 'No tide solution avaliable.'
      endif

      enddo

      stop
      end 

c...The output data is:
c   -65.44040, 277.80185, 63281053.43999, 24.68869, 24.68869 
c   -63.10150, 291.97490, 63288153.79200, 31.87873, 31.87873 


c  DW95.0/.1
c  ---- -  -
c       The DW95.0  ocean  tide  model  [Desai  and  Wahr,  1995]  is  an
c  empirical  ocean  tide  model  estimated  using data from exact repeat
c  cycles 10 to 78 of the T/P altimeter mission.  The data used  employed
c  the  JGM-3  orbit  computed at the University of Texas. In addition to
c  estimating a smooth response across the diurnal and semidiurnal  tidal
c  bands,  the  monthly,  fortnightly and termensual ocean tides are also
c  estimated. The model defines the diurnal admittance  function  as  one
c  with  the  expected  Free Core Nutation resonance removed.  The DW95.1
c  ocean tide model is a recent revision of the DW95.0 model and employed
c  additional cycles of T/P data.