3x3-deg region south of Southern tip of India, centered at 3.5S, 78.5E
extracted from llc_4320 (global simulation on 1/48th-degree lat-lon-cap grid)
to a depth of 1 km.

2-D output files are 144x156 Real*4 IEEE big-endian binary
3-D output files are 144x156x52 Real*4 IEEE big-endian binary

============

Grid information is in subdirectory "grid"

thk90.mat        :: dpt90(1:52) and thk90(1:52) are vertical coordinates in m
                    of cell centers and thicknesses, respectively

XC_144x156       :: longitude East of center of grid cell
XG_144x156       :: longitude East of southwest corner of grid cell
YC_144x156       :: latitude North of center of grid cell
YG_144x156       :: latitude North of southwest corner of grid cell
DXC_144x156      :: Cell center separation in X across western cell wall (m)
DXG_144x156      :: Cell face separation in X along southern cell wall (m)
DYC_144x156      :: Cell center separation in Y across southern cell wall (m)
DYG_144x156      :: Cell face separation in Y along western cell wall (m)
Depth_144x156    :: Model bathymetry (m)

AngleCS_144x156  :: cosine(alpha) relative to geographic direction at grid cell center
AngleSN_144x156  :: sine(alpha) relative to geographic direction at grid cell center
                    alpha = angle of model uVel direction vs geographical East
                          = angle of model vVel direction vs geographical North
                    (AngleCS*uVelc - AngleSN*vVelc, AngleSN*uVelc + AngleCN*vVelc)
                    rotates model velocity to geographical coordinates, where
                    (uVelc,vVelc) is model velocity vector at center of grid cell
U2zonDir_144x156 :: cosine of grid orientation angle at U point location
V2zonDir_144x156 :: minus sine of  orientation angle at V point location

RAC_144x156      :: vertical face area of tracer cell (m^2)
RAS_144x156      :: vertical face area of v cell (m^2)
RAW_144x156      :: vertical face area of u cell (m^2)
RAZ_144x156      :: vertical face area of vorticity points (m^2)

hFacC_144x156x52 :: mask of tracer cell (0 is land, >0 is wet)
hFacS_144x156x52 :: mask of v cell (0 is land, >0 is wet)
hFacW_144x156x52 :: mask of u cell (0 is land, >0 is wet)

============

Available hourly output is:
Eta.*.data      sea surface height (m)
KPPhbl.*.data   mixing layer depth (m)
PhiBot.*.data   bottom pressure
SIarea.*.data   fractional ice-covered area [0 to 1]
SIheff.*.data   effective ice thickness (m)
SIhsalt.*.data  effective salinity (g/m^2)
SIhsnow.*.data  effective snow thickness (m)
SIuice.*.data   zonal (relative to grid) ice velocity, >0 from West to East (m/s)
SIvice.*.data   merid. (relative to grid) ice velocity, >0 from South to North (m/s)
Salt.*.data     salinity (g/kg)
Theta.*.data    potential temperature (deg C)
U.*.data        zonal (relative to grid) velocity, >0 from West to East (m/s)
V.*.data        merid. (relative to grid) velocity, >0 from South to North (m/s)
W.*.data        vertical velocity (m/s)
oceQnet.*.data  net upward surface heat flux (including shortwave), >0 decreases theta (W/m^2)
oceQsw.*.data   net upward shortwave radiation, >0 decreases theta (W/m^2)
oceSflux.*.data net upward salt flux, >0 decreases salinity (g/m^2/s)
oceTAUX.*.data  zonal (relative to grid) surface wind stress, >0 increases uVel (N/m^2)
oceTAUY.*.data  meridional (relative to grid) surf. wind stress, >0 increases vVel (N/m^2)

Please note that U, V, oceTAUX, oceTAUY, SIuice, and SIvice are aligned
relative to model grid, not geographical coordinates, and that they are
specified at the SouthWest C-grid velocity points.  All other scalar fields
are specified at the tracer point, i.e., the center of each grid box.

============