# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 42*20.,
 sRef = 42*35.,
 bottomDragQuadratic = 0.002,
 viscAr=1.5E-4,
# This is the minimum viscosity used total=viscah+viscleith
# Often U*deltaX/2 for grid noise (Griffies and Hallberg,2000)
#  viscAh=1.1E4,
 viscAh=1.E3,
# viscA4=3.E2,
# viscA4=1.E4,
# This is the standard Leith Viscosity, prop to grad(curl v_h)
 viscC2leith=1.5,
# viscC4leith=.2,
# 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=1.5,
# viscC4leithD=.2,
# These are the minimum viscosities used with Smagorinsky
# total=viscastrain+viscsmag, etc.
# viscAstrain=3.E1,
# viscAtension=3.E1,
# This is harmonic Smagorinsky Coefficient
# 0.15=value from 3-d turbulence theory; 
# 3-4=stable oceanic values (Griffies and Hallberg, 2000)
# viscC2Smag=4.,
# This is the max Viscosity used, be it Smag or Leith, and it scales 
# with grid size and timestep
# Often 0.5 for stability (Griffies and Hallberg,2000)
# viscAhGridMax=0.5,
# viscA4GridMax=.5,
# viscA4max = 1.E11,
 useFullLeith=.TRUE.,
# viscA4=1.4E10,
 viscA4=1.5E10,
#- diffKhT=.5E3, 
# ********** DIFFUSION *************
#- diffKhS=.5E3,
 diffK4T=1.5E10,
 diffK4S=1.5E10,
 diffKrNrT=3*3.E-5,12*3.E-5,1*1.E-5,1*7.E-6,1*2.E-6,1*8.E-7,1*4.E-7,22*1.E-7,
 diffKrNrS=3*3.E-5,12*3.E-5,1*1.E-5,1*7.E-6,1*2.E-6,1*8.E-7,1*4.E-7,22*1.E-7,
 staggerTimeStep = .TRUE.
 tempAdvScheme = 33,
 saltAdvScheme = 33,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 nonLinFreeSurf = 0, 
# exactConserv = .TRUE.,
 rhonil=1035.,
 gravity=9.81,
# eosType = 'MDJWF',
 eosType = 'JMD95Z',
 ivdc_kappa=1.,
 useOldFreezing=.TRUE.,
# useRealFreshWaterFlux=.TRUE.,
 useNHMTerms=.TRUE.,
# turn on looped cells
 hFacMin=.1,
 hFacMindr=30.,
# set precision of data files
 readBinaryPrec=32,
 useSingleCpuIO=.TRUE.,
#debugMode  = .TRUE.,
 &

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

# Time stepping parameters
 &PARM03
  nIter0 =      0,
  nTimeSteps = 240,
#- 366.days ; 32.days (1 yr of integration will yield a reasonable flow field)
#endTime = 31622400.,
#endTime =  2764800.,
 deltaTmom = 1800.0,
 deltaTtracer= 1800.0,
 deltaTClock = 1800.0,
 abEps = 0.1,
#pChkptFreq=  2628000., 
 dumpFreq= 864000.,
 taveFreq= 2628000.,
# taveFreq=  31536000., 
 monitorFreq= 86400.,
 monitorFreq=  3600.,
# 1 months restoring timescale for temperature
# tauThetaClimRelax =  2592000.,
# 5 days restoring timescale for salinity
# tauSaltClimRelax =   432000.,
# periodicExternalForcing=.TRUE.,
# externForcingPeriod=86400.,
# externForcingCycle=31104000.,
 nTimeSteps = 12,
 monitorFreq= 1.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.TRUE.,
 delR = 10., 10., 10., 20., 20.,
        20., 30., 30., 30., 40.,
        40., 40., 40., 40., 40.,
        40., 40., 50., 50., 50.,
        60., 60., 70., 70., 80.,
        80., 90., 90.,100.,120.,
        130.,150.,150.,200.,200.,
        200.,200.,200.,200.,300.,
        300.,300.,
 ygOrigin = 30.4375,
 xgOrigin =-13.6875,
 dySpacing=.125,
 dxSpacing=.125,
 &

# Input datasets
 &PARM05
 bathyFile=      'bathy_MITgcm.bin',
 hydrogThetaFile='init_temp.bin',
 hydrogSaltFile= 'init_sal.bin',
# zonalWindFile=  'wind_stress_x.bin',
# meridWindFile=  'wind_stress_y.bin',
# thetaClimFile=  'lev_sst.bin',
# saltClimFile=   'lev_sss.bin',
# surfQFile=      'ncep_qnet.bin',
# EmPmRFile=      'ncep_emp.bin',
 &