# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef= 21.0 , 20.0 , 19.0 , 18.0 , 17.0 ,
       16.0 , 15.2 , 14.5 , 13.9 , 13.3 ,
       12.4 , 11.3 ,  9.9 ,  8.4 ,  6.7 ,
        5.2 ,  3.8 ,  2.9 ,  2.3 ,  1.8 ,
        1.5 ,  1.1 ,  0.8 ,  0.6,   0.6,
 sRef= 34.65, 34.75, 34.82, 34.87, 34.90,
       34.90, 34.86, 34.78, 34.69, 34.60,
       34.58, 34.62, 34.68, 34.72, 34.73,
       34.74, 34.73, 34.73, 34.72, 34.72,
       34.73, 34.73, 34.74, 34.74, 34.75,
 viscAr= 1.77e-5
 no_slip_sides=.TRUE.,
 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.
 hFacMinDr=50.,
 hFacMin=0.3,
 hFacInf=0.1,
 hFacSup=5.,
# z* - begin
    select_rStar=2,
    nonlinFreeSurf=4,
# * - end
 implicitDiffusion=.true.,
 implicitViscosity=.true.,
 allowFreezing=.false.,
 viscA4     = 0.,
 viscA4Grid = 2e-3,
 viscAh     = 0.,
 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
 startTime=0.0,
 endTime=639998400
 abEps=0.1,
 deltaT=2400.0,
 cAdjFreq=0.,
 pChkptFreq=31557600.
 tavefreq =2629800
 tauSaltClimRelax = 15552000.,
 forcing_In_AB = .FALSE.
 monitorFreq=240000
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.FALSE.,
 usingCurvilinearGrid=.TRUE.,
 delX=130*0.35979,
 delY=102*0.35979,
 delR= 1.000000e+01, 1.000000e+01, 1.000000e+01, 1.000000e+01,
       1.000000e+01, 1.143000e+01, 1.850000e+01, 2.747000e+01,
       3.883000e+01, 5.283000e+01, 7.000000e+01, 9.050000e+01,
       1.147300e+02, 1.428300e+02, 1.749300e+02, 2.108700e+02,
       2.503300e+02, 2.929000e+02, 3.418300e+02, 4.018300e+02,
       4.651700e+02, 5.335000e+02, 6.093300e+02, 6.901700e+02,
       7.885000e+02,
 phimin=-19.78818,
 thetamin=-22.66646,
 rSphere         = 6371.D3
 &

# Input datasets
 &PARM05
 bathyFile       = 'topog.bin_130_102',
 hydrogThetaFile = 'theta.phc_130_102_25_12',
 hydrogSaltFile  = 'salt.phc_130_102_25_12',
 &

# Open-boundaries
 &PARM06
 &