$ -------------------------------------------------------------------- $ $ WAVEWATCH III Grid preprocessor input file $ $ -------------------------------------------------------------------- $ $ Grid name (C*30, in quotes) $ 'REF-UG' $ $ 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.0373 32 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 T 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). $ $ 600. 200. 200. 10. $ $ 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. $ $ &UNST UGOBCAUTO = F , EXPFSN = T, UGOBCFILE='../input/ref1.mshb' / $ $ Bottom friction - - - - - - - - - - - - - - - - - - - - - - - - - - $ JONSWAP : Namelist SBT1 $ GAMMA : As it says. $ &SBT1 GAMMA = 0.15 / $ $ Propagation schemes ------------------------------------------------ $ $ First order : Namelist PRO1 $ CFLTM : Maximum CFL number for refraction. $ $ UQ with diffusion : Namelist PRO2 $ CFLTM : Maximum CFL number for refraction. $ FLSOFT : Flag for 'soft' land boundaries. $ 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 with averaging : Namelist PRO3 $ CFLTM : Maximum CFL number for refraction. $ FLSOFT : Flag for 'soft' land boundaries. $ WDTHCG : Tuning factor propag. direction. $ WDTHTH : Tuning factor normal direction. $ $ UQ with divergence : Namelist PRO4 $ CFLTM : Maximum CFL number for refraction. $ FLSOFT : Flag for 'soft' land boundaries. $ QTFAC : Tuning factor Eq. (3.41). $ RSFAC : Tuning factor Eq. (3.42). $ RNFAC : Tuning factor Eq. (3.43). $ $ Miscellaneous ------------------------------------------------------ $ $ Misc. parameters : Namelist MISC $ CICE0 : Ice concentration cut-off. $ CICEN : Ice concentration cut-off. $ 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. &REF1 REFCOAST=0.1 / $ $ In the 'Out of the box' test setup we run with sub-grid obstacles $ and with continuous ice treatment. $ $ $ Mandatory string to identify end of namelist input section. $ END OF NAMELISTS $ $ FLAG for grid features $ 1 Type of grid 'UNST' 'RECT' 'CURV' $ 2 Flag for geographical coordinates (LLG) $ 3 Flag for periodic grid $ $ Define grid -------------------------------------------------------- $ $ Four records containing : $ 1 NX, NY. As the outer grid lines are always defined as land $ points, the minimum size is 3x3. $ 2 Grid increments SX, SY (degr.or m) and scaling (division) factor. $ If NX*SX is 360., latitudinal closure is applied. $ 3 Coordinates of (1,1) (degr.) and scaling (division) factor. $ 4 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. $ $ Example for longitude-latitude grid (switch !/LLG), for Cartesian $ grid the unit is meters (NOT km). $ $ 'UNST' T F $ 4.0 0.30 20 -1. 4 1 '(20f10.2)' 'NAME' './../input/ref1.msh' $ $ If the above unit number equals 10, the bottom data is read from $ this file and follows below (no intermediate comment lines allowed). $ $ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 $ $ If sub-grid information is avalaible 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 bottem 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 *** $ 10 3 1 '(....)' 'PART' 'mapsta.inp' $ Input boundary points ---------------------------------------------- $ $ 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. $ 1 1 F 19 1 T $ 70 1 F 53 1 F 107 1 F 62 1 F $ 41 1 F 110 1 F 90 1 F 88 1 F $ 86 1 F 87 1 F 79 1 F 109 1 F 100 1 F $ 101 1 F 105 1 F 50 1 F 42 1 F $ $ 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 os sea points. Also close by point (0,0) $ 0 0 $ $ 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. $ $ -2.5312 48.5 0.00 0.008738 102 $ -2.5312 49.3850 0.013554 0.00 51 $ $ Close list by defining line with 0 points (mandatory) $ 0. 0. 0. 0. 0 $ $ -------------------------------------------------------------------- $ $ End of input file $ $ -------------------------------------------------------------------- $