Annotation of /MITgcm_contrib/heimbach/admtlm_setup/input_admtlm/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 15*20.,
 sRef = 15*35.,
 viscAr=1.E-3,
 viscAh=5.E5,
 diffKhT=0.0, 
 diffKrT=3.E-5,
 diffKhS=0.0,  
 diffKrS=3.E-5,
 rhonil=1035.,
 gravity=9.81,
#ph( new test
# eosType = 'JMD95Z',
 eosType = 'POLY3',
#ph)
 ivdc_kappa=100.,
 implicitDiffusion=.TRUE.,
#ph( new test
 vectorInvariantMomentum=.FALSE.,
 multiDimAdvection=.FALSE.,
 useRealFreshWaterFlux=.FALSE.,
 allowFreezing=.FALSE.,
 staggerTimeStep=.FALSE.,
 exactConserv = .FALSE.,
 useCDscheme=.TRUE.,
 tempAdvScheme=3,
 saltAdvScheme=3,
# bottomDragQuadratic=1.E-3,
# bottomDragLinear=1.E-3,
#ph)
# turn on looped cells
 hFacMin=.2,
 hFacMindr=50.,
# set precision of data files
 readBinaryPrec=32,
 &

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

# Time stepping parameters
 &PARM03
 tauCD =     321428.,
# nIter0 =      0,
# nTimeSteps = 120,
# 100 years of integration will yield a reasonable flow field
 startTime  = 0.,
# 1 year
 nTimesteps = 180,
# o 20 years
# endTime    = 622080000.,
# o 2 years
# endTime    = 62208000.,
# o 40 years
 deltaTmom      = 1200.0,
 deltaTtracer   = 172800.0,
 deltaTClock    = 172800.0,
 deltaTfreesurf = 172800.0,
# if you are using a version later than checkpoint45d on the main branch
# you can uncomment the following line and increase the time step 
# deltaTtracer and deltaTClock to 172800.0 as well to speed up the 
# asynchronous time stepping
 abEps = 0.1,
# pChkptFreq= 311040000.,
 taveFreq    = 2592000.,
 adjDumpFreq = 2592000.,
 monitorFreq = 31104000.,
# 2 months restoring timescale for temperature
 tauThetaClimRelax =  5184000.0,
# 6 months restoring timescale for salinity
 tauSaltClimRelax = 15552000.0,
 periodicExternalForcing=.TRUE.,
 externForcingPeriod=2592000.,
 externForcingCycle=31104000.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.TRUE.,
 delR= 50., 70., 100., 140., 190., 
       240., 290., 340., 390., 440., 
       490., 540., 590., 640., 690.,
 phiMin=-80.,
 dySpacing=4.,
 dxSpacing=4.,
 &

# Input datasets
 &PARM05
 bathyFile=      'bathymetry.bin',
 hydrogThetaFile='lev_t.bin',
 hydrogSaltFile= 'lev_s.bin',
 zonalWindFile=  'trenberth_taux.bin',
 meridWindFile=  'trenberth_tauy.bin',
 thetaClimFile=  'lev_sst.bin',
 saltClimFile=   'lev_sss.bin',
 surfQFile=      'ncep_qnet.bin',
# fresh water flux is turned off, uncomment next line to turn on
# (not recommened together with surface salinity restoring)
# EmPmRFile=      'ncep_emp.bin',
 &
1	heimbach	1.1	# ====================
2			# \| Model parameters \|
3			# ====================
4			#
5			# Continuous equation parameters
6			&PARM01
7			tRef = 15*20.,
8			sRef = 15*35.,
9			viscAr=1.E-3,
10			viscAh=5.E5,
11			diffKhT=0.0,
12			diffKrT=3.E-5,
13			diffKhS=0.0,
14			diffKrS=3.E-5,
15			rhonil=1035.,
16			gravity=9.81,
17			#ph( new test
18			# eosType = 'JMD95Z',
19			eosType = 'POLY3',
20			#ph)
21			ivdc_kappa=100.,
22			implicitDiffusion=.TRUE.,
23			#ph( new test
24			vectorInvariantMomentum=.FALSE.,
25			multiDimAdvection=.FALSE.,
26			useRealFreshWaterFlux=.FALSE.,
27			allowFreezing=.FALSE.,
28			staggerTimeStep=.FALSE.,
29			exactConserv = .FALSE.,
30			useCDscheme=.TRUE.,
31			tempAdvScheme=3,
32			saltAdvScheme=3,
33			# bottomDragQuadratic=1.E-3,
34			# bottomDragLinear=1.E-3,
35			#ph)
36			# turn on looped cells
37			hFacMin=.2,
38			hFacMindr=50.,
39			# set precision of data files
40			readBinaryPrec=32,
41			&
42
43			# Elliptic solver parameters
44			&PARM02
45			cg2dMaxIters=500,
46			cg2dTargetResidual=1.E-13,
47			&
48
49			# Time stepping parameters
50			&PARM03
51			tauCD = 321428.,
52			# nIter0 = 0,
53			# nTimeSteps = 120,
54			# 100 years of integration will yield a reasonable flow field
55			startTime = 0.,
56			# 1 year
57			nTimesteps = 180,
58			# o 20 years
59			# endTime = 622080000.,
60			# o 2 years
61			# endTime = 62208000.,
62			# o 40 years
63			deltaTmom = 1200.0,
64			deltaTtracer = 172800.0,
65			deltaTClock = 172800.0,
66			deltaTfreesurf = 172800.0,
67			# if you are using a version later than checkpoint45d on the main branch
68			# you can uncomment the following line and increase the time step
69			# deltaTtracer and deltaTClock to 172800.0 as well to speed up the
70			# asynchronous time stepping
71			abEps = 0.1,
72			# pChkptFreq= 311040000.,
73			taveFreq = 2592000.,
74			adjDumpFreq = 2592000.,
75			monitorFreq = 31104000.,
76			# 2 months restoring timescale for temperature
77			tauThetaClimRelax = 5184000.0,
78			# 6 months restoring timescale for salinity
79			tauSaltClimRelax = 15552000.0,
80			periodicExternalForcing=.TRUE.,
81			externForcingPeriod=2592000.,
82			externForcingCycle=31104000.,
83			&
84
85			# Gridding parameters
86			&PARM04
87			usingCartesianGrid=.FALSE.,
88			usingSphericalPolarGrid=.TRUE.,
89			delR= 50., 70., 100., 140., 190.,
90			240., 290., 340., 390., 440.,
91			490., 540., 590., 640., 690.,
92			phiMin=-80.,
93			dySpacing=4.,
94			dxSpacing=4.,
95			&
96
97			# Input datasets
98			&PARM05
99			bathyFile= 'bathymetry.bin',
100			hydrogThetaFile='lev_t.bin',
101			hydrogSaltFile= 'lev_s.bin',
102			zonalWindFile= 'trenberth_taux.bin',
103			meridWindFile= 'trenberth_tauy.bin',
104			thetaClimFile= 'lev_sst.bin',
105			saltClimFile= 'lev_sss.bin',
106			surfQFile= 'ncep_qnet.bin',
107			# fresh water flux is turned off, uncomment next line to turn on
108			# (not recommened together with surface salinity restoring)
109			# EmPmRFile= 'ncep_emp.bin',
110			&