# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 sRef=150*0.0,
## -----------------------------------------
## Horizontal viscosity and diffusivity:
## -------------------------------------------
 viscAhGrid = 0.002,
 diffKhT=1.E-5
#--------Heat (default is 3994 J/(Kg K)). The code uses a rhoConst = 999.8 kg m-3
 heatCapacity_Cp = 4210,
# -------------3D smagorinski
 useSmag3D=.TRUE.,
 smag3D_coeff=0.0005,
# -----------------------------------
 usejamartwetpoints=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 highOrderVorticity=.TRUE.
 useEnergyConservingCoriolis=.TRUE.,
## -----------------------------------------
## Vertical viscosity and diffusivity:
## -----------------------------------------
 viscAz=1.0E-6, 
 diffKzT=1.4E-7,
# Warning if using implicit viscosity with NH
 implicitViscosity=.FALSE.,
 implicitDiffusion=.FALSE.,
## --------------------------------------------
## Salinity?
## ---------------------------------------------
# do not use salt for now
 saltStepping=.FALSE.,
 diffKhS=0,
 diffKzS=1.43E-9,
## --------------------------------------------
## Hydrostatic pressure only?
## ---------------------------------------------
 nonHydrostatic=.TRUE.,
## --------------------------------------------
## Free Surface
## --------------------------------------------
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 hFacMin=1,
 exactConserv=.TRUE.,
## --------------------------------------------
## Drag
## --------------------------------------------
 no_slip_bottom=.TRUE.,
 no_slip_sides=.TRUE.,
## --------------------------------------------
## Equation of state:
## --------------------------------------------
 eosType='MDJWF',
## --------------------------------------------
## Others:
## --------------------------------------------
 allowFreezing=.FALSE.,
 tempAdvScheme=33,
 saltAdvScheme=33,  
 f0 = 1.065e-4,
 beta=0.E-11,
 readBinaryPrec=64,
 staggerTimeStep=.TRUE.,
## --------------------------------------------
## Global files?
## --------------------------------------------
#- not safe to use globalFiles in multi-processors runs
# globalFiles=.TRUE.,
 plotLevel=0,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=300,
 cg2dTargetResidual=1.E-13,
 cg3dMaxIters=20,
 cg3dTargetResidual=1.E-8,
 &

# Time stepping parameters
 &PARM03
 nIter0=0,
#-------36 hours
 nTimeSteps = 259200,
 deltaT=0.5,
 abEps=0.01,
 pChkptFreq=3600.0,
 chkptFreq=0.0,
 dumpFreq=600.0,
 monitorSelect=2,
 monitorFreq=3600.,
 periodicExternalForcing = .TRUE.,
 externForcingPeriod = 600.0,
 externForcingCycle  = 172800,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 dYspacing=2.0,
 dXspacing=2.0,
 delZ=40*0.05,80*0.1,30*0.2,
 &

# Input datasets
 &PARM05
 bathyFile='bathyBasin.bin',
 hydrogThetaFile='T.init',
 surfQfile='Qnet.forcing',
 zonalWindFile = 'ustress_run1.forcing',
#  zonalWindFile = 'ustress_run2.forcing',
#  zonalWindFile = 'ustress_run3.forcing',
#  zonalWindFile = 'ustress_run4.forcing',
#  zonalWindFile = 'ustress_run5.forcing',
#  zonalWindFile = 'ustress_run6.forcing',
 meridWindFile = 'vstress.forcing',
 &
