global_ocean.90x40x15/input_ad/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 15*20.,
 sRef = 15*35.,
 viscAr=1.E-3,
 diffKhT=0.0,
 diffKrT=3.E-5,
 diffKhS=0.0,
 diffKrS=3.E-5,
#ph(
#ph Many of the following params. are
#ph just for AD testing; not necesarily
#ph meaningful in a 4x4 deg. setting.
 viscC4Leith=1.5,
 viscC4Leithd=1.5,
 viscA4GridMax=0.5,
 useAreaViscLength=.TRUE.,
 bottomDragQuadratic = 0.0021,
#ph)
 rhonil=1035.,
 rotationPeriod=86400.,
 gravity=9.81,
 eosType = 'JMD95Z',
 ivdc_kappa=100.,
 implicitDiffusion=.TRUE.,
 allowFreezing=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
 useCDscheme=.TRUE.,
 exactConserv=.TRUE.
# turn on looped cells
 hFacMin=.05,
 hFacMindr=50.,
# set precision of data files
 readBinaryPrec=32,
 debugLevel=1,
#ph(
 vectorInvariantMomentum=.TRUE.,
 hFacInf=0.2,
 hFacSup=2.0,
 nonlinFreeSurf=2,
 exactConserv=.TRUE.
 useRealFreshWaterFlux=.TRUE.,
 staggerTimeStep=.TRUE.,
 multiDimAdvection=.TRUE,
 tempAdvScheme=30,
 saltAdvScheme=30,
 tempVertAdvScheme=30,
 saltVertAdvScheme=30,
 tempImplVertAdv=.TRUE.,
 saltImplVertAdv=.TRUE.,
 inAdExact=.TRUE.,
#ph)
 &

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

# Time stepping parameters
 &PARM03
 nIter0 =      0,
 nTimeSteps = 10,
# 100 years of integration will yield a reasonable flow field
# startTime  =          0.,
# endTime    = 3110400000.,
 deltaTmom = 1200.0,
 tauCD =     321428.,
 deltaTtracer= 43200.0,
 deltaTClock = 43200.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
# deltaTfreesurf = 172800.0,
 abEps = 0.1,
 pChkptFreq= 311040000.,
 dumpFreq=   311040000.,
 dumpFreq=   864000.,
 adjDumpFreq= 432000.,
 adjMonitorFreq=1.,
# taveFreq=   311040000.,
# taveFreq=   864000.,
#monitorFreq=31104000.,
 monitorFreq=1.,
# 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
 usingSphericalPolarGrid=.TRUE.,
 delR= 50., 70., 100., 140., 190.,
       240., 290., 340., 390., 440.,
       490., 540., 590., 640., 690.,
 ygOrigin=-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.2	# ====================
