# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01

 viscA4=1e-5,
 viscAz=1e-5,
 diffK4T=1e-5,
 diffKzT=1e-5,
 diffK4S=1e-5,
 diffKzS=1e-5,

# viscAhGrid=2.E-2,
# viscA4Grid=1.E-3,
# viscAz=1e-5,
# diffK4T=1e-5,
# diffK4S=1e-5,

# useSmag3D = .TRUE.,

# bottomDragQuadratic = 0.0025,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.FALSE.,
 rhonil=1027.5,
 rhoConstFresh=999.8,
 f0=0.0,
 beta=0.0,
 eosType='JMD95Z',
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 hFacMin=0.02,
 nonHydrostatic=.TRUE.,
 tempAdvScheme=77,
 saltAdvScheme=77,
 staggerTimeStep=.TRUE.,
 buoyancyRelation = 'OCEANIC',

 useSingleCpuIO=.TRUE.,
 selectAddFluid = 1.,
 readBinaryPrec=64,
 writeBinaryPrec=32,
 debugLevel=-1,
 
 &
 
# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=300,
 cg2dTargetResidual=1.E-13,
 cg3dMaxIters=20,
 cg3dTargetResidual=1.E-8,
 &

# Time stepping parameters
 &PARM03
 nIter0=0,
 endTime=21600,
 deltaT=3,
 abEps=0.01,
 pChkptFreq=86400,
 chkptFreq=0,
 dumpFreq=0,
 monitorSelect=1,
 monitorFreq=900.,
 periodicExternalForcing=.FALSE.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 delXfile='dx.bin',
 delYfile='dy.bin',
 delRfile= 'dz.bin',
 &

# Input datasets
 &PARM05
 bathyFile='bathy.bin',
 hydrogThetaFile='temp_rink.bin',
 hydrogSaltFile ='salt_rink.bin',
 &