Weddell Sea regions extracted for Andrew Stewart on August 18, 2014

Region1
1/24 degree resolution:
lat -78 to -60
lat -62 to -15
daily-averaged u,v,w,s,t for 365 days

Region2
1/24 degree resolution:
lat -71.5 to -61.52
lon -15 to +50
daily-averaged u,v,w,s,t for 365 days

Region3
1/48 degree resolution:
lat -78 to -63
lon -62 to -20
daily-averaged u,v,w,s,t for 80 days

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

2-D output files are NXxNY Real*4 IEEE big-endian binary
3-D output files are NXxNYxNZ Real*4 IEEE big-endian binary

Grid information is in subdirectories "grid"

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

XC_NXxNY       :: longitude East of center of grid cell
XG_NXxNY       :: longitude East of southwest corner of grid cell
YC_NXxNY       :: latitude North of center of grid cell
YG_NXxNY       :: latitude North of southwest corner of grid cell
DXC_NXxNY      :: Cell center separation in X across western cell wall (m)
DXG_NXxNY      :: Cell face separation in X along southern cell wall (m)
DYC_NXxNY      :: Cell center separation in Y across southern cell wall (m)
DYG_NXxNY      :: Cell face separation in Y along western cell wall (m)
Depth_NXxNY    :: Model bathymetry (m)

AngleCS_NXxNY  :: cosine(alpha) relative to geographic direction at grid cell center
AngleSN_NXxNY  :: 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_NXxNY :: cosine of grid orientation angle at U point location
V2zonDir_NXxNY :: minus sine of  orientation angle at V point location

RAC_NXxNY      :: vertical face area of tracer cell (m^2)
RAS_NXxNY      :: vertical face area of v cell (m^2)
RAW_NXxNY      :: vertical face area of u cell (m^2)
RAZ_NXxNY      :: vertical face area of vorticity points (m^2)

hFacC_NXxNYxNZ :: mask of tracer cell (0 is land, >0 is wet)
hFacS_NXxNYxNZ :: mask of v cell (0 is land, >0 is wet)
hFacW_NXxNYxNZ :: 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)
oceFWflx.*.data net upward freshwater flux, >0 increases salinity (kg/m^2/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.

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