2			# \| Model parameters \|
3			# ====================
4			#
5			# Continuous equation parameters
6			&PARM01
7	heimbach	1.3	tRef = 15*20.,
8			sRef = 15*35.,
9			viscAr=1.E-3,
10	jmc	1.7	diffKhT=0.0,
11	heimbach	1.3	diffKrT=3.E-5,
12	jmc	1.7	diffKhS=0.0,
13	heimbach	1.3	diffKrS=3.E-5,
14	heimbach	1.12	#ph(
15			#ph Many of the following params. are
16			#ph just for AD testing; not necesarily
17			#ph meaningful in a 4x4 deg. setting.
18			viscC4Leith=1.5,
19			viscC4Leithd=1.5,
20			viscA4GridMax=0.5,
21			useAreaViscLength=.TRUE.,
22			bottomDragQuadratic = 0.0021,
23			#ph)
24	heimbach	1.3	rhonil=1035.,
25	jmc	1.7	rotationPeriod=86400.,
26	heimbach	1.2	gravity=9.81,
27	heimbach	1.9	eosType = 'JMD95Z',
28	heimbach	1.2	ivdc_kappa=100.,
29	heimbach	1.3	implicitDiffusion=.TRUE.,
30	heimbach	1.2	allowFreezing=.TRUE.,
31	heimbach	1.3	useRealFreshWaterFlux=.TRUE.,
32	jmc	1.5	useCDscheme=.TRUE.,
33	heimbach	1.12	exactConserv=.TRUE.
34	heimbach	1.3	# turn on looped cells
35			hFacMin=.05,
36			hFacMindr=50.,
37			# set precision of data files
38	heimbach	1.2	readBinaryPrec=32,
39	heimbach	1.11	debugLevel=1,
40	heimbach	1.13	#ph(
41			vectorInvariantMomentum=.TRUE.,
42			hFacInf=0.2,
43			hFacSup=2.0,
44			nonlinFreeSurf=2,
45			exactConserv=.TRUE.
46			useRealFreshWaterFlux=.TRUE.,
47			staggerTimeStep=.TRUE.,
48			multiDimAdvection=.TRUE,
49			tempAdvScheme=30,
50			saltAdvScheme=30,
51	heimbach	1.15	tempVertAdvScheme=30,
52			saltVertAdvScheme=30,
53			tempImplVertAdv=.TRUE.,
54			saltImplVertAdv=.TRUE.,
55	heimbach	1.13	inAdExact=.TRUE.,
56			#ph)
57	heimbach	1.2	&
58
59			# Elliptic solver parameters
60			&PARM02
61	heimbach	1.3	cg2dMaxIters=500,
62			cg2dTargetResidual=1.E-13,
63	heimbach	1.2	&
64
65			# Time stepping parameters
66			&PARM03
67	heimbach	1.3	nIter0 = 0,
68	heimbach	1.9	nTimeSteps = 10,
69	heimbach	1.3	# 100 years of integration will yield a reasonable flow field
70			# startTime = 0.,
71			# endTime = 3110400000.,
72			deltaTmom = 1200.0,
73			tauCD = 321428.,
74			deltaTtracer= 43200.0,
75			deltaTClock = 43200.0,
76			# if you are using a version later than checkpoint45d on the main branch
77	jmc	1.7	# you can uncomment the following line and increase the time step
78			# deltaTtracer and deltaTClock to 172800.0 as well to speed up the
79	heimbach	1.3	# asynchronous time stepping
80			# deltaTfreesurf = 172800.0,
81			abEps = 0.1,
82			pChkptFreq= 311040000.,
83			dumpFreq= 311040000.,
84			dumpFreq= 864000.,
85	heimbach	1.6	adjDumpFreq= 432000.,
86	heimbach	1.8	adjMonitorFreq=1.,
87	heimbach	1.6	# taveFreq= 311040000.,
88			# taveFreq= 864000.,
89	heimbach	1.3	#monitorFreq=31104000.,
90			monitorFreq=1.,
91			# 2 months restoring timescale for temperature
92			tauThetaClimRelax = 5184000.0,
93			# 6 months restoring timescale for salinity
94			tauSaltClimRelax = 15552000.0,
95			periodicExternalForcing=.TRUE.,
96			externForcingPeriod=2592000.,
97			externForcingCycle=31104000.,
98	heimbach	1.2	&
99
100			# Gridding parameters
101			&PARM04
102			usingSphericalPolarGrid=.TRUE.,
103	jmc	1.7	delR= 50., 70., 100., 140., 190.,
104			240., 290., 340., 390., 440.,
105	heimbach	1.3	490., 540., 590., 640., 690.,
106	jmc	1.14	ygOrigin=-80.,
107	heimbach	1.3	dySpacing=4.,
108			dxSpacing=4.,
109	heimbach	1.2	&
110
111			# Input datasets
112			&PARM05
113	heimbach	1.3	bathyFile= 'bathymetry.bin',
114			hydrogThetaFile='lev_t.bin',
115			hydrogSaltFile= 'lev_s.bin',
116			zonalWindFile= 'trenberth_taux.bin',
117			meridWindFile= 'trenberth_tauy.bin',
118			thetaClimFile= 'lev_sst.bin',
119			saltClimFile= 'lev_sss.bin',
120			surfQFile= 'ncep_qnet.bin',
121			# fresh water flux is turned off, uncomment next line to turn on
122			# (not recommened together with surface salinity restoring)
123			# EmPmRFile= 'ncep_emp.bin',
124	heimbach	1.2	&