Contents of /MITgcm_contrib/ifenty/ECCO_v4/input_ad/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.,
 adjDumpFreq = 31536000.0,
 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_eccollc_90x50.bin',
 bathyFile      ='bathy_eccollc_90x50_min2pts.bin',
 hydrogThetaFile='T_OWPv1_M_eccollc_90x50.bin',
 hydrogSaltFile ='S_OWPv1_M_eccollc_90x50.bin',
 viscA4Dfile    ='viscA4Dfld_eccollc_90x50.bin',
 viscA4Zfile    ='viscA4Zfld_eccollc_90x50.bin',
#
 /
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	viscAr=0.5E-4,
14	#
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	# balanceEmPmR=.TRUE.,
35	# balanceQnet=.TRUE.,
36	allowFreezing=.FALSE.,
37	### hFacInf=0.2,
38	### hFacSup=2.0,
39	hFacMin=.2,
40	hFacMinDr=5.,
41	select_rStar=2,
42	nonlinFreeSurf=4,
43	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	/
67
68	# Elliptic solver parameters
69	&PARM02
70	cg2dMaxIters=300,
71	cg2dTargetResWunit=1.E-12,
72	/
73
74	# Time stepping parameters
75	&PARM03
76	nIter0=1,
77	#2 lev2 for testing:
78	#nTimeSteps=8,
79	#20y:
80	nTimeSteps=175295,
81	#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	# monitorFreq = 7200.0,
105	monitorFreq = 6307200.0,
106	dumpInitAndLast = .TRUE.,
107	adjDumpFreq = 31536000.0,
108	adjMonitorFreq = 864000.0,
109	pickupStrictlyMatch=.FALSE.,
110	/
111
112	# Gridding parameters
113	&PARM04
114	usingCurvilinearGrid=.TRUE.,
115	delR =
116	10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
117	10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04, 19.82, 24.85,
118	31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
119	93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
120	139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
121	341.50,364.50,387.50,410.50,433.50,456.50,
122	/
123
124	# Input datasets
125	&PARM05
126	adTapeDir='tapes',
127	#bathyFile ='bathy_eccollc_90x50.bin',
128	bathyFile ='bathy_eccollc_90x50_min2pts.bin',
129	hydrogThetaFile='T_OWPv1_M_eccollc_90x50.bin',
130	hydrogSaltFile ='S_OWPv1_M_eccollc_90x50.bin',
131	viscA4Dfile ='viscA4Dfld_eccollc_90x50.bin',
132	viscA4Zfile ='viscA4Zfld_eccollc_90x50.bin',
133	#
134	/