program info_ORI_NAO 

      implicit double precision (a-h,o-z)
      logical in_sol,pseudo,rad
      dimension time(2),Itime(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 after Jan 1,1958, desired by the model
      Itime(i) = INT((time(i)-36204.d0)*86400.d0)

c...Call ORI_NAO model to compute the tide height
c       in_sol - logical denoting whether tide data exists
c                at desired location
c       tpd    - geocentric ocean tide height in cm
c       tld    - long period tide in cm        
c       tod    - loading tide in cm
c       tide   - predicted ocean tide in cm

      call tidetp(rlat(i),rlon(i),Itime(i),in_sol,tpd)

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

c...Compute the loading tide height
      call tptide(rlat(i),rlon(i),time(i)*86400.d0,tod,in_sol,
     & u,v,pseudo,rad)

       tide = tpd/10.d0 + (tld+tod)

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

      enddo

      stop
      end 

c...The output data is:
c   -65.44040, 277.80185, 1136196252, 23.24488, 27.57874, -2.30446, -2.02940
c   -63.10150, 291.97490, 1136203353, 28.92261, 33.71347, -2.19313  -2.59773


c  ORI
c  ---
c       The ORI model was developed by the Ocean  Research  Institute  at
c  the  University  of  Tokyo  and  the National Astronomical Observatory
c  [Matsumoto et al., 1995].  This model is based on harmonic analysis of
c  data from crossover points from the first 77 cycles of TOPEX altimetry
c  processed  using  the  JGM-2  orbits.   A  hydrodynamic  interpolation
c  scheme,  similar  to  that  of  Schwiderski,  was  employed  to  allow
c  interpolation between the  crossovers  and  for  high  latitude  areas
c  outside  the  T/P  domain,  although  admittedly  with less precision.
c  Tidal constants were computed for the eight major constituents  (!O1!,
c  !Q1!,  !P1!,  !K1!, !M2!, !S2!, !K2! and !N2!), while eight additional
c  terms (2!N2!, !Mu2!, !Nu2!, !L2!, !T2!, !M1!,  !J1!  and  O!O1!)  were
c  computed  by  interpolations  of  semidiurnal and diurnal admittances.
c  The model spatial resolution is !1^times^1! degree.