-
Notifications
You must be signed in to change notification settings - Fork 6
/
ww3_grid_namelist.inp
381 lines (381 loc) · 21.7 KB
/
ww3_grid_namelist.inp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
$ Start of namelist input section ------------------------------------ $
$ Starting with WAVEWATCH III version 2.00, the tunable parameters
$ for source terms, propagation schemes, and numerics are read using
$ namelists. Any namelist found in the folowing sections up to the
$ end-of-section identifier string (see below) is temporarily written
$ to ww3_grid.scratch, and read from there if necessary. Namelists
$ not needed for the given switch settings will be skipped
$ automatically, and the order of the namelists is immaterial.
$ As an example, namelist input to change SWELLF and ZWND in the
$ Tolman and Chalikov input would be
$
$ &SIN2 SWELLF = 0.1, ZWND = 15. /
$
$ Define constants in source terms ----------------------------------- $
$
$ Stresses - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
$ TC 1996 with cap : Namelist FLX3
$ CDMAX : Maximum allowed CD (cap)
$ CTYPE : Cap type :
$ 0: Discontinuous (default).
$ 1: Hyperbolic tangent.
$ Hwang 2011 : Namelist FLX4
$ CDFAC : re-scaling of drag
$
$ Linear input - - - - - - - - - - - - - - - - - - - - - - - - - - - -
$ Cavaleri and M-R : Namelist SLN1
$ CLIN : Proportionality constant.
$ RFPM : Factor for fPM in filter.
$ RFHF : Factor for fh in filter.
$
$ Exponential input - - - - - - - - - - - - - - - - - - - - - - - - -
$ WAM-3 : Namelist SIN1
$ CINP : Proportionality constant.
$
$ Tolman and Chalikov : Namelist SIN2
$ ZWND : Height of wind (m).
$ SWELLF : swell factor in (n.nn).
$ STABSH, STABOF, CNEG, CPOS, FNEG :
$ c0, ST0, c1, c2 and f1 in . (n.nn)
$ through (2.65) for definition of
$ effective wind speed (!/STAB2).
$ WAM4 and variants : Namelist SIN3
$ ZWND : Height of wind (m).
$ ALPHA0 : minimum value of Charnock coefficient
$ Z0MAX : maximum value of air-side roughness z0
$ BETAMAX : maximum value of wind-wave coupling
$ SINTHP : power of cosine in wind input
$ ZALP : wave age shift to account for gustiness
$ TAUWSHELTER : sheltering of short waves to reduce u_star
$ SWELLFPAR : choice of swell attenuation formulation
$ (1: TC 1996, 3: ACC 2008)
$ SWELLF : swell attenuation factor
$ Extra parameters for SWELLFPAR=3 only
$ SWELLF2, SWELLF3 : swell attenuation factors
$ SWELLF4 : Threshold Reynolds number for ACC2008
$ SWELLF5 : Relative viscous decay below threshold
$ Z0RAT : roughness for oscil. flow / mean flow
$ BYDRZ input : Namelist SIN6
$ SINA0 : factor for negative input
$ SINU10 : wind speed scaling option
$
$ Nonlinear interactions - - - - - - - - - - - - - - - - - - - - - - -
$ Discrete I.A. : Namelist SNL1
$ LAMBDA : Lambda in source term.
$ NLPROP : C in sourc term. NOTE : default
$ value depends on other source
$ terms selected.
$ KDCONV : Factor before kd in Eq. (n.nn).
$ KDMIN, SNLCS1, SNLCS2, SNLCS3 :
$ Minimum kd, and constants c1-3
$ in depth scaling function.
$ Exact interactions : Namelist SNL2
$ IQTYPE : Type of depth treatment
$ 1 : Deep water
$ 2 : Deep water / WAM scaling
$ 3 : Shallow water
$ TAILNL : Parametric tail power.
$ NDEPTH : Number of depths in for which
$ integration space is established.
$ Used for IQTYPE = 3 only
$ Namelist ANL2
$ DEPTHS : Array with depths for NDEPTH = 3
$ Gen. Multiple DIA : Namelist SNL3
$ NQDEF : Number of quadruplets.
$ MSC : Scaling constant 'm'.
$ NSC : Scaling constant 'N'.
$ KDFD : Deep water relative filter depth,
$ KDFS : Shallow water relative filter depth,
$ Namelist ANL3
$ QPARMS : 5 x NQDEF paramaters describing the
$ quadruplets, repeating LAMBDA, MU, DT12.
$ Cdeep and Cshal. See examples below.
$ Two Scale Approx. : Namelist SNL4
$ INDTSA : Index for TSA/FBI computations
$ (0 = FBI ; 1 = TSA)
$ ALTLP : Index for alternate looping
$ (1 = no ; 2 = yes)
$
$ Traditional DIA setup (default):
$
$ &SNL3 NQDEF = 1, MSC = 0.00, NSC = -3.50 /
$ &ANL3 QPARMS = 0.250, 0.000, -1.0, 0.1000E+08, 0.0000E+00 /
$
$ GMD3 from 2010 report (G13d in later paper) :
$
$ &SNL3 NQDEF = 3, MSC = 0.00, NSC = -3.50 /
$ &ANL3 QPARMS = 0.126, 0.000, -1.0, 0.4790E+08, 0.0000E+00 ,
$ 0.237, 0.000, -1.0, 0.2200E+08, 0.0000E+00 ,
$ 0.319, 0.000, -1.0, 0.1110E+08, 0.0000E+00 /
$
$ G35d from 2010 report:
$
$ &SNL3 NQDEF = 5, MSC = 0.00, NSC = -3.50 /
$ &ANL3 QPARMS = 0.066, 0.018, 21.4, 0.170E+09, 0.000E+00 ,
$ 0.127, 0.069, 19.6, 0.127E+09, 0.000E+00 ,
$ 0.228, 0.065, 2.0, 0.443E+08, 0.000E+00 ,
$ 0.295, 0.196, 40.5, 0.210E+08, 0.000E+00 ,
$ 0.369, 0.226, 11.5, 0.118E+08, 0.000E+00 /
$
$ Nonlinear filter based on DIA - - - - - - - - - - - - - - - - - - -
$ Namelist SNLS
$ A34 : Relative offset in quadruplet
$ FHFC : Proportionality constants.
$ DMN : Maximum relative change.
$ FC1-3 : Constants in frequency filter.
$
$ Whitecapping dissipation - - - - - - - - - - - - - - - - - - - - -
$ WAM-3 : Namelist SDS1
$ CDIS, APM : As in source term.
$
$ Tolman and Chalikov : Namelist SDS2
$ SDSA0, SDSA1, SDSA2, SDSB0, SDSB1, PHIMIN :
$ Constants a0, a1, a2, b0, b1 and
$ PHImin.
$
$ WAM4 and variants : Namelist SDS3
$ SDSC1 : WAM4 Cds coeffient
$ MNMEANP, WNMEANPTAIL : power of wavenumber
$ for mean definitions in Sds and tail
$ SDSDELTA1, SDSDELTA2 : relative weights
$ of k and k^2 parts of WAM4 dissipation
$ SDSLF, SDSHF : coefficient for activation of
$ WAM4 dissipation for unsaturated (SDSLF) and
$ saturated (SDSHF) parts of the spectrum
$ SDSC2 : Saturation dissipation coefficient
$ SDSC4 : Value of B0=B/Br for wich Sds is zero
$ SDSBR : Threshold Br for saturation
$ SDSP : power of (B/Br-B0) in Sds
$ SDSBR2 : Threshold Br2 for the separation of
$ WAM4 dissipation in saturated and non-saturated
$ SDSC5 : coefficient for turbulence dissipation
$ SDSC6 : Weight for the istropic part of Sds_SAT
$ SDSDTH: Angular half-width for integration of B
$
$ BYDRZ : Namelist SDS6
$ SDSET : Select threshold normalization spectra
$ SDSA1, SDSA2, SDSP1, SDSP2 :
$ Coefficients for dissipation terms T1 and T2
$ : Namelist SWL6
$ SWLB1 : Coefficient for swell dissipation
$
$ Bottom friction - - - - - - - - - - - - - - - - - - - - - - - - - -
$ JONSWAP : Namelist SBT1
$ GAMMA : Bottom friction emprical constant
$
$
$ Surf breaking - - - - - - - - - - - - - - - - - - - - - - - - - - -
$ Battjes and Janssen : Namelist SDB1
$ BJALFA : Dissipation constant (default = 1)
$ BJGAM : Breaking threshold (default = 0.73)
$ BJFLAG : TRUE - Use Hmax/d ratio only (default)
$ FALSE - Use Hmax/d in Miche formulation
$
$ Dissipation in the ice - - - - - - - - - - - - - - - - - - - - - -
$ Generalization of Liu et al. : Namelist SIC2
$ IC2DISPER : If true uses Liu formulation with eddy viscosity
$ If false, uses the generalization with turbulent
$ to laminar transition
$ IC2TURB : empirical factor for the turbulent part
$ IC2ROUGH : under-ice roughness length
$ IC2REYNOLDS: Re number for laminar to turbulent transition
$ IC2SMOOTH : smoothing of transition reprensenting random waves
$ IC2VISC : empirical factor for viscous part
$
$
$ Scattering in the ice & creep dissipations- - - - - - - - - - - - -
$ Generalization of Wiliams et al. : Namelist SIS2
$ ISC1 : scattering coefficient (default = 1)
$ IS2BACKSCAT : fraction of energy back-scattered (default = 1 )
$ IS2BREAK : TRUE - changes floe max diameter
$ : FALSE - does not change floe max diameter
$ IS2C1 : scattering in pack ice
$ IS2C2 : frequency dependance of scattering in pack ice
$ IS2C3 : frequency dependance of scattering in pack ice
$ ISBACKSCAT : fraction of scattered energy actualy redistributed
$ IS2DISP : use of ice-specific dispersion relation (T/F)
$ FRAGILITY : parameter between 0 and 1 that gives the shape of FSD
$ IS2DMIN : minimum floe diameter in meters
$ IS2DAMP : multiplicative coefficient for dissipation term from RP
$ IS2UPDATE : TRUE - updates the max floe diameter with forcing only
$ : FALSE - updates the max floe diameter at each time step
$
$ Triad nonlinear interactions - - - - - - - - - - - - - - - - - - - -
$ Lumped Triad Interaction (LTA) : Namelist STR1 (To be implemented)
$ PTRIAD1 : Proportionality coefficient (default 0.05)
$ PTRIAD2 : Multiple of Tm01 up to which interaction
$ is computed (2.5)
$ PTRIAD3 : Ursell upper limit for computing
$ interactions (not used, default 10.)
$ PTRIAD4 : Shape parameter for biphase
$ computation (0.2)
$ PTRIAD5 : Ursell number treshold for computing
$ interactions (0.01)
$
$ Shoreline reflections - - - - - - - - - - - - - - - - - - - - - - - -
$ ref. parameters : Namelist REF1
$ REFCOAST : Reflection coefficient at shoreline
$ REFFREQ : Activation of freq-dependent ref.
$ REFMAP : Scale factor for bottom slope map
$ REFRMAX : maximum ref. coeffient (default 0.8)
$ REFFREQPOW: power of frequency
$ REFICEBERG: Reflection coefficient for icebergs
$ REFSUBGRID: Reflection coefficient for islands
$ REFCOSP_STRAIGHT: power of cosine used for
$ straight shoreline
$
$ Bound 2nd order spectrum and free IG - - - - - - - - - - - - - - - - -
$ IG1 parameters : Namelist SIG1
$ IGMETHOD : 1: Hasselmann, 2: Krasitskii-Janssen
$ IGADDOUTP : activation of bound wave correction
$ in ww3_outp / ww3_ounp
$ IGSOURCE : 1: uses bound waves, 2: empirical
$ IGSTERMS : > 0 : no source term in IG band
$ IGMAXFREQ : maximum frequency of IG band
$ IGEMPIRICAL: constant in empirical free IG source
$ IGBCOVERWRITE: T: Replaces IG spectrum, does not add
$ IGSWELLMAX: T: activates free IG sources for all freq.
$
$
$ Propagation schemes ------------------------------------------------ $
$ First order : Namelist PRO1
$ CFLTM : Maximum CFL number for refraction.
$
$ UQ/UNO with diffusion : Namelist PRO2
$ CFLTM : Maximum CFL number for refraction.
$ DTIME : Swell age (s) in garden sprinkler
$ correction. If 0., all diffusion
$ switched off. If small non-zero
$ (DEFAULT !!!) only wave growth
$ diffusion.
$ LATMIN : Maximum latitude used in calc. of
$ strength of diffusion for prop.
$
$ UQ/UNO with averaging : Namelist PRO3
$ CFLTM : Maximum CFL number for refraction.
$ WDTHCG : Tuning factor propag. direction.
$ WDTHTH : Tuning factor normal direction.
$
$ Note that UQ and UNO schemes have no tunable parameters.
$ All tuneable parameters are associated with the refraction
$ limitation and the GSE alleviation.
$
$ Unstructured grids ------------------------------------------------ $
$ UNST parameters : Namelist UNST
$ UGOBCAUTO : TRUE: OBC points are taken from type 15 elements
$ FALSE: OBC points must be listed in ww3_grid.inp
$ UGOBCDEPTH: Threshold ( < 0) depth for OBC points if UGOBCAUTO is TRUE
$ EXPFSN : Activation of N scheme
$ EXPFSPSI : Activation of PSI scheme
$ EXPFSFCT : Activation of FCT scheme
$ IMPFSN : Activation of N implicit scheme
$
$ SMC grid propagation : Namelist PSMC and default values
$ CFLTM : Maximum CFL no. for propagation, 0.7
$ DTIME : Swell age for diffusion term (s), 0.0
$ LATMIN : Maximum latitude (deg) for GCT, 86.0
$ RFMAXD : Maximum refraction turning (deg), 80.0
$ LvSMC : No. of refinement level, default 1
$ ISHFT : Shift number of i-index, default 0
$ JEQT : Shift number of j-index, default 0
$ NBISMC : No. of input boundary points, 0
$ UNO3 : Use 3rd order advection scheme, .FALSE.
$ UNO3 : Add extra spatial averaging, .FALSE.
$ &PSMC DTIME = 39600.0, LATMIN=85.0, RFMAXD = 36.0, LvSMC=3, JEQT=1344 /
$
$ Output of 3D arrays------------------------------------------------- $
$ In order to limit the use of memory, arrays for 3D output fiels (i.e.
$ variables that are a function of both space and frequency, are not
$ declared, and thus cannot be used, unless specified by namelists.
$ NB: Output of 'first 5' moments E, th1m, sth1m, th2, sth2m allows to estimate the full
$ directional spectrum using, e.g. MEM (Lygre&Krogstad 1986).
$
$ Parameters (integers) : Namelist OUTS
$ For the frequency spectrum E(f)
$ E3D : <=0: not declared, > 0: declared
$ I1E3D : First frequency index of output (default is 1)
$ I2E3D : Last frequency index of output (default is NK)
$ For the mean direction th1m(f), and spread sth1m(f)
$ TH1MF, STH1MF : <=0: not declared, > 0: declared
$ I1TH1MF, I1STH1MF: First frequency index of output (default is 1)
$ I2TH1MF, I2STH1MF: First frequency index of output (default is 1)
$ For the mean direction th2m(f), and spread sth2m(f)
$ TH2MF, STH2MF : <=0: not declared, > 0: declared
$ I1TH2MF, I1STH2MF: First frequency index of output (default is 1)
$ I2TH2MF, I2STH2MF: First frequency index of output (default is 1)
$ For 2nd order pressure at K=0 (source of microseisms & microbaroms)
$ P2SF : <=0: not declared, > 0: declared
$ I1P2SF : First frequency index of output (default is 1)
$ I2P2SF : Last frequency index of output (default is NK)
$
$ Miscellaneous ------------------------------------------------------ $
$ Misc. parameters : Namelist MISC
$ CICE0 : Ice concentration cut-off.
$ CICEN : Ice concentration cut-off.
$ PMOVE : Power p in GSE aleviation for
$ moving grids in Eq. (D.4).
$ XSEED : Xseed in seeding alg. (!/SEED).
$ FLAGTR : Indicating presence and type of
$ subgrid information :
$ 0 : No subgrid information.
$ 1 : Transparancies at cell boun-
$ daries between grid points.
$ 2 : Transp. at cell centers.
$ 3 : Like 1 with cont. ice.
$ 4 : Like 2 with cont. ice.
$ XP, XR, XFILT
$ Xp, Xr and Xf for the dynamic
$ integration scheme.
$ IHMAX : Number of discrete levels in part.
$ HSPMIN : Minimum Hs in partitioning.
$ WSM : Wind speed multiplier in part.
$ WSC : Cut of wind sea fraction for
$ identifying wind sea in part.
$ FLC : Flag for combining wind seas in
$ partitioning.
$ NOSW : Number of partitioned swell fields
$ in field output.
$ FMICHE : Constant in Miche limiter.
$ STDX : Space-Time Extremes X-Length
$ STDY : Space-Time Extremes Y-Length
$ STDT : Space-Time Extremes Duration
$ P2SF : ......
$
$ Diagnostic Sea-state Dependent Stress- - - - - - - - - - - - - - - - -
$ Reichl et al. 2014 : Namelist FLD1
$ TAILTYPE : High Frequency Tail Method
$ 0: Constant value (prescribed)
$ 1: Wind speed dependent
$ (Based on GFDL Hurricane
$ Model Z0 relationship)
$ TAILLEV : Level of high frequency tail
$ (if TAILTYPE==0)
$ Valid choices:
$ Capped min: 0.001, max: 0.02
$ TAILT1 : Tail transition ratio 1
$ TAILT1*peak input frequency
$ is the first transition point of
$ the saturation specturm
$ Default is 1.25
$ TAILT1 : Tail transition ratio 2
$ TAILT2*peak input frequency
$ is the second transition point of
$ the saturation specturm
$ Default is 3.00
$ Donelan et al. 2012 : Namelist FLD2
$ TAILTYPE : See above (FLD1)
$ TAILLEV : See above (FLD1)
$ TAILT1 : See above (FLD1)
$ TAILT2 : See above (FLD1)
$
$ In the 'Out of the box' test setup we run with sub-grid obstacles
$ and with continuous ice treatment.
$
&MISC CICE0 = 0.25, CICEN = 0.75, FLAGTR = 4 /
&FLX3 CDMAX = 3.5E-3 , CTYPE = 0 /
$ &SDB1 BJGAM = 1.26, BJFLAG = .FALSE. /
$
$ Mandatory string to identify end of namelist input section.
$
END OF NAMELISTS
$