global_ocean.90x40x15/input_oad/data

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

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