# ====================
# | 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.
 globalfiles=.FALSE.
 hFacMinDr=50.,
 hFacMin=0.3,
#ph(
 hFacInf=0.1,
 hFacSup=5.,
# z* - begin
#    select_rStar=2,
    nonlinFreeSurf=2,
# * - end
#ph)
 implicitDiffusion=.true.,
 implicitViscosity=.true.,
 allowFreezing=.false.,
 viscA4     = 0.,
 viscA4Grid = 2e-2,
 viscAh     = 1.E2,
 bottomDragQuadratic = 0.002,
 vectorInvariantMomentum=.TRUE.,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
#ph(
 multiDimAdvection=.FALSE.,
 StaggerTimeStep=.FALSE.,
 useCDscheme=.TRUE.,
#
 inAdExact=.TRUE.,
 tempAdvScheme=30,
 saltAdvScheme=30,

#ph)
 debuglevel=1
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=300,
 cg2dTargetResidual=1.E-10,
 &

# Time stepping parameters
 &PARM03
## startTime=0.0,
## endTime=864000.,
 niter0=0,
#1 model year
 nTimeSteps=12960,
 abEps=0.1,
 deltaT=2400.0,
 cAdjFreq=0.,
 pChkptFreq=31536000.,
 chkptFreq= 31536000.,
 dumpfreq=2592000.,
 monitorFreq=2592000.,
 adjDumpFreq=2592000.,
#
 tauSaltClimRelax = 15552000.,
 forcing_In_AB = .FALSE.
 tauCD=172800.0,
 &

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