Annotation of /MITgcm_contrib/gael/verification/ECCO_v4_r2/input/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef               = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13.,
                      12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1.,
 sRef               = 50*34.5,
 no_slip_sides  = .TRUE.,
 no_slip_bottom = .TRUE.,
#
 viscAr=0.5E-4,
#
 viscAh=1.E0,
 viscAhGrid=2.E-2,
# viscAh=2.0e4,
#
 diffKhT=1.E1,
 diffKrT=1.E-5,
 diffKhS=1.E1,
 diffKrS=1.E-5,
#
### diffKrBL79surf=0.1E-4,
### diffKrBL79deep=1.0E-4,
 bottomDragQuadratic = 0.001,
#when using ggl90
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
# balanceThetaClimRelax=.TRUE.,
 balanceSaltClimRelax=.TRUE.,
# balanceEmPmR=.TRUE.,
# balanceQnet=.TRUE.,
 allowFreezing=.FALSE.,
### hFacInf=0.2,
### hFacSup=2.0,
 hFacMin=.2,
 hFacMinDr=5.,
 select_rStar=2,
 nonlinFreeSurf=4,
 gravity=9.81,
 rhonil=1029.,
 rhoConst=1029.,
 rhoConstFresh=1000.,
 convertFW2Salt=-1.,
 eosType='JMD95Z',
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 useSingleCpuIO=.TRUE.,
 tempAdvScheme=30,
 saltAdvScheme=30,
 tempVertAdvScheme=3,
 saltVertAdvScheme=3,
 tempImplVertAdv=.TRUE.,
 saltImplVertAdv=.TRUE.,
 staggerTimeStep=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
#when using the cd scheme:
# useCDscheme=.TRUE.,
 useJamartWetPoints=.TRUE.,
 readBinaryPrec=32,
 writeBinaryPrec=32,
 debugLevel=1,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=300,
 cg2dTargetResWunit=1.E-12,
 &

# Time stepping parameters
 &PARM03
 nIter0=1,
#2 lev2 for testing:
#nTimeSteps=8,
#20y:
 nTimeSteps=175295,
#60 years
#nTimeSteps=525985,
#
 forcing_In_AB=.FALSE.,
 momDissip_In_AB=.FALSE.,
#when using the cd scheme:
# epsAB_CD = 0.25,
# tauCD=172800.0,
 deltaTmom   =3600.,
 deltaTtracer=3600.,
 deltaTfreesurf=3600.,
 deltaTClock =3600.,
#when using ab2:
# abEps = 0.1,
#when using ab3:
 doAB_onGtGs=.FALSE.,
 alph_AB=0.5,
 beta_AB=0.281105,
#
 pChkptFreq  =31536000.0,
 chkptFreq   =31536000.0,
# taveFreq    =31536000.0,
# dumpFreq    =31536000.0,
# monitorFreq = 7200.0,
 monitorFreq = 6307200.0,
 dumpInitAndLast = .TRUE.,
 adjMonitorFreq = 864000.0,
 pickupStrictlyMatch=.FALSE.,
 &

# Gridding parameters
 &PARM04
 usingCurvilinearGrid=.TRUE.,
 delR = 
     10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
     10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04, 19.82, 24.85,
     31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
     93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
     139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
     341.50,364.50,387.50,410.50,433.50,456.50,
 &

# Input datasets
 &PARM05
 adTapeDir='tapes',
 bathyFile      ='bathy.bin',
 hydrogThetaFile='T_OWPv1_M_eccollc_90x50.bin',
 hydrogSaltFile ='S_OWPv1_M_eccollc_90x50.bin',
 viscA4Dfile    ='viscA4Dfld.bin',
 viscA4Zfile    ='viscA4Zfld.bin',
 geothermalFile='geothermalFlux.bin',
#
 &
