# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef=15*20.,
 sRef=15*35.,
 viscAh =3.E5,
 viscAr =1.E-3,
 diffKhT=0., 
 diffK4T=0., 
 diffKrT=3.E-5,
 diffKhS=0.,  
 diffK4S=0., 
 diffKrS=3.E-5,
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 gravity=9.81,
 rhonil=1035.,
 eosType='JMD95Z',
#staggerTimeStep=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 select_rStar=2,
 nonlinFreeSurf=4,
 hFacInf=0.2,
 hFacSup=2.0,
 useRealFreshWaterFlux=.TRUE.,
 allowFreezing=.TRUE.,
 hFacMin=.1,
 hFacMinDr=20.,
 readBinaryPrec=64,
 globalFiles=.TRUE.
 useSingleCPUio=.TRUE.
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=200,
#cg2dTargetResidual=1.E-9,
 cg2dTargetResWunit=1.E-14,
 &

# Time stepping parameters
 &PARM03
 nIter0=0,
 nTimeSteps=20,
 deltaTmom   =1200.,
 deltaTtracer=1200.,
 deltaTfreesurf=1200.,
 deltaTClock =1200.,
 abEps = 0.1,
 forcing_In_AB=.FALSE.,
 dumpFreq    =24000.,
# 2 months restoring timescale for temperature
 tauThetaClimRelax = 5184000.,
# 2yrs restoring timescale for salinity
 tauSaltClimRelax = 62208000.,
 monitorFreq =1.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.FALSE.,
 usingCurvilinearGrid=.TRUE.,
 delR= 50., 70., 100., 140., 190.,
       240., 290., 340., 390., 440.,
       490., 540., 590., 640., 690.,
 &

# Input datasets
 &PARM05
 bathyFile      ='bathy_Hmin50.bin',
 hydrogThetaFile='lev_T_cs_15k.bin',
 hydrogSaltFile ='lev_S_cs_15k.bin',
 &