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