global_ocean.90x40x15/input/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 15*20.,
 sRef = 15*35.,
 viscAr=1.E-3,
 viscAh=5.E5,
 viscA4=1.E14,
 diffKhT=0.,
 diffKrT=3.E-5,
 diffKhS=0.,
 diffKrS=3.E-5,
 rhoConst=1035.,
 rhoConstFresh=1000.,
 gravity=9.81,
 eosType='JMD95P',
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 allowFreezing=.TRUE.,
 exactConserv=.TRUE.,
 select_rStar=2,
 nonlinFreeSurf=4,
 hFacInf=0.2,
 hFacSup=2.0,
 useRealFreshWaterFlux=.TRUE.,
 useCDscheme=.TRUE.,
 quasiHydrostatic=.TRUE.,
 useNHMTerms=.TRUE.,
# turn on looped cells
 hFacMin=.05,
 hFacMindr=50.,
#- to check that it has no impact:
 doResetHFactors=.TRUE.,
# set precision of data files
 readBinaryPrec=32,
#writeBinaryPrec=64,
#debugMode=.TRUE.,
 &

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

# Time stepping parameters
 &PARM03
 nIter0=36000,
 nTimeSteps=10,
# 100 years of integration will yield a reasonable flow field
# startTime =          0.,
# endTime   = 3110400000.,
 deltaTmom = 1800.,
 tauCD =     321428.,
 deltaTtracer= 86400.,
 deltaTClock = 86400.,
 deltaTfreesurf= 86400.,
 abEps = 0.1,
 pChkPtFreq= 311040000.,
 dumpFreq=   311040000.,
 dumpFreq=   864000.,
#monitorFreq=31104000.,
 monitorFreq=1.,
# 2 months restoring timescale for temperature
 tauThetaClimRelax=  5184000.,
# 6 months restoring timescale for salinity
 tauSaltClimRelax = 15552000.,
 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',
 the_run_name=   'global_ocean.90x40x15',
 EmPmRFile=      'ncep_emp.bin',
 &
1	adcroft	1.1	# ====================
2			# \| Model parameters \|
3			# ====================
4			#
5			# Continuous equation parameters
6			&PARM01
7			tRef = 15*20.,
8			sRef = 15*35.,
9	mlosch	1.2	viscAr=1.E-3,
10	adcroft	1.1	viscAh=5.E5,
11	jmc	1.18	viscA4=1.E14,
12	jmc	1.12	diffKhT=0.,
13	mlosch	1.2	diffKrT=3.E-5,
14	jmc	1.12	diffKhS=0.,
15	mlosch	1.2	diffKrS=3.E-5,
16	jmc	1.15	rhoConst=1035.,
17	jmc	1.12	rhoConstFresh=1000.,
18	adcroft	1.1	gravity=9.81,
19	jmc	1.13	eosType='JMD95P',
20	jmc	1.17	ivdc_kappa=10.,
21	adcroft	1.1	implicitDiffusion=.TRUE.,
22	jmc	1.12	allowFreezing=.TRUE.,
23			exactConserv=.TRUE.,
24	jmc	1.17	select_rStar=2,
25			nonlinFreeSurf=4,
26			hFacInf=0.2,
27			hFacSup=2.0,
28	mlosch	1.2	useRealFreshWaterFlux=.TRUE.,
29	jmc	1.4	useCDscheme=.TRUE.,
30	jmc	1.12	quasiHydrostatic=.TRUE.,
31	jmc	1.5	useNHMTerms=.TRUE.,
32	mlosch	1.2	# turn on looped cells
33			hFacMin=.05,
34			hFacMindr=50.,
35	jmc	1.19	#- to check that it has no impact:
36			doResetHFactors=.TRUE.,
37	mlosch	1.2	# set precision of data files
38			readBinaryPrec=32,
39	jmc	1.17	#writeBinaryPrec=64,
40			#debugMode=.TRUE.,
41	adcroft	1.1	&
42
43			# Elliptic solver parameters
44			&PARM02
45	mlosch	1.2	cg2dMaxIters=500,
46			cg2dTargetResidual=1.E-13,
47	adcroft	1.1	&
48
49			# Time stepping parameters
50			&PARM03
51	jmc	1.12	nIter0=36000,
52			nTimeSteps=10,
53	mlosch	1.2	# 100 years of integration will yield a reasonable flow field
54	jmc	1.12	# startTime = 0.,
55			# endTime = 3110400000.,
56			deltaTmom = 1800.,
57	adcroft	1.1	tauCD = 321428.,
58	jmc	1.12	deltaTtracer= 86400.,
59			deltaTClock = 86400.,
60			deltaTfreesurf= 86400.,
61	adcroft	1.1	abEps = 0.1,
62	jmc	1.17	pChkPtFreq= 311040000.,
63	mlosch	1.2	dumpFreq= 311040000.,
64			dumpFreq= 864000.,
65	adcroft	1.3	#monitorFreq=31104000.,
66			monitorFreq=1.,
67	mlosch	1.2	# 2 months restoring timescale for temperature
68	jmc	1.12	tauThetaClimRelax= 5184000.,
69	mlosch	1.2	# 6 months restoring timescale for salinity
70	jmc	1.12	tauSaltClimRelax = 15552000.,
71	adcroft	1.1	periodicExternalForcing=.TRUE.,
72			externForcingPeriod=2592000.,
73			externForcingCycle=31104000.,
74			&
75
76			# Gridding parameters
77			&PARM04
78			usingSphericalPolarGrid=.TRUE.,
79	jmc	1.12	delR= 50., 70., 100., 140., 190.,
80			240., 290., 340., 390., 440.,
81			490., 540., 590., 640., 690.,
82	jmc	1.16	ygOrigin=-80.,
83	mlosch	1.2	dySpacing=4.,
84			dxSpacing=4.,
85	adcroft	1.1	&
86
87			# Input datasets
88			&PARM05
89	mlosch	1.2	bathyFile= 'bathymetry.bin',
90			hydrogThetaFile='lev_t.bin',
91			hydrogSaltFile= 'lev_s.bin',
92			zonalWindFile= 'trenberth_taux.bin',
93			meridWindFile= 'trenberth_tauy.bin',
94			thetaClimFile= 'lev_sst.bin',
95			saltClimFile= 'lev_sss.bin',
96			surfQFile= 'ncep_qnet.bin',
97	edhill	1.7	the_run_name= 'global_ocean.90x40x15',
98	jmc	1.12	EmPmRFile= 'ncep_emp.bin',
99	adcroft	1.1	&