program info_SR95 

      implicit double precision (a-h,o-z)
      logical in_sol,select(8)
      dimension time(2),rlat(2),rlon(2),u(3,2),v(3,2)

      data select /8*.true./

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.41959

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

      do i = 1,2

c...Call SR95.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) in cm
c       tld    - long period tide in cm        
c       tide   - predicted ocean tide in cm

      call perth(rlat(i),rlon(i),time(i),select,tpd,in_sol)

c...Compute the long period tide height
      call lpeqmt (time(i)*86400.d0,rlat(i),tld)

       tide = tpd + tld

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

      enddo

      stop
      end 

c...The output data is:
c   -65.44040, 277.80185, 49354.41959, 18.82200, 21.12646, -2.30446
c   -63.10150, 291.97490, 49354.50178, 33.27564, 35.46877, -2.19313


c  SR95.0/.1
c  ---- -  -
c       The SR95.0 model  (version  950308)  is  an  update  of  the  one
c  described  by  Schrama and Ray [1994].  Most of the details concerning
c  the data processing and other aspects can be found in that paper.  The
c  two  major  changes  are  that  the  present  model  was  derived as a
c  correction to the Grenoble FES94.1 model, whereas the published  paper
c  shows corrections to the Schwiderski and Cartwright-Ray models.  As in
c  FES94.1, the model is given on a 0.5   geographical  grid.   Secondly,
c  the load tide required for deriving the ocean tide from the geocentric
c  altimetric tide was computed in a rigorous manner, following  Appendix
c  A of Cartwright and Ray [1991].  Another small improvement was a minor
c  correction to the orbit,  which  was  implemented  in  a  semianalytic
c  fashion by Bettadpur and Eanes [1994]; the orbit therefore corresponds
c  to the one computed from the JGM-3 gravity model with a T/P-based tide
c  model  (as  in  the  second  release  of the GDR's).  Additionally, of
c  course,  more altimeter data have been used: both TOPEX  and  POSEIDON
c  altimetry from cycles 9 through 71 were processed for this version.
c       As in the original paper, a simple harmonic method was  used  for
c  deriving  the  tidal  solution.   Five  constituents  were solved for:
c  !M2!, !S2!, !N2!, !O1!, and !K1!.  The !Q1! and !K2! constituents were
c  adopted  (courtesy  of C. Le Provost) directly from the FES94.1 model.
c  For computing tidal height predictions (e.g., for  use  in  correcting
c  altimetry),  some  16  additional  minor constituents were included by
c  linear inference from the major constituents.

     The SR95.1 model is actually identical to SR95.0,  except  for  a
change  in the supplied tidal prediction software.  The newer programs
provide for the 16 minor tides that had been neglected in  the  SR95.0
software.