# Model parameters
# Continuous equation parameters
 &PARM01
#tRef and sRef should be properly initialized and account for the number of layers 
 tRef=10., 8.,  4.,  2.,  2.,  2.,  2.,  2.,  2.,  2.,
 sRef=35., 34.5, 34., 34., 34., 34., 34., 34., 34., 34.,

 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.,
 useRealFreshWaterFlux=.TRUE.,
# ---------------------------------------------------------

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

 eosType='LINEAR',
 tAlpha=2.E-4,
 sBeta =7.4E-04,

# 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 =       1800.,
#  tauCD =         321428.,
 deltaTtracer   = 1800.,
#  deltaTClock    = 1800.,
#  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= 837.09,    837.09,    837.09,    837.09,    837.09,    837.09,    837.09,    837.09,    837.09,    837.09,
#  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_301x181x10.bin',
 thetaClimFile=  'lev_sst_301x181.bin',

# Salinity equation
 hydrogSaltFile= 'lev_initial_salinity_301x181x10.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', 
 &