# ====================
# | 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,
 gravity=9.81,
 rhoConst=1030.,
 rhoConstFresh=1000.,
 eosType='JMD95Z',
#allowFreezing=.TRUE.,
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 select_rStar=2,
 nonlinFreeSurf=4,
 hFacInf=0.2,
#hFacInf=0.1,
 hFacSup=2.0,
 useRealFreshWaterFlux=.TRUE.,
 temp_EvPrRn=0.,
 hFacMin=.1,
 hFacMinDr=20.,
 vectorInvariantMomentum=.TRUE.,
 staggerTimeStep=.TRUE.,
 readBinaryPrec=64,
 writeBinaryPrec=64,
 debugLevel=0,
 plotLevel=0,
 useSingleCpuIO=.TRUE.,
# tempAdvScheme=77,
# saltAdvScheme=77,
# multiDimAdvection=.TRUE.,
 &

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

# Time stepping parameters
 &PARM03
# pickupSuff='ckptA',
#0yr
# nIter0= 0,
 startTime=26282880000.0,
 endTime=26438400000.0,
# deltaT=    3600.,
 deltaTmom   =3600.,
 deltaTtracer=3600.,
 deltaTfreesurf= 3600.,
 deltaTClock =86400.,
#must be same as cpl_atmSendFrq?
 abEps = 0.1,
# ChkptFreq  =15552000.,
 pChkptFreq=155520000.,
# pChkptFreq=86400.,
# taveFreq=3110400000.,
# taveFreq    =933120000.,
 dumpFreq=0.,
 monitorFreq =0.,
# monitorFreq =1.,
 forcing_In_AB=.FALSE.,
# periodicExternalForcing=.TRUE.,
# externForcingPeriod=2592000.,
# externForcingCycle=31104000.,
# 2 years restoring timescale for temperature
# tauThetaClimRelax=  62208000.,
# 10 years restoring timescale for salinity
# tauSaltClimRelax = 311040000.,
# 2 months restoring timescale for temperature
# tauThetaClimRelax=  5184000.,
# 2 years restoring timescale for salinity
# tauSaltClimRelax = 62208000.,

 &

# Gridding parameters
 &PARM04
 usingCurvilinearGrid=.TRUE.,
 horizGridFile='grid_cs32',
 radius_fromHorizGrid=6370.E3,
# delR= 10., 40., 80., 130., 190.,
#       250., 250., 250., 250., 250.,
#       250., 250., 250., 250., 250.,
#^2950m
#
 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',
# uVelInitFile ='iniOcnU.Ex10.cpl_CTRL_CS32_15k.bin',
# vVelInitFile ='iniOcnV.Ex10.cpl_CTRL_CS32_15k.bin',
# pSurfInitFile ='iniOcnEta.Ex10.cpl_CTRL_CS32.bin',
# zonalWindFile=  'J-OFURO_TauX_steady.CS32.bin',
# meridWindFile=  'J-OFURO_TauY_steady.CS32.bin',
# thetaClimFile=  'SST.cos_steady.CS32.bin',
# saltClimFile=   'SSS.WOA2023_steady.CS32.bin',
# surfQnetFile   ='',
# EmPmRFile      ='',
 &