# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
# f0 = 1.409e-4,
# beta = 1.0e-11,
 &PARM01
 Tref = 102*1.0,
 Sref = 102*34.4,
# useSmag3D = .TRUE.,
# smag3D_coeff = 1e-2,
 viscAz=1.E-4,
 viscAh=1.E-3,
 diffKhT=1E-3,
 diffKzT=1E-4,
 diffKhS=1E-3,
 diffKzS=1E-4,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.FALSE.,
 selectBotDragQuadr = 2,
 bottomDragQuadratic=2.5E-3,
 eosType='LINEAR', 
 HeatCapacity_cp = 3974.0,
 rhoNil=1030.,
 rhoConst=1030.,
 gravity=9.81,
# convertFW2Salt = 33.4,
 convertFW2Salt = -1,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 hFacMin=0.05,
 hFacMindz = 0.0,
 nonHydrostatic=.FALSE.,
 readBinaryPrec=64,
 useCDScheme = .FALSE.,
 debugLevel = -1,
 f0 = -1.409e-4,
 beta = -1.0e-11,
# f0 = 0.0,
# beta = 0.0,
 implicitDiffusion = .TRUE.,
 implicitViscosity = .TRUE.,
 staggerTimeStep = .TRUE.,
 tempAdvScheme = 33,
 saltAdvScheme = 33,
 useSingleCPUio = .TRUE.,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-13,
 cg3dMaxIters=400,
 cg3dTargetResidual=1.E-13,
 &

# Time stepping parameters
# tauCD = 400000.,
 &PARM03
 nIter0=0,
 nTimeSteps=100000000,
 deltaT=30.0,
 abEps=0.1,
 cAdjFreq = -1.,
 pChkptFreq=0.0,
 chkptFreq=86400.0,
 dumpFreq=0.0,
 taveFreq=0.0,
 monitorFreq=3600.0,
 monitorSelect=1,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 ygOrigin = 0.0,
 delX=1*200.0,
 delY=200*50.0,
 delZ=102*1.0,
 seaLev_Z = -498.0,
 top_Pres = 5.03e6,
 &

# Input datasets
 &PARM05
 bathyFile='bathy.box',
 &