1	gforget	1.1	# ====================
2			# \| Model parameters \|
3			# ====================
4			#
5			# Continuous equation parameters
6			&PARM01
7			tRef = 323.,322.,21.,220.,19.,218.,17.,2*16.,15.,14.,13.,
8			12.,11.,29.,8.,7.,26.,25.,34.,33.,42.,12*1.,
9			sRef = 50*34.5,
10			no_slip_sides = .TRUE.,
11			no_slip_bottom = .TRUE.,
12			#
13	gforget	1.6	viscAr=0.5E-4,
14	gforget	1.1	#
15			viscAh=1.E0,
16			viscAhGrid=2.E-2,
17			# viscAh=2.0e4,
18			#
19			diffKhT=1.E1,
20			diffKrT=1.E-5,
21			diffKhS=1.E1,
22			diffKrS=1.E-5,
23			#
24			### diffKrBL79surf=0.1E-4,
25			### diffKrBL79deep=1.0E-4,
26			bottomDragQuadratic = 0.001,
27			#when using ggl90
28			ivdc_kappa=10.,
29			implicitDiffusion=.TRUE.,
30			implicitViscosity=.TRUE.,
31			useRealFreshWaterFlux=.TRUE.,
32			# balanceThetaClimRelax=.TRUE.,
33			balanceSaltClimRelax=.TRUE.,
34	gforget	1.6	# balanceEmPmR=.TRUE.,
35	gforget	1.1	# balanceQnet=.TRUE.,
36			allowFreezing=.FALSE.,
37			### hFacInf=0.2,
38			### hFacSup=2.0,
39			hFacMin=.2,
40			hFacMinDr=5.,
41			select_rStar=2,
42	gforget	1.3	nonlinFreeSurf=4,
43	gforget	1.1	gravity=9.81,
44			rhonil=1029.,
45			rhoConst=1029.,
46			rhoConstFresh=1000.,
47			convertFW2Salt=-1.,
48			eosType='JMD95Z',
49			implicitFreeSurface=.TRUE.,
50			exactConserv=.TRUE.,
51			useSingleCpuIO=.TRUE.,
52			tempAdvScheme=30,
53			saltAdvScheme=30,
54			tempVertAdvScheme=3,
55			saltVertAdvScheme=3,
56			tempImplVertAdv=.TRUE.,
57			saltImplVertAdv=.TRUE.,
58			staggerTimeStep=.TRUE.,
59			vectorInvariantMomentum=.TRUE.,
60			#when using the cd scheme:
61			# useCDscheme=.TRUE.,
62			useJamartWetPoints=.TRUE.,
63			readBinaryPrec=32,
64			writeBinaryPrec=32,
65			debugLevel=1,
66	gforget	1.5	&
67	gforget	1.1
68			# Elliptic solver parameters
69			&PARM02
70			cg2dMaxIters=300,
71			cg2dTargetResWunit=1.E-12,
72	gforget	1.5	&
73	gforget	1.1
74			# Time stepping parameters
75			&PARM03
76	gforget	1.2	nIter0=1,
77	gforget	1.1	#2 lev2 for testing:
78	gforget	1.6	#nTimeSteps=8,
79			#20y:
80			nTimeSteps=175295,
81	gforget	1.1	#60 years
82			#nTimeSteps=525985,
83			#
84			forcing_In_AB=.FALSE.,
85			momDissip_In_AB=.FALSE.,
86			#when using the cd scheme:
87			# epsAB_CD = 0.25,
88			# tauCD=172800.0,
89			deltaTmom =3600.,
90			deltaTtracer=3600.,
91			deltaTfreesurf=3600.,
92			deltaTClock =3600.,
93			#when using ab2:
94			# abEps = 0.1,
95			#when using ab3:
96			doAB_onGtGs=.FALSE.,
97			alph_AB=0.5,
98			beta_AB=0.281105,
99			#
100			pChkptFreq =31536000.0,
101			chkptFreq =31536000.0,
102			# taveFreq =31536000.0,
103			# dumpFreq =31536000.0,
104	gforget	1.6	# monitorFreq = 7200.0,
105			monitorFreq = 6307200.0,
106	gforget	1.1	dumpInitAndLast = .TRUE.,
107			adjMonitorFreq = 864000.0,
108			pickupStrictlyMatch=.FALSE.,
109	gforget	1.5	&
110	gforget	1.1
111			# Gridding parameters
112			&PARM04
113			usingCurvilinearGrid=.TRUE.,
114			delR =
115			10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
116			10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04, 19.82, 24.85,
117			31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
118			93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
119			139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
120			341.50,364.50,387.50,410.50,433.50,456.50,
121	gforget	1.5	&
122	gforget	1.1
123			# Input datasets
124			&PARM05
125			adTapeDir='tapes',
126	gforget	1.7	bathyFile ='bathy.bin',
127	gforget	1.6	hydrogThetaFile='T_OWPv1_M_eccollc_90x50.bin',
128			hydrogSaltFile ='S_OWPv1_M_eccollc_90x50.bin',
129	gforget	1.7	viscA4Dfile ='viscA4Dfld.bin',
130			viscA4Zfile ='viscA4Zfld.bin',
131	gforget	1.8	geothermalFile='geothermalFlux.bin',
132	gforget	1.1	#
133	gforget	1.5	&