# Model parameters
# Continuous equation parameters
 &PARM01
#tRef and sRef should be properly initialized and account for the number of layers 
#  tRef=75*20.0,
#  sRef=75*35.0,
 tRef=20.,10.,8.,6.,
 sRef=10.,10.,10.,10.,
 viscAh=1.E2,
 viscAz=1.E-2,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
# ---------------------------------------------------------

# Equation for temperature (theta)
 tempStepping=.TRUE.,
 diffKhT=4.E2,
 diffKzT=1.E-2,
# ---------------------------------------------------------

# Equation for salinity (S)
 saltStepping=.TRUE.,
 diffKhS=4.E2,
 diffKzS=1.E-2,
#  rhoConst=1035.,
#  rhoConstFresh=1000.,
# ---------------------------------------------------------

#  rotationPeriod=86400.,
#  beta=1.E-11,

 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,

 eosType='LINEAR',
 tAlpha=2.E-4,
 sBeta =7.4E-04,
 exactConserv=.TRUE.,
#  useRealFreshWaterFlux=.TRUE.,

# set precision of data files
 readBinaryPrec=64,
 writeBinaryPrec=64,
 &

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

# Time stepping parameters
 &PARM03
#  startTime=0.,
#  endTime=1200000., 
#  deltaTmom=1800.0,
#  deltaTtracer=1800.0,
#  abEps=0.1,
#  pChkptFreq=17000.0,
#  chkptFreq=0.0,
#  dumpFreq=2592000.0,
#  monitorFreq=1.,

# nIter0=      0,
# nTimeSteps = 100,
startTime =          0.,
# endTime   =   86400.,
endTime   =   31104000.,

 deltaTmom =       3600.,
#  tauCD =         321428.,
 deltaTtracer   = 3600.,
#  deltaTClock    = 1800.,
#  deltaTClock    = 86400.,
#  deltaTfreesurf = 86400.,
 abEps = 0.1,
 pChkptFreq= 1728000.,
 dumpFreq=   86400.,
#  taveFreq=   864000.,
 monitorFreq=1.,

# # 2 months restoring timescale for temperature
#  tauThetaClimRelax=  0,
# #  tauThetaClimRelax=  5184000.,

# # 6 months restoring timescale for salinity
#  tauSaltClimRelax = 0,
# #  tauSaltClimRelax = 15552000.,

#  periodicExternalForcing=.TRUE.,

# # 30 days in seconds
#  externForcingPeriod=2592000.,

# # 365 year in seconds
#  externForcingCycle=31104000.,   
 &

# Gridding parameters
 &PARM04
 usingSphericalPolarGrid=.TRUE.,
#  xgOrigin=-100.,
#  ygOrigin=6.9938, 
#  delX=301*0.2667,
#  delY=181*0.2667,
 ygOrigin=0.,
 delX=60*1.,
 delY=60*1.,
 delZ=500.,500.,500.,500.,
#  delZ=    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,
#           111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,
#           111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,
#           111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,
#           111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,
#           111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,
#           111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,    111.61,
#           111.61,    111.61,    111.61,    111.61,    111.61,
 &

 &PARM05
 bathyFile='topog.box',
 hydrogThetaFile=,
 hydrogSaltFile=,
 zonalWindFile='windx.sin_y',
 meridWindFile=,

#  bathyFile='bathymetry_301x181_closed_lakes_seas__closed_domain.bin',
#  bathyFile='bathymetry_301x181_closed_lakes_seas.bin',

# Momentum equations 
#  zonalWindFile=  'wind_taux_301x181.bin',
#  meridWindFile=  'wind_tauy_301x181.bin',
#  zonalWindFile=  'wind_taux_301x181_box.bin',
#  meridWindFile=  'wind_tauy_301x181_box.bin',

# Temperature equation
#  hydrogThetaFile='lev_initial_temperature_301x181x75.bin',
#  thetaClimFile=  'lev_sst_301x181.bin',

# Salinity equation
#  hydrogSaltFile= 'lev_initial_salinity_301x181x75.bin',
#  saltClimFile=   'lev_sss_301x181.bin',
#  surfQnetFile=   'qnet_301x181.bin',
# # fresh water flux is turned on, comment next line to it turn off
# # (maybe better with surface salinity restoring)
#  EmPmRFile=      'E-P-R_301x181.bin',   

 the_run_name=   'natl_latlon', 
 &