# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 sRef=89*0.0,
## -----------------------------------------
## Horizontal viscosity and diffusivity:
## -------------------------------------------
  viscAh=1.E-4,
  diffKhT=1.E-5
#--------Heat (default is 3994 J/(Kg K)). The code uses a rhoConst = 999.8 kg m-3
 heatCapacity_Cp = 4200,
# -------------3D smagorinski
 useSmag3D=.TRUE.,
 smag3D_coeff=0.0001,
#- to use Smag3D diffusivity, set non-dimensional Smag3D diffusion
#  coeff: smag3D_diffCoeff (default = 0 --> unused) as shown below:
# smag3D_diffCoeff= 1.E-4,
# -----------------------------------
 usejamartwetpoints=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 highOrderVorticity=.TRUE.
 useEnergyConservingCoriolis=.TRUE.,
## -----------------------------------------
## Vertical viscosity and diffusivity:
## -----------------------------------------
 viscAz=1.0E-6, 
 diffKzT=1.4E-7,
# -------implicit viscosity? 
 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=0.5,
 exactConserv=.TRUE.,
## --------------------------------------------
## Drag
## --------------------------------------------
 no_slip_bottom=.TRUE.,
 no_slip_sides=.TRUE.,
## --------------------------------------------
## Equation of state:
## --------------------------------------------
 eosType='JMD95P',
## --------------------------------------------
## Others:
## --------------------------------------------
 allowFreezing=.FALSE.,
 tempAdvScheme=33,
# tempAdvScheme=77,
# tempAdvScheme=7,
 saltAdvScheme=33,  
 f0=0.e-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,
#-------24 hours
 nTimeSteps = 172800,
 deltaT=0.5,
 abEps=0.01,
 pChkptFreq=3600.0,
 chkptFreq=0.0,
 dumpFreq=300.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,49*0.1,
 &

# Input datasets
 &PARM05
 bathyFile='bathyBasin.bin',
 hydrogThetaFile='T.init',
 surfQfile='Qnet.forcing',
 &
