#:::====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
# momstepping=.FALSE.,
# tempstepping=.FALSE.,
# saltstepping=.FALSE.,
 Tref = 50*-1.9,
 Sref = 50*34.4,
 viscAr=1.E-3,
#viscAh=600.0,
 viscAhGrid=0.2,
 viscA4Grid=0.02,
 no_slip_sides=.TRUE.,
 no_slip_bottom=.FALSE.,
 diffKhT= 100.0,
#diffKhT= 600.0,
 diffKrT=5.E-5,
 diffKhS= 100.0, 
#diffKhS= 600.0,
 diffKrS=5.E-5,
 bottomDragQuadratic=2.5E-3,
 staggerTimestep = .true.,
 tempAdvScheme=30,
 saltAdvScheme=30,
 vectorInvariantMomentum = .true.,
### momImplVertAdv=.true.,
 tempImplVertAdv=.true.,
 saltImplVertAdv=.true.,
 eosType='JMD95Z', 
#tAlpha = 2.E-4,
#sBeta  = 7.4E-4,
 HeatCapacity_cp = 3974.0,
 rhoConst=1000.,
 rhoNil=1000.,
 gravity=9.81,
 gBaro=9.81,
 implicitDiffusion = .true.,
 implicitViscosity = .true.,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 useRealFreshWaterFlux = .true.,
 exactConserv=.TRUE.,
# start nlfs with rstar
#select_rStar=2,
nonlinFreeSurf=4,
hFacInf=0.2,
hFacSup=2.0,
# end nlfs with rstar
 hFacMin=0.2,
#integr_GeoPot = 1,
#ivdc_kappa = 0.005.,
#implicitDiffusion=.true.,
#implicitViscosity=.true.,
#useCDscheme=.true.,
#nonHydrostatic=.true.,
#allowFreezing = .true.,
 readBinaryPrec=64,
#writeBinaryPrec=64,
 debuglevel = 1,
 selectCoriMap = 0,
 f0 = 0.0,
# useSingleCpuIO=.TRUE.,
 globalFiles=.TRUE.,
 &

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

#Time stepping parameters
 &PARM03
 nTimeSteps=518400,
 startTime=0.,
 deltaT=600.0,
 forcing_In_AB  = .false.,
 abEps=0.1,
 cAdjFreq = -1.,
 pChkptFreq=864000.0,
 chkptFreq=864000.0,
 monitorFreq=2592000.,
 dumpFreq=864000.,
 &

# Gridding parameters
 &PARM04
# usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.TRUE.,
 delR=100*10.,
 delX=3*.125,
 delY=200*.0078125
 xgOrigin = -105.5,
 ygOrigin =  -75.4457,
# delX = 200*
 &

# Input datasets
 &PARM05
 bathyFile='bathymetry.pig.bin',
 hydrogThetaFile='theta.init',	
 hydrogSaltFile='salt.init',
 pSurfInitFile='etainit.round.bin'
 
 &