global_ocean.90x40x15/input_ad/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,
 eosType = 'POLY3',
 ivdc_kappa=100.,
 implicitDiffusion=.TRUE.,
 allowFreezing=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
# turn on looped cells
 hFacMin=.05,
 hFacMindr=50.,
# set precision of data files
 readBinaryPrec=32,
 &

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

# Time stepping parameters
 &PARM03
 nIter0 =      0,
 nTimeSteps = 60,
#ph should perform at least 60 timesteps (30 days) for gradient
#ph checks to allow signal to propagate away from cost function.
#
# 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.,
 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
 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.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	heimbach	1.2	viscAh=5.E5,
11	heimbach	1.3	diffKhT=0.0,
12			diffKrT=3.E-5,
13			diffKhS=0.0,
14			diffKrS=3.E-5,
15			rhonil=1035.,
16	heimbach	1.2	gravity=9.81,
17	heimbach	1.3	eosType = 'POLY3',
18	heimbach	1.2	ivdc_kappa=100.,
19	heimbach	1.3	implicitDiffusion=.TRUE.,
20	heimbach	1.2	allowFreezing=.TRUE.,
21	heimbach	1.3	useRealFreshWaterFlux=.TRUE.,
22			# turn on looped cells
23			hFacMin=.05,
24			hFacMindr=50.,
25			# set precision of data files
26	heimbach	1.2	readBinaryPrec=32,
27			&
28
29			# Elliptic solver parameters
30			&PARM02
31	heimbach	1.3	cg2dMaxIters=500,
32			cg2dTargetResidual=1.E-13,
33	heimbach	1.2	&
34
35			# Time stepping parameters
36			&PARM03
37	heimbach	1.3	nIter0 = 0,
38			nTimeSteps = 60,
39			#ph should perform at least 60 timesteps (30 days) for gradient
40			#ph checks to allow signal to propagate away from cost function.
41	heimbach	1.2	#
42	heimbach	1.3	# 100 years of integration will yield a reasonable flow field
43			# startTime = 0.,
44			# endTime = 3110400000.,
45			deltaTmom = 1200.0,
46			tauCD = 321428.,
47			deltaTtracer= 43200.0,
48			deltaTClock = 43200.0,
49			# if you are using a version later than checkpoint45d on the main branch
50			# you can uncomment the following line and increase the time step
51			# deltaTtracer and deltaTClock to 172800.0 as well to speed up the
52			# asynchronous time stepping
53			# deltaTfreesurf = 172800.0,
54			abEps = 0.1,
55			pChkptFreq= 311040000.,
56			dumpFreq= 311040000.,
57			dumpFreq= 864000.,
58			taveFreq= 311040000.,
59			taveFreq= 864000.,
60			#monitorFreq=31104000.,
61			monitorFreq=1.,
62			# 2 months restoring timescale for temperature
63			tauThetaClimRelax = 5184000.0,
64			# 6 months restoring timescale for salinity
65			tauSaltClimRelax = 15552000.0,
66			periodicExternalForcing=.TRUE.,
67			externForcingPeriod=2592000.,
68			externForcingCycle=31104000.,
69	heimbach	1.2	&
70
71			# Gridding parameters
72			&PARM04
73			usingCartesianGrid=.FALSE.,
74			usingSphericalPolarGrid=.TRUE.,
75	heimbach	1.3	delR= 50., 70., 100., 140., 190.,
76			240., 290., 340., 390., 440.,
77			490., 540., 590., 640., 690.,
78			phiMin=-80.,
79			dySpacing=4.,
80			dxSpacing=4.,
81	heimbach	1.2	&
82
83			# Input datasets
84			&PARM05
85	heimbach	1.3	bathyFile= 'bathymetry.bin',
86			hydrogThetaFile='lev_t.bin',
87			hydrogSaltFile= 'lev_s.bin',
88			zonalWindFile= 'trenberth_taux.bin',
89			meridWindFile= 'trenberth_tauy.bin',
90			thetaClimFile= 'lev_sst.bin',
91			saltClimFile= 'lev_sss.bin',
92			surfQFile= 'ncep_qnet.bin',
93			# fresh water flux is turned off, uncomment next line to turn on
94			# (not recommened together with surface salinity restoring)
95			# EmPmRFile= 'ncep_emp.bin',
96	heimbach	1.2	&