#
# ******************************
# 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=.TRUE.,
 viscAz=1.E-3,
 viscAh=1.E4,
 diffKhT=1.E2,
 diffKzT=1.E-5,
 diffKhS=1.E2,
 diffKzS=1.E-5,
 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.,
 useNHMTerms=.TRUE.,
 exactConserv=.FALSE.,
# exactConserv=.TRUE.,
 tempAdvScheme = 3,
 saltAdvScheme = 3,
 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 = 0,
#  nTimeSteps=1,

# 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,

# 1 year
# 365 days
# nTimeSteps=8760,

# 366 days
# nTimeSteps=8785,

# what const. did
#  nTimeSteps=8784,

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

 deltaTmom=3600.0,
 deltaTtracer=3600.0,
 deltaTClock =3600.0,
 cAdjFreq=0.,
 abEps=0.1,

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

#  dumpFreq = 21600.,
#  taveFreq = 21600.,
  
# Daily
# dumpFreq = 86400.,
# taveFreq = 86400.,
# adjDumpFreq = 86400.,
# chkptFreq = 86400.,
# pChekptFreq = 86400.,
# monitorFreq = 86400.,

# Weekly
# dumpFreq = 604800.,
# taveFreq = 604800.,
# adjDumpFreq = 604800.,
# chkptFreq  = 604800.,
# pChekptFreq = 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=.TRUE.,
 usingCurvilinearGrid=.FALSE.,

 delX=41*1.E0,
 delY=33*1.E0,
 delZ=10.,10.,15.,20.,20.,25.,35.,50.,75.,100.,150.,200.,275.,350.,415.,450.,500.,500.,500.,500.,500.,500.,500.,

 phimin=45.0,
 thetamin=279.0,
 &

/

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

# IAN
  bathyFile='LABSEA_bathy_fl_closed_boundary.bin',

# hydrogThetaFile = 'trial_t',
# hydrogSaltFile = 'trial_s',
# uVelInitFile = 'trial_u',
# vVelInitFile = 'trial_u',
#
#pSurfInitFile = 'trial_u',

 hydrogThetaFile = 'LABSEA_LEVITUS_1x1_ptmp_NEW_corK',
 hydrogSaltFile = 'LABSEA_LEVITUS_1x1_salt_NEW_corK',

#  hydrogThetaFile = 'ECCO_ITERATION_134_TBAR_2000_LS',
#  hydrogSaltFile  = 'ECCO_ITERATION_134_SBAR_2000_LS',
#  uVelInitFile    = 'ECCO_ITERATION_134_UBAR_2000_LS',
#  vVelInitFile    = 'ECCO_ITERATION_134_VBAR_2000_LS',
#  pSurfInitFile   = 'ECCO_ITERATION_134_PSBAR_2000_LS',
# END IAN

 /