-
Notifications
You must be signed in to change notification settings - Fork 555
/
Copy pathww3_grid.inp
880 lines (879 loc) · 44.6 KB
/
ww3_grid.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
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
$ -------------------------------------------------------------------- $
$ WAVEWATCH III Grid preprocessor input file $
$ -------------------------------------------------------------------- $
$ Grid name (C*30, in quotes)
$
'TEST GRID (GULF OF NOWHERE) '
$
$ Frequency increment factor and first frequency (Hz) ---------------- $
$ number of frequencies (wavenumbers) and directions, relative offset
$ of first direction in terms of the directional increment [-0.5,0.5].
$ In versions 1.18 and 2.22 of the model this value was by definiton 0,
$ it is added to mitigate the GSE for a first order scheme. Note that
$ this factor is IGNORED in the print plots in ww3_outp.
$
1.1 0.04118 25 24 0.
$
$ Set model flags ---------------------------------------------------- $
$ - FLDRY Dry run (input/output only, no calculation).
$ - FLCX, FLCY Activate X and Y component of propagation.
$ - FLCTH, FLCK Activate direction and wavenumber shifts.
$ - FLSOU Activate source terms.
$
F T T T F T
$
$ Set time steps ----------------------------------------------------- $
$ - Time step information (this information is always read)
$ maximum global time step, maximum CFL time step for x-y and
$ k-theta, minimum source term time step (all in seconds).
$
900. 950. 900. 300.
$
$ 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
$ CAPCHA : enable functional form for Cd (1=enabled),
$ only valid when SINTABLE=0.
$ UCAP : Wind speed threshold for cap (CAPCHA=1)
$ SIGMAUCAP : Width of functional form (CAPCHA=1)
$ CHA0 : Initial Charnock coefficient (CAPCHA=1)
$ CHAMIN : Minimum Charnock value (CAPCHA=1)
$ 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
$ SINWS : wind speed scaling option
$ SINFC : high-frequency extent of the
$ prognostic frequency region
$
$ 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.
$ IQTYPE : Type of depth treatment
$ -2 : Deep water GQM with scaling
$ 1 : Deep water DIA
$ 2 : Deep water DIA with scaling
$ 3 : Shallow water DIA
$ TAILNL : Parametric tail power.
$ GQMNF1 : number of points along the locus
$ GQMNT1 : idem
$ GQMNQ_OM2 : idem
$ GQMTHRSAT : threshold on saturation for SNL calculation
$ GQMTHRCOU : threshold for filter on coupling coefficient
$ GQAMP1, GQAMP2, GQAMP3, GQAMP4 : amplification factor
$ 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
$
$ Dissipation by sea ice
$ Empirical/parametric representations : Namelist SIC4
$ IC4METHOD : integer 1 to 7
$ : In most cases, additional input
$ : is required.
$ : See examples in /regtests/ww3_tic1.1/
$ : See also: 1) description in manual
$ : and 2) inline documentation in
$ w3sic4md.ftn
$
$ 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
$ UGBCCFL : Turns on/off (TRUE/FALSE) the computation of the CFL number on the physical domain boundary.
$ If FALSE the explicit scheme can be much faster though stability is not guaranteed
$ (default TRUE)
$ UGOBCAUTO : TRUE: OBC points are taken from type 15 elements (default)
$ FALSE: OBC points must be listed in ww3_grid.inp
$ UGOBCDEPTH : Threshold ( < 0) depth (default -10) for OBC points if UGOBCAUTO is TRUE
$ UGOBCFILE : File name of file for reading boudary (default 'unset')
$
$ The following are TRUE/FALSE variables and only one can be TRUE.
$ By default, EXPFSN is TRUE and must be set to false to activate another option
$ EXPFSN : Activation of N scheme (default option)
$ EXPFSPSI : Activation of PSI scheme
$ EXPFSFCT : Activation of FCT scheme
$ IMPFSN : Activation of N implicit scheme
$ EXPTOTAL : Activation of Block explicit N scheme solver
$ IMPTOTAL : Activation of fully implicit scheme, non splitting
$
$ The following TRUE/FALSE variables are only for IMPTOTAL=TRUE
$ IMPREFRACTION : Turn on implicit freq. shift (default FALSE)
$ IMPFREQSHIFT : Turn on implicit freq. shift terms (default FALSE)
$ IMPSOURCE : Turn on implicit source terms (default FALSE)
$
$ JGS_TERMINATE_MAXITER : Terminate based on max number of iterations (TRUE/FALSE, default TRUE)
$ JGS_TERMINATE_DIFFERENCE : Terminate based on the total change of the unweightet sum of wave action (TRUE/FALSE, default TRUE)
$ JGS_TERMINATE_NORM : Terminate based on the norm of the solution (TRUE/FALSE, default FALSE)
$ JGS_USE_JACOBI : Use Jacobi solver family (TRUE/FALSE, default TRUE)
$ JGS_BLOCK_GAUSS_SEIDEL : Use Block Gauss Seidel method for imptotal instead of the conservative jacobi iterator. (TRUE/FALSE, default TRUE)
$ JGS_MAXITER : Max. Number of solver iterations for JGS_TERMINATE_MAXITER (integer, default 100)
$ JGS_PMIN : % of grid points that do not need to converge during solver iteration (real, default 1)
$ JGS_DIFF_THR : Implicit solver threshold for JGS_TERMINATE_DIFFERENCE (real, default 1.0e-10)
$ JGS_NORM_THR : Norm of the solution for JGS_TERMINATE_NORM (real, default 1.0e-20)
$ JGS_LIMITER : TRUE: Use total (quasi-steady: limits whole equation) instead of local limiter (un-steady: limits only source terms)
$ FALSE: default
$ JGS_LIMITER_FUNC : 1 - old limiter (default)
$ 2 - alternatnive limiter
$ SETUP_APPLY_WLV : Compute wave setup (TRUE/FALSE, default TRUE)
$ SOLVERTHR_SETUP : Solver threshold for setup computations (default 1E-6)
$ CRIT_DEP_SETUP : Critical depth for setup computations (default 0.1)
$
$ SMC grid propagation : Namelist PSMC and default values
$ CFLSM : Maximum CFL no. for propagation, 0.7
$ DTIMS : Swell age for diffusion term (s), 360.0
$ RFMAXD : Maximum refraction turning (deg), 36.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
$ AVERG : Add extra spatial averaging, .TRUE.
$ UNO3 : Use 3rd order advection scheme, .FALSE.
$ SEAWND : Use sea-point only wind input. .FALSE.
$ Arctic : Include Arctic part cell array. .FALSE.
$ &PSMC DTIMS = 39600.0, LvSMC=3, JEQT=1344, SEAWND=.TRUE., Arctic=.TRUE. /
$
$ Rotated pole ------------------------------------------------------ $
$ Pole parameters : Namelist ROTD
$ PLAT : Rotated pole latitude
$ PLON : Rotated pole longitude
$ UNROT : Logical, un-rotate directions to
$ true north
$
$ Compile switch /RTD required.
$
$ These will be used to derive rotation angle corrections in the
$ model. The corrections are used for rotation of boundary spectra
$ and for restoring conventional lat/lon orientation of the
$ output spectra, mean direction or any related variables.
$ The PLAT/LON example below is a standard setting for Met
$ Office UK regional models.
$
$ &ROTD PLAT = 37.5, PLON = 177.5, UNROT = .TRUE. /
$
$ The default values for ROTD represent a non-rotated lat/lon grid:
$
$ &ROTD PLAT = 90.0, PLON = -180.0, UNROT = .FALSE. /
$
$ Output boundary conditions to rotated pole grids ------------------ $
$
$ Pole parameters : Namelist ROTB
$ BPLAT(1:9) : Pole latitude of each destination grid
$ BPLON(1:9) : Pole longitude of each destination grid
$
$ Compile switch /RTD required.
$
$ In the section 'Output boundary points' later in this file, the points must
$ be given as straight lines in the reference system of each destination grid.
$ The index J of BPLAT(J),BPLON(J) is incremented by one where a new output
$ file is started by specifying a negative value for the number of points.
$
$ Each destination may be either a rotated grid or a standard (non-rotated)
$ lat/lon grid. The default values for ROTB represent standard lat/lon grids:
$
$ &ROTB BPLAT(1)=90., BPLON(1)=-180., BPLAT(2)=90. / $ (etc.)
$
$ 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)
$ For the surface Stokes drift partitions (USP)
$ USSP : First index (default is 1, should always be 1)
$ IUSSP : Last index (must be <= than NK and should be
$ between 3 and ~10 with the tradeoff
$ between accuracy and resources)
$ STK_WN : List of wavenumbers (size of IUSSP)
$ e.g.: USSP = 1, IUSSP=3, STK_WN = 0.04, 0.11, 0.33
$ provides 3 partitions of both x & y component,
$ with a reasonable accuracy for using in
$ a climate model.
$
$ 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.
$ TRCKCMPR : Logical variable (T/F). Set to F to
$ disable "compression" of track output.
$ This simplifies post-processing.
$ Default is T and will create track
$ output in the traditional manner
$ (WW3 v3, v4, v5).
$ 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.
$ PTM : Partioning method:
$ 1 : Default WW3
$ 2 : Watershedding + wind cutoff
$ 3 : Watershedding only
$ 4 : Wind speed cutoff only
$ 5 : High/Low band cutoff (see PTFC)
$ PTFC : Cutouf frequency for High/Low band
$ partioning (PTM=5). Default = 0.1Hz
$ FMICHE : Constant in Miche limiter.
$ STDX : Space-Time Extremes X-Length
$ STDY : Space-Time Extremes Y-Length
$ STDT : Space-Time Extremes Duration
$ P2SF : ......
$ CALTYPE: Calendar type. The only accepted
$ values are 'standard' (default),
$ '365_day', or '360_day'.
$
$ 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
$
$ Define grid -------------------------------------------------------- $
$
$ Five records containing :
$
$ 1 Type of grid, coordinate system and type of closure: GSTRG, FLAGLL,
$ CSTRG. Grid closure can only be applied in spherical coordinates.
$ GSTRG : String indicating type of grid :
$ 'RECT' : rectilinear
$ 'CURV' : curvilinear
$ 'UNST' : unstructured (triangle-based)
$ 'SMCG' : Spherical Multiple-Cell grid.
$ FLAGLL : Flag to indicate coordinate system :
$ T : Spherical (lon/lat in degrees)
$ F : Cartesian (meters)
$ CSTRG : String indicating the type of grid index space closure :
$ 'NONE' : No closure is applied
$ 'SMPL' : Simple grid closure : Grid is periodic in the
$ : i-index and wraps at i=NX+1. In other words,
$ : (NX+1,J) => (1,J). A grid with simple closure
$ : may be rectilinear or curvilinear.
$ 'TRPL' : Tripole grid closure : Grid is periodic in the
$ : i-index and wraps at i=NX+1 and has closure at
$ : j=NY+1. In other words, (NX+1,J<=NY) => (1,J)
$ : and (I,NY+1) => (NX-I+1,NY). Tripole
$ : grid closure requires that NX be even. A grid
$ : with tripole closure must be curvilinear.
$ 2 NX, NY. As the outer grid lines are always defined as land
$ points, the minimum size is 3x3.
$
$ Branch here based on grid type
$
$ IF ( RECTILINEAR GRID ) THEN
$
$ 3 Grid increments SX, SY (degr.or m) and scaling (division) factor.
$ If CSTRG='SMPL', then SX is set to 360/NX.
$ 4 Coordinates of (1,1) (degr.) and scaling (division) factor.
$
$ ELSE IF ( CURVILINEAR GRID ) THEN
$
$ 3 Unit number of file with x-coordinate.
$ Scale factor and add offset: x <= scale_fac * x_read + add_offset.
$ IDLA, IDFM, format for formatted read, FROM and filename.
$ IDLA : Layout indicator :
$ 1 : Read line-by-line bottom to top.
$ 2 : Like 1, single read statement.
$ 3 : Read line-by-line top to bottom.
$ 4 : Like 3, single read statement.
$ IDFM : format indicator :
$ 1 : Free format.
$ 2 : Fixed format with above format descriptor.
$ 3 : Unformatted.
$ FROM : file type parameter
$ 'UNIT' : open file by unit number only.
$ 'NAME' : open file by name and assign to unit.
$
$ If the above unit number equals 10, then the x-coord is read from this
$ file. The x-coord must follow the above record. No comment lines are
$ allowed within the x-coord input.
$
$ 4 Unit number of file with y-coordinate.
$ Scale factor and add offset: y <= scale_fac * y_read + add_offset.
$ IDLA, IDFM, format for formatted read, FROM and filename.
$ IDLA : Layout indicator :
$ 1 : Read line-by-line bottom to top.
$ 2 : Like 1, single read statement.
$ 3 : Read line-by-line top to bottom.
$ 4 : Like 3, single read statement.
$ IDFM : format indicator :
$ 1 : Free format.
$ 2 : Fixed format with above format descriptor.
$ 3 : Unformatted.
$ FROM : file type parameter
$ 'UNIT' : open file by unit number only.
$ 'NAME' : open file by name and assign to unit.
$
$ If the above unit number equals 10, then the y-coord is read from this
$ file. The y-coord must follow the above record. No comment lines are
$ allowed within the y-coord input.
$
$ ELSE IF ( UNSTRUCTURED GRID ) THEN
$ Nothing to declare: all the data will be read from the GMESH file
$ END IF ( CURVILINEAR GRID )
$
$ 5 Limiting bottom depth (m) to discriminate between land and sea
$ points, minimum water depth (m) as allowed in model, unit number
$ of file with bottom depths, scale factor for bottom depths (mult.),
$ IDLA, IDFM, format for formatted read, FROM and filename.
$ IDLA : Layout indicator :
$ 1 : Read line-by-line bottom to top.
$ 2 : Like 1, single read statement.
$ 3 : Read line-by-line top to bottom.
$ 4 : Like 3, single read statement.
$ IDFM : format indicator :
$ 1 : Free format.
$ 2 : Fixed format with above format descriptor.
$ 3 : Unformatted.
$ FROM : file type parameter
$ 'UNIT' : open file by unit number only.
$ 'NAME' : open file by name and assign to unit.
$
$ If the above unit number equals 10, then the bottom depths are read from
$ this file. The depths must follow the above record. No comment lines are
$ allowed within the depth input. In the case of unstructured grids, the file
$ is expected to be a GMESH grid file containing node and element lists.
$
$ ------------------------------------------------------------------------
$ Example for rectilinear grid with spherical (lon/lat) coordinate system.
$ Note that for Cartesian coordinates the unit is meters (NOT km).
$
'RECT' T 'NONE'
12 12
1. 1. 4.
-1. -1. 4.
-0.1 2.50 10 -10. 3 1 '(....)' 'NAME' 'bottom.inp'
$
6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 5 4 2 0 2 4 5 6 6
6 6 6 5 4 2 0 2 4 5 6 6
6 6 6 5 4 2 0 2 4 5 6 6
6 6 6 5 4 2 0 0 4 5 6 6
6 6 6 5 4 4 2 2 4 5 6 6
6 6 6 6 5 5 4 4 5 6 6 6
6 6 6 6 6 6 5 5 6 6 6 6
6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 6 6 6 6 6 6 6 6 6
$
$ ------------------------------------------------------------------------
$ Example for curvilinear grid with spherical (lon/lat) coordinate system.
$ Same spatial grid as preceding rectilinear example.
$ Note that for Cartesian coordinates the unit is meters (NOT km).
$
$ 'CURV' T 'NONE'
$ 12 12
$
$ 10 0.25 -0.5 3 1 '(....)' 'NAME' 'x.inp'
$
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$ 1 2 3 4 5 6 7 8 9 10 11 12
$
$ 10 0.25 0.5 3 1 '(....)' 'NAME' 'y.inp'
$
$ 1 1 1 1 1 1 1 1 1 1 1 1
$ 2 2 2 2 2 2 2 2 2 2 2 2
$ 3 3 3 3 3 3 3 3 3 3 3 3
$ 4 4 4 4 4 4 4 4 4 4 4 4
$ 5 5 5 5 5 5 5 5 5 5 5 5
$ 6 6 6 6 6 6 6 6 6 6 6 6
$ 7 7 7 7 7 7 7 7 7 7 7 7
$ 8 8 8 8 8 8 8 8 8 8 8 8
$ 9 9 9 9 9 9 9 9 9 9 9 9
$ 10 10 10 10 10 10 10 10 10 10 10 10
$ 11 11 11 11 11 11 11 11 11 11 11 11
$ 12 12 12 12 12 12 12 12 12 12 12 12
$
$ -0.1 2.50 10 -10. 3 1 '(....)' 'NAME' 'bottom.inp'
$
$ 6 6 6 6 6 6 6 6 6 6 6 6
$ 6 6 6 5 4 2 0 2 4 5 6 6
$ 6 6 6 5 4 2 0 2 4 5 6 6
$ 6 6 6 5 4 2 0 2 4 5 6 6
$ 6 6 6 5 4 2 0 0 4 5 6 6
$ 6 6 6 5 4 4 2 2 4 5 6 6
$ 6 6 6 6 5 5 4 4 5 6 6 6
$ 6 6 6 6 6 6 5 5 6 6 6 6
$ 6 6 6 6 6 6 6 6 6 6 6 6
$ 6 6 6 6 6 6 6 6 6 6 6 6
$ 6 6 6 6 6 6 6 6 6 6 6 6
$ 6 6 6 6 6 6 6 6 6 6 6 6
$
$ -------------------------------------------------------------
$ SMC grid use the same spherical lat-lon grid parameters
$ 'SMCG' T 'SMPL'
$ 1024 704
$ SMC grid base level resolution dlon dlat and start lon lat
$ 0.35156250 0.23437500 1.
$ 0.17578125 -78.6328125 1.
$
$ Normal depth input line is used to passing the minimum depth
$ though the depth file is not read for SMC grid.
$ -0.1 10.0 30 -1. 1 1 '(....)' 'NAME' 'SMC25Depth.dat'
$ SMC cell and face arrays and obstruction ratio:
$ 32 1 1 '(....)' 'S6125MCels.dat'
$ 33 1 1 '(....)' 'S6125ISide.dat'
$ 34 1 1 '(....)' 'S6125JSide.dat'
$ 31 1.0 1 1 '(...)' 'NAME' 'SMC25Subtr.dat'
$ The input boundary cell file is only needed when NBISMC > 0.
$ 35 1 1 '(....)' 'S6125Bundy.dat'
$ Extra cell and face arrays for Arctic part if Arctic=.TRUE. is selected.
$ 36 1 1 '(....)' 'S6125MBArc.dat'
$ 37 1 1 '(....)' 'S6125AISid.dat'
$ 38 1 1 '(....)' 'S6125AJSid.dat'
$ Normal land-sea mask file input line is kept but file is not used.
$ 39 1 1 '(....)' 'NAME' 'S6125Masks.dat'
$ Boundary cell id list file (unit 35) is only required if boundary
$ cell number entered above is non-zero. The cell id number should be
$ the sequential number in the cell array (unit 32) S625MCels.dat.
$
$ If sub-grid information is available as indicated by FLAGTR above,
$ additional input to define this is needed below. In such cases a
$ field of fractional obstructions at or between grid points needs to
$ be supplied. First the location and format of the data is defined
$ by (as above) :
$ - Unit number of file (can be 10, and/or identical to bottom depth
$ unit), scale factor for fractional obstruction, IDLA, IDFM,
$ format for formatted read, FROM and filename
$
10 0.2 3 1 '(....)' 'NAME' 'obstr.inp'
$
$ *** NOTE if this unit number is the same as the previous bottom
$ depth unit number, it is assumed that this is the same file
$ without further checks. ***
$
$ If the above unit number equals 10, the bottom data is read from
$ this file and follows below (no intermediate comment lines allowed,
$ except between the two fields).
$
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 5 0 0 0 0 0
0 0 0 0 0 0 5 0 0 0 0 0
0 0 0 0 0 0 4 0 0 0 0 0
0 0 0 0 0 0 4 0 0 0 0 0
0 0 0 0 0 0 5 0 0 0 0 0
0 0 0 0 0 0 5 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
$
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 5 5 5 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
$
$ *** NOTE size of fields is always NX * NY ***
$
$ Input boundary points and excluded points -------------------------- $
$ The first line identifies where to get the map data, by unit number
$ IDLA and IDFM, format for formatted read, FROM and filename
$ if FROM = 'PART', then segmented data is read from below, else
$ the data is read from file as with the other inputs (as INTEGER)
$
10 3 1 '(....)' 'PART' 'mapsta.inp'
$
$ Read the status map from file ( FROM != PART ) --------------------- $
$
$ 3 3 3 3 3 3 3 3 3 3 3 3
$ 3 2 1 1 1 1 0 1 1 1 1 3
$ 3 2 1 1 1 1 0 1 1 1 1 3
$ 3 2 1 1 1 1 0 1 1 1 1 3
$ 3 2 1 1 1 1 0 0 1 1 1 3
$ 3 2 1 1 1 1 1 1 1 1 1 3
$ 3 2 1 1 1 1 1 1 1 1 1 3
$ 3 2 1 1 1 1 1 1 1 1 1 3
$ 3 2 1 1 1 1 1 1 1 1 1 3
$ 3 2 1 1 1 1 1 1 1 1 1 3
$ 3 2 1 1 1 1 1 1 1 1 1 3
$ 3 3 3 3 3 3 3 3 3 3 3 3
$
$ The legend for the input map is :
$
$ 0 : Land point.
$ 1 : Regular sea point.
$ 2 : Active boundary point.
$ 3 : Point excluded from grid.
$
$ Input boundary points from segment data ( FROM = PART ) ------------ $
$ An unlimited number of lines identifying points at which input
$ boundary conditions are to be defined. If the actual input data is
$ not defined in the actual wave model run, the initial conditions
$ will be applied as constant boundary conditions. Each line contains:
$ Discrete grid counters (IX,IY) of the active point and a
$ connect flag. If this flag is true, and the present and previous
$ point are on a grid line or diagonal, all intermediate points
$ are also defined as boundary points.
$
2 2 F
2 11 T
$
$ Close list by defining point (0,0) (mandatory)
$
0 0 F
$
$ Excluded grid points from segment data ( FROM != PART )
$ First defined as lines, identical to the definition of the input
$ boundary points, and closed the same way.
$
0 0 F
$
$ Second, define a point in a closed body of sea points to remove
$ the entire body of sea points. Also close by point (0,0)
$
0 0
$
$ Sedimentary bottom map if namelist &SBT4 SEDMAPD50 = T
$
$ 22 1. 1 1 '(f10.6)' 'NAME' 'SED.txt'
$
$ Output boundary points --------------------------------------------- $
$ Output boundary points are defined as a number of straight lines,
$ defined by its starting point (X0,Y0), increments (DX,DY) and number
$ of points. A negative number of points starts a new output file.
$ Note that this data is only generated if requested by the actual
$ program. Example again for spherical grid in degrees. Note, these do
$ not need to be defined for data transfer between grids in the multi
$ grid driver.
$
1.75 1.50 0.25 -0.10 3
2.25 1.50 -0.10 0.00 -6
0.10 0.10 0.10 0.00 -10
$
$ Close list by defining line with 0 points (mandatory)
$
0. 0. 0. 0. 0
$
$ -------------------------------------------------------------------- $
$ End of input file $
$ -------------------------------------------------------------------- $