program info_CSR 

      implicit double precision (a-h,o-z)
      logical in_sol,pseudo,rad
      dimension time(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.41959

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

      do i = 1,2
c...Convert time to seconds desired by the model
      time(i) = time(i)*86400.d0

c...Call CSR3.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 csrtptide(rlat(i),rlon(i),time(i),tpd,in_sol,u,v,pseudo,rad)

c...Compute the long period tide height
      call lpeqmt (time(i),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, 4264221852.576004, 19.87621, 22.18067, -2.30446
c   -63.10150, 291.97490, 4264228953.792007, 30.73671, 32.92984, -2.19313

c  CSR3.0
c  ---- -
c      This model, developed by Richard Eanes and colleagues at the
c  University of Texas, is basically a long wavelength adjustment to the
c  Grenoble FES94.1 hydrodynamic model.  Thereby, a tide model product is
c  produced  which preserves the long wavelength accuracy of T/P with the
c  detailed spatial resolution of the Grenoble model.
c     The model is based upon 89 cycles (2.4 years) of  T/P  altimetry.
c  First  diurnal  orthoweights were fit to the !Q1!, !O1!, !P1! and !K1!
c  constituents of the Grenoble hydrodynamical model FES94.1 [Le  Provost
c  et  al.,  1994],  and  semidiurnal  orthoweights were fit to the !N2!,
c  !M2!, !S2! and !K2!  constituents  of  Andersen's  "Adjusted  Grenoble
c  Model"  [Andersen  et  al.,  1995].   Tides  in the Mediterranean from
c  Canceil et al. [1995] were used in both tidal bands as  they  appeared
c  in  the Andersen Adjusted Grenoble model as well as in FES94.1 itself.
c  Radial ocean loading tides from the previous CSR2.0 model  were  added
c  to the Grenoble ocean tides to convert them to geocentric tides.  Then
c  T/P altimetry was used to solve for corrections to these  orthoweights
c  in  !3^times^3!  degree  spatial bins.  The orthoweight corrections so
c  obtained were then smoothed by convolution with  a  2-d  gaussian  for
c  which   the  full-width-half-maximum  (FWHM)  was  7.0  degrees.   The
c  smoothed orthoweight corrections were output  on  the  !0.5^times^0.5!
c  degree  grid  of  the  Grenoble  model  and then added to the Grenoble
c  values to obtain the new model with a global domain.   The  T/P  orbit
c  used  for  this  tide model development was computed at Texas and used
c  the JGM-3 gravity field and a dynamical ocean tide model based upon an
c  even earlier Texas solution (CSR1.6).