PROGRAM TP20GR 
C
      PARAMETER(IXS=112,IYS=81,KS=27)
      INTEGER MSKSP(IXS,IYS)
      DIMENSION XST(IXS),XSV(IXS),YST(IYS),YSV(IYS),ZS(KS)
      dimension t(ixs,ks),s(ixs,ks),u(ixs,ks),v(ixs,ks),tcor(ixs,ks)
      dimension taux(ixs),tauy(ixs),spc(6),ps(ixs)
C
         read(iun,end=799) imtsp,jmtsp,kmsp,xst,xsv,yst,ysv,zs,msksp,spc
         nflds = spc(1)
         ndays = spc(2)
         iyr   = spc(3) - 1900
         imo   = spc(4)
         iday  = spc(5)
         ihr   = spc(6)
C
C---- READ DATA 
C
         DO 70 J=1,jmtsp
           read (iun,end=2000) u,v,t,s,tcor,taux,tauy,ps
 70      continue
      stop
      END

xst --- Temp. grid of e-w direction;
xsv --- vel.  grid of e-w direction;
yst --- Temp. grid of n-s direction;
ysv --- vel.  grid of n-s direction;
zs  --- vertical levels (in cm.)
msksp --- ocean mask: 0 = land; n= no. of active ocean model levels;
          where 0 < n < 28
imtsp,jmtsp,kmsp --- dimensions of the model domain.
nflds --- no. of 2-d fields in each record (=5), i.e. u,v,t,s, and tcor

u,v --- zonal, meridional velocity components;
t,s --- temperature, and salinity;
tcor --- average temperature correction due to data assim. per hour.
taux,tauy --- zonal meridional stress components;
ps --- surface presure deviation at rigid lid, 
  note: sea level deviation = ps/980.0 in cm.

each read of u,v,t,s read in a e-w slab of model fields (x-z);
each read of taux,tauy,ps read in a one dimensional data in e-w
direction for the jth latitude.

the ocean re-analysis files are in ieee binary, fortran readable for 
most of Unix work stations (Sun, Sgi, Apollo, Dec station etc.)


References:
Ji, M., A. Leetmaa and J. Derber, 1995: An Ocean Analysis System for Seasonal
      to Interannual Climate Studies. Mon. Wea. Rev., 123, 460-481.

Ji, M., T.M. Smith, 1995: Ocean model response to temperature data assimilation
      and varying surface wind stress: Intercomparisons and implications for
      climate forecast. Mon. Wea. Rev, 123, 1811-1821.

Derber, J.D., and A. Rosati, 1989: A Global oceanic data assimilation system. 
      J. Phys. Oceanogr., 19, 1333-1347.