#
# ******************************
# Model parameters
# Continuous equation parameters
# ******************************
 &PARM01
 tRef= 24.0 , 23.0 , 22.0 , 21.0 , 20.0 ,
       19.0 , 18.0 , 17.0 , 16.0 , 15.0 ,
       14.0 , 13.0 , 12.0 , 11.0 , 10.0 ,
        9.0 ,  8.0 ,  7.0 ,  6.0,   5.0 ,
        4.0 ,  3.0 ,  2.0 ,

 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.71, 34.70, 34.69,

 no_slip_sides=.false.,
 no_slip_bottom=.false.,

# Harmonic Values From ECCO
# viscAz=1.E-3,
# viscAh=1.E4,
 
# diffKhT=1.E2,
# diffKzT=1.E-5,
# diffKhS=1.E2,
# diffKzS=1.E-5,

# Baylor Suggestion #1 about an order of mag. less
# and diffusion off.
# actually ~ 6 times less in the horizontal..

 viscAz=1.E-3,
 viscAh=1.5E3,
 
# diffKhT=0,
# diffKzT=0,
# diffKhS=0,
# diffKzS=0,

# Baylor Suggestion #2 

# BIHARMONIC
#viscA4     = 0.,
#viscA4 = 1.39e-7,

# the same damping on the smallest scales. 
# viscA4 = 3.225e11,

# viscA4Grid = 2e-3,

# a bit more than that to really get rid of that noise.
 viscA4Grid = 1e-1,
 
 diffK4T=0.,
 diffK4S=0.,

# Advection Scheme 33 is itself a bit diffsive.
 tempAdvScheme=33,
 saltAdvScheme=33,

 StaggerTimeStep=.true.
 multiDimAdvection=.true.

 beta=1.E-11,
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 gravity=9.81,
 gBaro=9.81,
 rigidLid=.FALSE.,
 implicitFreeSurface=.true.,
 eosType='JMD95Z',
 readBinaryPrec=32,
 writeBinaryPrec=32,
 saltStepping=.TRUE.,
 tempStepping=.TRUE.,
 momStepping=.TRUE.,
 implicitDiffusion=.true.,
 implicitViscosity=.true.,
 globalFiles=.false.,
 useSingleCpuIO=.TRUE.,
 useCDscheme=.FALSE.,
 exactConserv=.false.,
 debugLevel=0, 
 vectorInvariantMomentum=.TRUE.,
 /

# **************************
# Elliptic solver parameters
# **************************
 &PARM02

# cg2dMaxIters=1000,
 cg2dMaxIters=100,
# cg2dTargetResidual=1.E-11,
 cg2dTargetResidual=1.E-5,
 /

# ************************
# Time stepping parameters
# ************************
 &PARM03
#  calendarDumps=.TRUE.,
calendarDumps=.FALSE.,

#
#  nIter0 = 1,

 nIter0 = 0,
# nTimeSteps = 1,
# one week = 168
# nTimeSteps=168,

# just 12 time stpes
# nTimeSteps = 12,


# 1 month 30.5*24,
#  nTimeSteps=732,
# 1 month 31*24,
#  nTimeSteps=750

# 2 months 30.5*24*2
# nTimeSteps=1464,

#nTimeSteps = 7300,

# about 10 months.
# nTimeSteps=7300,

# 1 year
# 365 days
  nTimeSteps=8760,

# 366 days
# nTimeSteps=8785,

# what const. did
#  nTimeSteps=8784,

# 5 years
# nTimeSteps=43920,

# 2 years
# nTimeSteps=17440,
# 11 years
# nTimeSteps=96398,
# 12 years
# nTimeSteps=105120,
# 13 years
# nTimeSteps=113880,
#
 forcing_In_AB = .FALSE.,

# deltaTtracer=600.0,
# deltaTClock =600.0,
# deltaTmom =  600.0,
 
 deltaTmom  =  3600.0,
 deltaTtracer= 3600.0,
 deltaTClock = 3600.0,

# deltaTmom  = 200.0,
# deltaTtracer=200.0,
# deltaTClock =200.0,


 cAdjFreq=0.,
 abEps=0.1,

# IAN
# every time step
#  dumpFreq = 1.,
 
#  dumpFreq = 3600.,
#  taveFreq = 3600.,
#  adjDumpFreq = 3600.,
#  chkptFreq = 3600.,
#  pChkptFreq = 3600.,
#  monitorFreq = 3600.,

# Daily

# dumpFreq = 86400.,
# dumpFreq = 1.,
# monitorFreq=1.,

# taveFreq = 86400.,
# adjDumpFreq = 86400.,
# chkptFreq = 86400.,
# pChekptFreq = 86400.,
# monitorFreq = 86400.,

# Weekly
# dumpFreq = 604800.,
# taveFreq = 604800.,
# adjDumpFreq = 604800.,
# chkptFreq  = 604800.,
# pChkptFreq = 604800.,

# two weeks
# dumpFreq = 54900.,
# taveFrew = 54900.,

# monthly. 30.5 days
 dumpFreq = 2635200.,
 taveFreq = 2635200.,

# monthly 31 days
# dumpFreq = 2678400.,
# TAVeFreq = 2678400.,

# ENDIAN

# year?
# monitorFreq = 31536000.0,
 /

# *******************
# Gridding parameters
# *******************
 &PARM04
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.FALSE.,
 usingCurvilinearGrid=.TRUE.,
 delZ=10.,10.,15.,20.,20.,25.,35.,50.,75.,100.,150.,200.,275.,350.,415.,450.,500.,500.,500.,500.,500.,500.,500.,
 horizGridFile = 'a',
 &
/

# **********
# Data Files
# **********
 &PARM05

# IAN

#  bathyFile='ib_closed.bin',
bathyFile='ib_open_15.bin',

#hydrogSaltFile = 'salt_init.from_pickup',
#hydrogThetaFile = 'theta_init.from_pickup',

#uVelInitFile = 'u_init.from_pickup',
#vVelInitFile = 'v_init.from_pickup',
#pSurfInitFile = 'eta_init.from_pickup',

# 22nd Jan 
hydrogSaltFile = 'salt_init.ecco_month_9',
hydrogThetaFile = 'theta_init.ecco_month_9',

uVelInitFile = 'u_init.ecco_month_9',
vVelInitFile = 'v_init.ecco_month_9',


# END IAN

 /