-
Notifications
You must be signed in to change notification settings - Fork 505
/
ww3_shel.inp
378 lines (378 loc) · 16.9 KB
/
ww3_shel.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
$ -------------------------------------------------------------------- $
$ WAVEWATCH III shel input file $
$ -------------------------------------------------------------------- $
$ Define input to be used with F/T/C flag for use or nor or coupling and
$ T/F flag for definition as a homogeneous field.
$
$ Include ice and mud parameters only if IC1/2/3/4 used :
F F Ice parameter 1
F F Ice parameter 2
F F Ice parameter 3
F F Ice parameter 4
F F Ice parameter 5
F F Mud parameter 1
F F Mud parameter 2
F F Mud parameter 3
F F Water levels
F F Currents
T T Winds
T F Ice concentrations
F F Atmospheric momentum
F F Air density
F Assimilation data : Mean parameters
F Assimilation data : 1-D spectra
F Assimilation data : 2-D spectra
$
$ Time frame of calculations ----------------------------------------- $
$ - Starting time in yyyymmdd hhmmss format.
$ - Ending time in yyyymmdd hhmmss format.
$
19680606 000000
19680606 060000
$
$ Define output data ------------------------------------------------- $
$
$ Define output server mode. This is used only in the parallel version
$ of the model. To keep the input file consistent, it is always needed.
$ IOSTYP = 1 is generally recommended. IOSTYP > 2 may be more efficient
$ for massively parallel computations. Only IOSTYP = 0 requires a true
$ parallel file system like GPFS.
$
$ IOSTYP = 0 : No data server processes, direct access output from
$ each process (requires true parallel file system).
$ 1 : No data server process. All output for each type
$ performed by process that performs computations too.
$ 2 : Last process is reserved for all output, and does no
$ computing.
$ 3 : Multiple dedicated output processes.
$
2
$
$ Five output types are available (see below). All output types share
$ a similar format for the first input line:
$ - first time in yyyymmdd hhmmss format, output interval (s), and
$ last time in yyyymmdd hhmmss format (all integers).
$ Output is disabled by setting the output interval to 0.
$
$ ------------------------------------------------------------------- $
$
$ Type 1 : Fields of mean wave parameters
$ Standard line and line with logical flags to activate output
$ fields as defined in section 2.4 of the manual. The logical
$ flags are not supplied if no output is requested. The logical
$ flags can be placed on multiple consecutive lines. However,
$ the total number and order of the logical flags is fixed.
$ The raw data file is out_grd.ww3,
$ see w3iogo.ftn for additional doc.
$
19680606 000000 3600 19680608 000000
$----------------------------------------------------------------
$ Output request flags identifying fields.
$
$ The table below provides a full definition of field output parameters
$ as well as flags indicating if they are available in different field
$ output output file types (ASCII, grib, NetCDF).
$ Further definitions are found in section 2.4 of the manual.
$
$ Selection of field outputs may be made in two ways:
$ F/T flags: first flag is set to F, requests made per group (1st line)
$ followed by parameter flags (total of 10 groups).
$ Namelists: first line is set to N, next line contains parameter
$ symbol as per table below.
$
$ Example of F/T flag use is given in this sample ww3_shel.inp, below.
$ For namelist usage, see the sample ww3_ounf.inp for an example.
$
$ ----------------------------------------
$ Output field parameter definitions table
$ ----------------------------------------
$
$ All parameters listed below are available in output file of the types
$ ASCII and NetCDF. If selected output file types are grads or grib,
$ some parameters may not be available. The first two columns in the
$ table below identify such cases by flags, cols 1 (GRB) and 2 (GXO)
$ refer to grib (ww3_grib) and grads (gx_outf), respectively.
$
$ Columns 3 and 4 provide group and parameter numbers per group.
$ Columns 5, 6 and 7 provide:
$ 5 - code name (internal)
$ 6 - output tags (names used is ASCII file extensions, NetCDF
$ variable names and namelist-based selection (see ww3_ounf.inp)
$ 7 - Long parameter name/definition
$
$ G G
$ R X Grp Param Code Output Parameter/Group
$ B O Numb Numbr Name Tag Definition
$ --------------------------------------------------
$ 1 Forcing Fields
$ -------------------------------------------------
$ T T 1 1 DW DPT Water depth.
$ T T 1 2 C[X,Y] CUR Current velocity.
$ T T 1 3 UA WND Wind speed.
$ T T 1 4 AS AST Air-sea temperature difference.
$ T T 1 5 WLV WLV Water levels.
$ T T 1 6 ICE ICE Ice concentration.
$ T T 1 7 IBG IBG Iceberg-induced damping.
$ T T 1 8 TAUA TAU Atm. momentum.
$ T T 1 9 RHOAIR RHO Air density.
$ T T 1 10 D50 D50 Median sediment grain size.
$ T T 1 11 IC1 IC1 Ice thickness.
$ T T 1 12 IC5 IC5 Ice flow diameter.
$ -------------------------------------------------
$ 2 Standard mean wave Parameters
$ -------------------------------------------------
$ T T 2 1 HS HS Wave height.
$ T T 2 2 WLM LM Mean wave length.
$ T T 2 3 T02 T02 Mean wave period (Tm0,2).
$ T T 2 4 T0M1 T0M1 Mean wave period (Tm0,-1).
$ T T 2 5 T01 T01 Mean wave period (Tm0,1).
$ T T 2 6 FP0 FP Peak frequency.
$ T T 2 7 THM DIR Mean wave direction.
$ T T 2 8 THS SPR Mean directional spread.
$ T T 2 9 THP0 DP Peak direction.
$ T T 2 10 HIG HIG Infragravity height
$ T T 2 11 STMAXE MXE Max surface elev (STE)
$ T T 2 12 STMAXD MXES St Dev of max surface elev (STE)
$ T T 2 13 HMAXE MXH Max wave height (STE)
$ T T 2 14 HCMAXE MXHC Max wave height from crest (STE)
$ T T 2 15 HMAXD SDMH St Dev of MXC (STE)
$ T T 2 16 HCMAXD SDMHC St Dev of MXHC (STE)
$ F T 2 17 WBT WBT Dominant wave breaking probability bT
$ F F 2 18 FP0 TP Peak period (from peak freq)
$ F F 2 19 WNMEAN WNM Mean wavenumber
$ -------------------------------------------------
$ 3 Spectral Parameters (first 5)
$ -------------------------------------------------
$ F F 3 1 EF EF Wave frequency spectrum
$ F F 3 2 TH1M TH1M Mean wave direction from a1,b2
$ F F 3 3 STH1M STH1M Directional spreading from a1,b2
$ F F 3 4 TH2M TH2M Mean wave direction from a2,b2
$ F F 3 5 STH2M STH2M Directional spreading from a2,b2
$ F F 3 6 WN WN Wavenumber array
$ -------------------------------------------------
$ 4 Spectral Partition Parameters
$ -------------------------------------------------
$ T T 4 1 PHS PHS Partitioned wave heights.
$ T T 4 2 PTP PTP Partitioned peak period.
$ T T 4 3 PLP PLP Partitioned peak wave length.
$ T T 4 4 PDIR PDIR Partitioned mean direction.
$ T T 4 5 PSI PSPR Partitioned mean directional spread.
$ T T 4 6 PWS PWS Partitioned wind sea fraction.
$ T T 4 7 PTHP0 PDP Peak wave direction of partition.
$ T T 4 8 PQP PQP Goda peakdedness parameter of partition.
$ T T 4 9 PPE PPE JONSWAP peak enhancement factor of partition.
$ T T 4 10 PGW PGW Gaussian frequency width of partition.
$ T T 4 11 PSW PSW Spectral width of partition.
$ T T 4 12 PTM1 PTM10 Mean wave period (m-1,0) of partition.
$ T T 4 13 PT1 PT01 Mean wave period (m0,1) of partition.
$ T T 4 14 PT2 PT02 Mean wave period (m0,2) of partition.
$ T T 4 15 PEP PEP Peak spectral density of partition.
$ T T 4 16 PWST TWS Total wind sea fraction.
$ T T 4 17 PNR PNR Number of partitions.
$ -------------------------------------------------
$ 5 Atmosphere-waves layer
$ -------------------------------------------------
$ T T 5 1 UST UST Friction velocity.
$ F T 5 2 CHARN CHA Charnock parameter
$ F T 5 3 CGE CGE Energy flux
$ F T 5 4 PHIAW FAW Air-sea energy flux
$ F T 5 5 TAUWI[X,Y] TAW Net wave-supported stress
$ F T 5 6 TAUWN[X,Y] TWA Negative part of the wave-supported stress
$ F F 5 7 WHITECAP WCC Whitecap coverage
$ F F 5 8 WHITECAP WCF Whitecap thickness
$ F F 5 9 WHITECAP WCH Mean breaking height
$ F F 5 10 WHITECAP WCM Whitecap moment
$ F F 5 11 FWS FWS Wind sea mean period
$ -------------------------------------------------
$ 6 Wave-ocean layer
$ -------------------------------------------------
$ F F 6 1 S[XX,YY,XY] SXY Radiation stresses.
$ F F 6 2 TAUO[X,Y] TWO Wave to ocean momentum flux
$ F F 6 3 BHD BHD Bernoulli head (J term)
$ F F 6 4 PHIOC FOC Wave to ocean energy flux
$ F F 6 5 TUS[X,Y] TUS Stokes transport
$ F F 6 6 USS[X,Y] USS Surface Stokes drift
$ F F 6 7 [PR,TP]MS P2S Second-order sum pressure
$ F F 6 8 US3D USF Spectrum of surface Stokes drift
$ F F 6 9 P2SMS P2L Micro seism source term
$ F F 6 10 TAUICE TWI Wave to sea ice stress
$ F F 6 11 PHICE FIC Wave to sea ice energy flux
$ F F 6 12 USSP USP Partitioned surface Stokes drift
$ F F 6 13 TAUOC[X,Y] TOC Total momentum to the ocean
$ -------------------------------------------------
$ 7 Wave-bottom layer
$ -------------------------------------------------
$ F F 7 1 ABA ABR Near bottom rms amplitides.
$ F F 7 2 UBA UBR Near bottom rms velocities.
$ F F 7 3 BEDFORMS BED Bedforms
$ F F 7 4 PHIBBL FBB Energy flux due to bottom friction
$ F F 7 5 TAUBBL TBB Momentum flux due to bottom friction
$ -------------------------------------------------
$ 8 Spectrum parameters
$ -------------------------------------------------
$ F F 8 1 MSS[X,Y] MSS Mean square slopes
$ F F 8 2 MSC[X,Y] MSC Spectral level at high frequency tail
$ F F 8 3 MSSD MSD Slope direction
$ F F 8 4 MSCD MCD Tail slope direction
$ F F 8 5 QP QP Goda peakedness parameter
$ F F 8 6 QKK QKK Wavenumber peakedness
$ F F 8 7 SKEW SKW Skewness of elevation for zero slopes
$ F F 8 8 EMBIA1 EMB Mean sea level at zero slopes / Hs
$ F F 8 9 EMBIA2 EMC Tracker bias for LRM least square altimetry
$ -------------------------------------------------
$ 9 Numerical diagnostics
$ -------------------------------------------------
$ T T 9 1 DTDYN DTD Average time step in integration.
$ T T 9 2 FCUT FC Cut-off frequency.
$ T T 9 3 CFLXYMAX CFX Max. CFL number for spatial advection.
$ T T 9 4 CFLTHMAX CFD Max. CFL number for theta-advection.
$ F F 9 5 CFLKMAX CFK Max. CFL number for k-advection.
$ -------------------------------------------------
$ 10 User defined
$ -------------------------------------------------
$ F F 10 1 U1 User defined #1. (requires coding ...)
$ F F 10 2 U2 User defined #1. (requires coding ...)
$ -------------------------------------------------
$
$ Section 4 consist of a set of fields, index 0 = wind sea, index
$ 1:NOSWLL are first NOSWLL swell fields.
$
$ Actual active parameter selection section
$
$ (1) Forcing Fields
T
$ DPT CUR WND AST WLV ICE IBG TAU RHO D50 IC1 IC5
T T T T T F F F F F F F
$ (2) Standard mean wave Parameters
T
$ HS LM T02 T0M1 T01 FP DIR SPR DP
T T T T T T T T T
$ (3) Frequency-dependent parameters
T
$ EF TH1M STH1M TH2M STH2M WN
T T T F F F
$ (4) Spectral Partition Parameters
T
$ PHS PTP PLP PDIR PSPR PNR PDP PQP PPE PGW PSW PTM10 PT01 PT02 PEP PWS TWS
T T T T T T T T T T T T T T T T T
$ (5) Atmosphere-waves layer
T
$ UST CHA CGE FAW TAW TWA WCC WCF WCH WCM FWS
T T T T T T T T T T T
$ (6) Wave-Ocean layer
T
$ SXY TWO BHD FOC TUS USS P2S USF P2L TWI FIC USP
T T T T T T T F F F F T
$ (7) Wave-bottom layer
T
$ ABR UBR BED FBB TBB
T T T T T
$ (8) Spectrum parameters
T
$ MSS MSC WL02 AXT AYT AXY
T T T T T T
$ (9) Numerical diagnostics
T
$ DTD FC CFX CFD CFK
T T T T T
$ (10) User defined (NOEXTR flags needed)
F
$ U1 U2
$ T T
$
$----------------------------------------------------------------
$
$ Type 2 : Point output
$ Standard line and a number of lines identifying the
$ longitude, latitude and name (C*40) of output points.
$ The list is closed by defining a point with the name
$ 'STOPSTRING'. No point info read if no point output is
$ requested (i.e., no 'STOPSTRING' needed).
$ Example for spherical grid.
$ The raw data file is out_pnt.ww3,
$ see w3iogo.ftn for additional doc.
$
$ NOTE : Spaces may be included in the name, but this is not
$ advised, because it will break the GrADS utility to
$ plots spectra and source terms, and will make it more
$ difficult to use point names in data files.
$
19680606 000000 900 19680608 000000
$
-0.25 -0.25 'Land '
0.0 0.0 'Point_1 '
2.0 1.0 'Point_2 '
1.8 2.2 'Point_3 '
2.1 0.9 'Point_4 '
5.0 5.0 'Outside '
$
0.0 0.0 'STOPSTRING'
$
$ Type 3 : Output along track.
$ Flag for formatted input file.
$ The data files are track_i.ww3 and
$ track_o.ww3, see w3iotr.ftn for ad. doc.
$
19680606 000000 1800 19680606 013000
T
$
$ Type 4 : Restart files (no additional data required).
$ The data file is restartN.ww3, see
$ w3iors.ftn for additional doc.
$ Flag for second restart (data server mode)
$ Flag for extra fields in the restart file (coupling restart)
$
19680606 030000 3600 19680607 030000 F F
$
$ Type 5 : Boundary data (no additional data required).
$ The data file is nestN.ww3, see
$ w3iobcmd.ftn for additional doc.
$
19680606 000000 3600 20010102 000000
$
$ Type 6 : Separated wave field data (dummy for now).
$ First, last step IX and IY, flag for formatted file
$
19680606 000000 3600 20010102 000000
0 999 1 0 999 1 T
$
$ Type 7 : Coupling. (must be fully commented if not used with switch COU)
$ Couple at T+0 flag
$ Namelist type selection is used here.
$ Diagnostic fields to exchange. (see namcouple for more information)
$
$ 19680606 000000 3600 20010102 000000 F
$ N
$
$ - Sent fields by ww3:
$ - Ocean model : T0M1 OCHA OHS DIR BHD TWO UBR FOC TAW TUS USS LM DRY
$ - Atmospheric model : ACHA AHS TP (or FP) FWS
$ - Ice model : IC5 TWI
$
$ CHA
$
$ - Received fields by ww3:
$ - Ocean model : SSH CUR
$ - Atmospheric model : WND TAU RHO
$ - Ice model : ICE IC1 IC5
$
$ WND
$
$ Homogeneous field data --------------------------------------------- $
$ Homogeneous fields can be defined by a list of lines containing an ID
$ string 'LEV' 'CUR' 'WND', date and time information (yyyymmdd
$ hhmmss), value (S.I. units), direction (current and wind, meteorol.
$ convention degrees)) and air-sea temperature difference (degrees C).
$ 'STP' is mandatory stop string.
$ Also defined here are the speed with which the grid is moved
$ continuously, ID string 'MOV', parameters as for 'CUR'.
$
'LEV' 19680606 010000 1.00
'CUR' 19680606 073125 2.0 25.
'WND' 19680606 000000 20. 145. 2.0
'MOV' 19680606 013000 4.0 25.
'STP'
$
$ -------------------------------------------------------------------- $
$ End of input file $
$ -------------------------------------------------------------------- $