# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 15*20.,
 sRef = 15*35.,
 viscAr=1.E-3,
 viscAh=5.E5,
 diffKhT=1.E3, 
 diffKrT=3.E-5,
 diffKhS=1.E3,  
 diffKrS=3.E-5,
 rhonil=1035.,
 gravity=9.81,
 eosType = 'JMD95Z',
 ivdc_kappa=0.,
 implicitDiffusion=.TRUE.,
 tempAdvScheme=3,
 saltAdvScheme=3,
# useOldFreezing=.TRUE.,
 useRealFreshWaterFlux=.FALSE.,
 useCDscheme=.TRUE.,
# useNHMTerms=.TRUE.,
# turn on looped cells
 hFacMin=.3,
 hFacMindr=50.,
# set precision of data files
 readBinaryPrec=32,
 globalFiles=.TRUE.,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=500,
 cg2dTargetResidual=1.E-13,
 &

# Time stepping parameters
 &PARM03
 nIter0 =      2160000, 
# 2 years
# nTimeSteps =  1460,
# 4 months
 nTimeSteps =  240,
# 20 years
# #_time_steps=#_yrs*60*60*24*365/43200  
# 100  (3110400000) years of integration will yield a reasonable flow field
# startTime  = 0.,
# endTime    = 43545600000.,
 deltaTmom = 1200.0,
 tauCD =     321428.,
 deltaTtracer= 43200.0,
 deltaTClock = 43200.0,
# if you are using a version later than checkpoint45d on the main branch
# you can uncomment the following line and increase the time step 
# deltaTtracer and deltaTClock to 172800.0 as well to speed up the 
# asynchronous time stepping
 deltaTfreesurf = 172800.0,
 abEps = 0.1,
 pChkptFreq= 31104000,
# dumpFreq=   311040000.,
 adjDumpFreq= 864000.,
# taveFreq= 63072000.,
#ph cAdjFreq =-1.,
 cAdjFreq =0.,
 monitorFreq=6307200.,
# 2 months restoring timescale for temperature
 tauThetaClimRelax =  5184000.0,
# 6 months restoring timescale for salinity
 tauSaltClimRelax = 15552000.0,
# periodicExternalForcing=.TRUE.,
# externForcingPeriod=2592000.,
# externForcingCycle=31104000.,
 &                  

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

# Input datasets
 &PARM05
 bathyFile=      'bathymetry.bin',
 hydrogThetaFile='lev_t.bin',
 hydrogSaltFile= 'lev_s.bin',
 zonalWindFile=  'trenberth_taux_march.bin',
 meridWindFile=  'trenberth_tauy_march.bin',
 thetaClimFile=  'lev_sst_march.bin',
 saltClimFile=   'lev_sss_march.bin',
# surfQFile=      'ncep_qnet.bin',
# fresh water flux is turned off, uncomment next line to turn on
# (not recommened together with surface salinity restoring)
# EmPmRFile=      'ncep_emp.bin',
 &