****************
Model variables
****************
Model variables are defined in .h Fortran 77 files :
Domain variables (*grid.h*)
===========================
**grid.h** : Environmental twodimensional arrays associated with curvilinear horizontal coordinate system
 h : Model topography (bottom depth [m] at RHOpoints.)
 dh : Topograhy increment in case of moving bathymetry
 f : Coriolis parameter [1/s].
 fomn : Compound term, f/[pm*pn] at RHO points.
 angler : Angle [radians] between XIaxis and the direction to the EAST at RHOpoints.

 latr : Latitude (degrees_north) at RHO, U, and Vpoints.
 latu
 latv
 lonr : Longitude (degrees_east) at RHO, U, and Vpoints.
 lonu
 lonv
 xp : XIcoordinates [m] at PSIpoints.
 xr : XIcoordinates [m] at RHOpoints.
 yp : ETAcoordinates [m] at PSIpoints.
 yr : ETAcoordinates [m] at RHOpoints.
 pm : Coordinate transformation metric "m" [1/meters] associated with the differential distances in XI.
 pn : Coordinate transformation metric "n" [1/meters]associated with the differential distances in ETA.
 om_u : Grid spacing [meters] in the XI direction at Upoints.
 om_v : Grid spacing [meters] in the XI direction at Vpoints.
 on_u : Grid spacing [meters] in the ETAdirection at Upoints.
 on_v : Grid spacing [meters] in the ETAdirection at Vpoints.
 dmde : ETAderivative of inverse metric factor "m", d(1/M)/d(ETA).
 dndx : XIderivative of inverse metric factor "n", d(1/N)/d(XI).
 pmon_p : Compound term, pm/pn at PSIpoints.
 pmon_r : Compound term, pm/pn at RHOpoints.
 pmon_u : Compound term, pm/pn at Upoints.
 pnom_p : Compound term, pn/pm at PSIpoints.
 pnom_r : Compound term, pn/pm at RHOpoints.
 pnom_v : Compound term, pn/pm at Vpoints.
 rmask : Landsea masking arrays at RHO,U,V and PSIpoints (rmask,umask,vmask) = (0=Land, 1=Sea)
 umask
 vmask
 pmask : pmask=(0=Land, 1=Sea, 1gamma2 =boundary).
 reducu : reduction coefficient along xaxis for rivers sections
 reducv : reduction coefficient along yaxis for rivers sections
Barotropic variables (*ocean2d.h*)
===================================
**ocean2d.h** : 2D dynamical variables for fast mode
 zeta,rzeta : Free surface elevation [m] and its time tendency;
 ubar,rubar : Vertically integrated 2D velocity components in
 vbar,rvbar : XI and ETAdirections and their time tendencies;
Tridimensionnal variables (*ocean3d.h*)
========================================
**ocean3d.h** : 3D tracers dynamical variables for baroclinic mode
 u,v : 3D velocity components in XI and ETAdirections
 t : tracer array (temperature, salinity, passive tracers, sediment)
 Hz : level thickness
 z_r : depth at rho point
 z_w : depth at w point
 Huon : transport a U point
 Hvon : transport at V point
 We, Wi : vertical velocity (explicit, implicit)
 rho : density anomaly
 rho1 : potential density at 1 atm
Surface forcing (forces.h)
==========================
**forces.h** :
*Surface momemtum flux (wind stress) :*
 sustr : XI and ETAcomponents of kinematic surface momentum flux
 svstr : wind stresses) defined at horizontal U and Vpoints.dimensioned as [m^2/s^2].
*Bottom mometum flux :*
 bustr : XI and ETAcomponents of kinematic bottom momentum flux
 bvstr : (drag) defined at horizontal U and Vpoints [m^2/s^2].!
*Surface tracers fluxes :*
 stflx : Kinematic surface fluxes of tracer type variables at horizontal RHOpoints. Physical dimensions [degC m/s]  temperature; [PSU m/s]  salinity.
 dqdt : Kinematic surface net heat flux sensitivity to SST [m/s].
 sst : Current sea surface temperature [degree Celsius].
 dqdtg : Twotimelevel grided data for net surface heat flux
 sstg : sensitivity to SST grided data [Watts/m^2/Celsius] and sea surface temperature [degree Celsius].
 dqdtp : Twotimelevel point data for net surface heat flux
 sstp : sensitivity to SST grided data [Watts/m^2/Celsius] and sea surface temperature [degree Celsius].
 tsst : Time of sea surface temperature data.
 sss : Current sea surface salinity [PSU].
 tair : surface air temperature at 2m [degree Celsius].
 wsp : wind speed at 10m [degree Celsius].
 rhum : surface air relative humidity 2m [fraction]
 prate : surface precipitation rate [cm day1]
 radlw : net terrestrial longwave radiation [Watts meter2]
 radsw : net solar shortwave radiation [Watts meter2]
 patm2d : atmospheric pressure above mean seal level
 paref : reference pressure to compute inverse barometer effect
 srflx : Kinematic surface shortwave solar radiation flux [degC m/s] at horizontal RHOpoints
*Wind induced waves everything is defined at rhopoint :*
 wfrq : windinduced wave frequency [rad/s]
 uorb : xicomponent of waveinduced bed orbital velocity [m/s]
 vorb : etacomponent of waveinduced bed orbital velocity [m/s]
 wdrx : cosine of wave direction [non dimension]
 wdre : sine of wave direction [non dimension]
 whrm : (RMS) wave height (twice the wave amplitude) [m]
 wepb : breaking dissipation rate (\epsilon_b term) [m3/s3]
 wepd : frictional dissipation rate (\epsilon_d term) [m3/s3]
 wepr :roller dissipation rate (\epsilon_r term) [m3/s3]
 wbst : frictional dissipation stress (e_d k/sigma) [m2/s2]
*Wave averaged quantities :*
 brk2dx : xidireciton 2D breaking dissipation (rho)
 brk2de : etadirection 2D breaking dissipation (rho)
 frc2dx : xidireciton 2D frictional dissipation (rho)
 frc2de : etadirection 2D frictional dissipation (rho)
 ust2d : xidireciton Stokes transport (upoint)
 vst2d : etadireciton Stokes transport (vpoint)
 sup : quasistatic wave setup (rhopoint)
 calP : pressure correction term (rhopoint)
 Kapsrf : Bernoulli head terrm at the surface (rhopoint)
 brk3dx : xidireciton 3D breaking dissipation (rho)
 brk3de : etadirection 3D breaking dissipation (rho)
 ust : xidireciton 3D Stokes drift velocity (upoint)
 vst : etadireciton 3D Stokes drift velocity (vpoint)
 wst : vertical 3D Stokes drift velocity (rhopoint)
 Kappa : 3D Bernoulli head term (rhopoint)
 kvf : vertical vortex force term (K term, 3D, rhopoint)
 Akb : breakingwaveinduced additional diffusivity (wpoint)
 Akw : waveinduced additional diffusivity (rhopoint)
 E_pre : previous timestep value for Akw estimation (rho)
 frc3dx : xidireciton 3D frictional dissipation (rho)
 frc3de : etadirection 3D frictional dissipation (rho)