# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 15*20.,
 sRef = 15*35.,
 viscAr= 1.77e-5,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
 diffKhT=0.,
 diffK4T=0.,
 diffKrT=1.51e-5,
 diffKhS=0.,
 diffK4S=0.,
 diffKrS=1.51e-5,
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 gravity=9.81,
 eosType='JMD95Z',
 readBinaryPrec=32,
 writeBinaryPrec=32,
 useSingleCPUio=.TRUE.
 globalFiles=.TRUE.
 hFacMinDr=50.,
 hFacMin=0.3,
 hFacInf=0.1,
 hFacSup=5.,
# z* - begin
    select_rStar=2,
    nonlinFreeSurf=4,
# * - end
 implicitDiffusion=.true.,
 implicitViscosity=.true.,
 allowFreezing=.false.,
  viscAh=0.,
# viscAhMax=unset
  viscAhGrid=0.,
  viscC2Leith=0.,
  viscA4=0.,
#  viscA4Max=0.,
  viscA4Grid=0.,
  viscA4GridMin=3.e-4,
  viscA4GridMax=4.e-3,
  viscC4Leith=0.2,
 bottomDragQuadratic = 0.002,
 tempAdvScheme=33,
 saltAdvScheme=33,
 StaggerTimeStep=.true.,
 multiDimAdvection=.true.,
 vectorInvariantMomentum=.TRUE.,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 debuglevel=0
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=100,
 cg2dTargetResidual=1.E-5,
 &

# Time stepping parameters
 &PARM03
 nIter0 =      0,
 nTimeSteps = 20,
 deltaT = 1200.0,
 abEps = 0.1,
 pChkptFreq= 311040000.,
 dumpFreq=   311040000.,
 dumpFreq=   864000.,
 taveFreq=   311040000.,
 taveFreq=   864000.,
 monitorFreq=1.,
# 2 months restoring timescale for temperature
 tauThetaClimRelax =  5184000.0,
# 6 months restoring timescale for salinity
 tauSaltClimRelax = 15552000.0,
 &

# 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',
 &