# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01

 viscAh=10.,
 viscAz=1.E-5,
 viscA4=0.,
 diffKhT=0.1,
 diffKzT=1.E-5,
 diffK4T=0.,
 diffKhS=0.1,
 diffKzS=1.E-5,
 diffK4S=0.,

 useFullLeith=.TRUE.,
## This uses the Strain-Tension form for calculating viscous terms
 useStrainTensionVisc=.TRUE.,

## This is the standard Leith Viscosity coefficient, 
## visc=(visc2leith/pi)**3 grad(curl v_h)
 viscC2leith=0.5,
 viscC4leith=0.05,

## This is an additional Leith Viscosity, prop to grad(div.v_h)
## It traps instabilities that have no signal in grad(curl v_h)
 viscC2leithD=0.5,
 viscC4leithD=0.05,

## for this grid
 viscC2Smag=0.5,
 viscC4Smag=0.05,

## This is the max Viscosity used, be it Smag or Leith, and it scales 
## with grid size and timestep
 viscAhGridMax=1.E-2,
 viscA4GridMax=1.E-9,
##
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
 tempAdvScheme=33,
 saltAdvScheme=33,
 staggerTimeStep=.TRUE.,
#
 tAlpha=2.E-4,
 sBeta =7.E-4,
 gravity=9.81,
 f0=0.0000101456,
 beta=2.2209E-11,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 nonHydrostatic=.FALSE.,
 bottomDragLinear=0.0002,
 eosType='MDJWF',
 useRealFreshWaterFlux=.TRUE.,
 exactConserv=.TRUE.,
## nonlinFreeSurf=4,
## select_rStar=2,
 hFacMin=0.5,
# hFacMinDz=5.,
 hFacInf=0.2,
 hFacSup=2.0,
 useSingleCpuIO=.TRUE
 readBinaryPrec=64,
 writeBinaryPrec=64,

 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-13,
 cg3dMaxIters=400,
 cg3dTargetResidual=1.E-13,

 &

# Time stepping parameters
 &PARM03
 startTime=0.,
 endTime=31536000.,
 deltaTmom=180.0,
 pChkptFreq=31536000.,	
 dumpFreq=432000.,				
# taveFreq=864000.,
 &

# Gridding parameters
 &PARM04
 usingCurvilinearGrid=.FALSE.,
 usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.FALSE.,
 dxSpacing=5055.960,
 dySpacing=5100.408,
 delZ=  5.0 , 5.0 , 5.0 , 5.0 , 5.0, 5.0 , 5.0 , 5.0 , 5.0 , 5.0,
        5.0 , 5.0 , 5.0 , 5.0 , 5.0, 5.0 , 5.0 , 5.0 , 5.0 , 5.0,
        5.0 , 5.0 , 5.0 , 5.0 , 5.0, 5.0 , 5.0 , 5.0 , 5.0 , 5.0,
        10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0,
        50.0 , 100.0 , 100.0 , 500.0 , 500.0 , 500.0 , 1000.0 , 1000.0 , 1000.0,

 radius_fromHorizGrid=6370.E3,
 &

# Input datasets
 &PARM05
 hydrogThetaFile='temp_dec_init_5km.bin',
 hydrogSaltFile='salt_dec_init_5km.bin',
 uVelInitFile='U_dec_init_5km.bin',
 vVelInitFile='V_dec_init_5km.bin',
 bathyFile='bathy_BoCo_5km_UPD.bin',
 checkIniTemp=.false.,
 checkIniSalt=.false.,
 &