1 |
# |
2 |
# ****************************** |
3 |
# Model parameters |
4 |
# Continuous equation parameters |
5 |
# ****************************** |
6 |
&PARM01 |
7 |
tRef= 24.0 , 23.0 , 22.0 , 21.0 , 20.0 , |
8 |
19.0 , 18.0 , 17.0 , 16.0 , 15.0 , |
9 |
14.0 , 13.0 , 12.0 , 11.0 , 10.0 , |
10 |
9.0 , 8.0 , 7.0 , 6.0, 5.0 , |
11 |
4.0 , 3.0 , 2.0 , |
12 |
|
13 |
sRef= 34.65, 34.75, 34.82, 34.87, 34.90, |
14 |
34.90, 34.86, 34.78, 34.69, 34.60, |
15 |
34.58, 34.62, 34.68, 34.72, 34.73, |
16 |
34.74, 34.73, 34.73, 34.72, 34.72, |
17 |
34.71, 34.70, 34.69, |
18 |
|
19 |
no_slip_sides=.false., |
20 |
no_slip_bottom=.false., |
21 |
|
22 |
# Harmonic Values From ECCO |
23 |
# viscAz=1.E-3, |
24 |
# viscAh=1.E4, |
25 |
|
26 |
# diffKhT=1.E2, |
27 |
# diffKzT=1.E-5, |
28 |
# diffKhS=1.E2, |
29 |
# diffKzS=1.E-5, |
30 |
|
31 |
# Baylor Suggestion #1 about an order of mag. less |
32 |
# and diffusion off. |
33 |
# actually ~ 6 times less in the horizontal.. |
34 |
|
35 |
viscAz=1.E-3, |
36 |
viscAh=1.5E3, |
37 |
|
38 |
# diffKhT=0, |
39 |
# diffKzT=0, |
40 |
# diffKhS=0, |
41 |
# diffKzS=0, |
42 |
|
43 |
# Baylor Suggestion #2 |
44 |
|
45 |
# BIHARMONIC |
46 |
#viscA4 = 0., |
47 |
#viscA4 = 1.39e-7, |
48 |
|
49 |
# the same damping on the smallest scales. |
50 |
# viscA4 = 3.225e11, |
51 |
|
52 |
# viscA4Grid = 2e-3, |
53 |
|
54 |
# a bit more than that to really get rid of that noise. |
55 |
viscA4Grid = 1e-1, |
56 |
|
57 |
diffK4T=0., |
58 |
diffK4S=0., |
59 |
|
60 |
# Advection Scheme 33 is itself a bit diffsive. |
61 |
tempAdvScheme=33, |
62 |
saltAdvScheme=33, |
63 |
|
64 |
StaggerTimeStep=.true. |
65 |
multiDimAdvection=.true. |
66 |
|
67 |
beta=1.E-11, |
68 |
tAlpha=2.E-4, |
69 |
sBeta =7.4E-4, |
70 |
gravity=9.81, |
71 |
gBaro=9.81, |
72 |
rigidLid=.FALSE., |
73 |
implicitFreeSurface=.true., |
74 |
eosType='JMD95Z', |
75 |
readBinaryPrec=32, |
76 |
writeBinaryPrec=32, |
77 |
saltStepping=.TRUE., |
78 |
tempStepping=.TRUE., |
79 |
momStepping=.TRUE., |
80 |
implicitDiffusion=.true., |
81 |
implicitViscosity=.true., |
82 |
globalFiles=.false., |
83 |
useSingleCpuIO=.TRUE., |
84 |
useCDscheme=.FALSE., |
85 |
exactConserv=.false., |
86 |
debugLevel=0, |
87 |
vectorInvariantMomentum=.TRUE., |
88 |
/ |
89 |
|
90 |
# ************************** |
91 |
# Elliptic solver parameters |
92 |
# ************************** |
93 |
&PARM02 |
94 |
|
95 |
# cg2dMaxIters=1000, |
96 |
cg2dMaxIters=100, |
97 |
# cg2dTargetResidual=1.E-11, |
98 |
cg2dTargetResidual=1.E-5, |
99 |
/ |
100 |
|
101 |
# ************************ |
102 |
# Time stepping parameters |
103 |
# ************************ |
104 |
&PARM03 |
105 |
# calendarDumps=.TRUE., |
106 |
calendarDumps=.FALSE., |
107 |
|
108 |
# |
109 |
# nIter0 = 1, |
110 |
|
111 |
nIter0 = 0, |
112 |
# nTimeSteps = 1, |
113 |
# one week = 168 |
114 |
# nTimeSteps=168, |
115 |
|
116 |
# just 12 time stpes |
117 |
# nTimeSteps = 12, |
118 |
|
119 |
|
120 |
# 1 month 30.5*24, |
121 |
# nTimeSteps=732, |
122 |
# 1 month 31*24, |
123 |
# nTimeSteps=750 |
124 |
|
125 |
# 2 months 30.5*24*2 |
126 |
# nTimeSteps=1464, |
127 |
|
128 |
#nTimeSteps = 7300, |
129 |
|
130 |
# about 10 months. |
131 |
# nTimeSteps=7300, |
132 |
|
133 |
# 1 year |
134 |
# 365 days |
135 |
nTimeSteps=8760, |
136 |
|
137 |
# 366 days |
138 |
# nTimeSteps=8785, |
139 |
|
140 |
# what const. did |
141 |
# nTimeSteps=8784, |
142 |
|
143 |
# 5 years |
144 |
# nTimeSteps=43920, |
145 |
|
146 |
# 2 years |
147 |
# nTimeSteps=17440, |
148 |
# 11 years |
149 |
# nTimeSteps=96398, |
150 |
# 12 years |
151 |
# nTimeSteps=105120, |
152 |
# 13 years |
153 |
# nTimeSteps=113880, |
154 |
# |
155 |
forcing_In_AB = .FALSE., |
156 |
|
157 |
# deltaTtracer=600.0, |
158 |
# deltaTClock =600.0, |
159 |
# deltaTmom = 600.0, |
160 |
|
161 |
deltaTmom = 3600.0, |
162 |
deltaTtracer= 3600.0, |
163 |
deltaTClock = 3600.0, |
164 |
|
165 |
# deltaTmom = 200.0, |
166 |
# deltaTtracer=200.0, |
167 |
# deltaTClock =200.0, |
168 |
|
169 |
|
170 |
cAdjFreq=0., |
171 |
abEps=0.1, |
172 |
|
173 |
# IAN |
174 |
# every time step |
175 |
# dumpFreq = 1., |
176 |
|
177 |
# dumpFreq = 3600., |
178 |
# taveFreq = 3600., |
179 |
# adjDumpFreq = 3600., |
180 |
# chkptFreq = 3600., |
181 |
# pChkptFreq = 3600., |
182 |
# monitorFreq = 3600., |
183 |
|
184 |
# Daily |
185 |
|
186 |
# dumpFreq = 86400., |
187 |
# dumpFreq = 1., |
188 |
# monitorFreq=1., |
189 |
|
190 |
# taveFreq = 86400., |
191 |
# adjDumpFreq = 86400., |
192 |
# chkptFreq = 86400., |
193 |
# pChekptFreq = 86400., |
194 |
# monitorFreq = 86400., |
195 |
|
196 |
# Weekly |
197 |
# dumpFreq = 604800., |
198 |
# taveFreq = 604800., |
199 |
# adjDumpFreq = 604800., |
200 |
# chkptFreq = 604800., |
201 |
# pChkptFreq = 604800., |
202 |
|
203 |
# two weeks |
204 |
# dumpFreq = 54900., |
205 |
# taveFrew = 54900., |
206 |
|
207 |
# monthly. 30.5 days |
208 |
dumpFreq = 2635200., |
209 |
taveFreq = 2635200., |
210 |
|
211 |
# monthly 31 days |
212 |
# dumpFreq = 2678400., |
213 |
# TAVeFreq = 2678400., |
214 |
|
215 |
# ENDIAN |
216 |
|
217 |
# year? |
218 |
# monitorFreq = 31536000.0, |
219 |
/ |
220 |
|
221 |
# ******************* |
222 |
# Gridding parameters |
223 |
# ******************* |
224 |
&PARM04 |
225 |
usingCartesianGrid=.FALSE., |
226 |
usingSphericalPolarGrid=.FALSE., |
227 |
usingCurvilinearGrid=.TRUE., |
228 |
delZ=10.,10.,15.,20.,20.,25.,35.,50.,75.,100.,150.,200.,275.,350.,415.,450.,500.,500.,500.,500.,500.,500.,500., |
229 |
horizGridFile = 'a', |
230 |
& |
231 |
/ |
232 |
|
233 |
# ********** |
234 |
# Data Files |
235 |
# ********** |
236 |
&PARM05 |
237 |
|
238 |
# IAN |
239 |
|
240 |
# bathyFile='ib_closed.bin', |
241 |
bathyFile='ib_open_15.bin', |
242 |
|
243 |
#hydrogSaltFile = 'salt_init.from_pickup', |
244 |
#hydrogThetaFile = 'theta_init.from_pickup', |
245 |
|
246 |
#uVelInitFile = 'u_init.from_pickup', |
247 |
#vVelInitFile = 'v_init.from_pickup', |
248 |
#pSurfInitFile = 'eta_init.from_pickup', |
249 |
|
250 |
# 22nd Jan |
251 |
hydrogSaltFile = 'salt_init.ecco_month_9', |
252 |
hydrogThetaFile = 'theta_init.ecco_month_9', |
253 |
|
254 |
uVelInitFile = 'u_init.ecco_month_9', |
255 |
vVelInitFile = 'v_init.ecco_month_9', |
256 |
|
257 |
|
258 |
# END IAN |
259 |
|
260 |
/ |