# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 Tref = -0.194175, -0.582524, -0.970874, -1.35922, -1.74757, 
        -2.13592, -2.52427, -2.91262, -3.30097, -3.68932, 
	-4.07767, -4.46602, -4.85437, -5.24272, -5.63107, 
	-6.01942, -6.40777, -6.79612, -7.18447, -7.57282, 
	-7.96117, -8.34951, -8.73786, -9.12621, -9.51456,
 sRef= 25*1.,
 viscAz=5.E-1,
 viscAh=0.0E0,
 no_slip_sides=.TRUE.,
 no_slip_bottom=.TRUE.,
#- a trick to recover old version of side-drag:
#sideDragFactor=0.,
#-
 bottomDragQuadratic=2.E-3,
 viscA4=0.E12,
 viscC4Leith=1.458198138065,
 useAreaViscLength=.TRUE.,
 diffKhT=0.0E0
 diffKzT=0.E-2,
 diffKhS=0.E0,
 diffKzS=0.E-2,
 f0=1.0e-4,
 beta=0.E-11,
 tAlpha=2.E-4,
 sBeta =0.E-4,
 gravity=9.81,
 rhonil=1030.0,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 eosType='LINEAR',
 hFacMin=0.2,
 nonHydrostatic=.FALSE.,
 multiDimAdvection=.TRUE.,
 tempAdvScheme=77,
 saltAdvScheme=77,
 readBinaryPrec=64,
 writeBinaryPrec=64,
 writeStatePrec=64,
 globalFiles=.TRUE.,
 staggerTimeStep=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 &

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

# Time stepping parameters
 &PARM03
 niter0=0,
#nTimeSteps=12000,
 nTimeSteps=20,
 deltaT=300.,
 abEps=0.1,
 pChkptFreq=0.0,
 chkptFreq=100000.0,
 dumpFreq=30000.,
#cAdjFreq=30.,
#monitorFreq=6000.,
 monitorFreq=300.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.FALSE.,
# delXfile='delXvar'
 delX=200*10000.0,
 delYfile='delYvar',
 delZ = 25*144.0, 
 &

# Input datasets
 &PARM05
 hydrogThetaFile='T.init',
 bathyFile='topog.slope',
 &