CMS-Wave Model Parameters

Model Control

CMS-Wave Model Control window.

To setup the model parameters for CMS-Wave:

1. Go to CMS-Wave, Model Control, and turn on Allow wetting and drying and Bed friction (Figure 38),

2. Users can also specify constant or varied forward and backward reflection coefficients in Settings,

3. Water level and wind information are optional source as specified under Wave Source in addition to the spectral input data,

4. File, Save As, Wave.sim (selecting the Save As Type as

CMS-Wave Model Control File

The *.std has a maximum of 24 parameters - the first 15 parameters are more the basic ones as described in the CMS-Wave Technical Report (CHL-TR-08-13) while the remaining 9 parameters are relatively new for advanced CMS-Wave features.

Table 1. CMS-Wave parameters in STD file

Number	Variable	Argument Type	Range	Description
1	iprp	REAL	0 - waves and wind input in *.eng, 1 - waves only, neglect wind input in *.eng, -1 - fast mode, 2 - forced grid internal rotation 3 - without lateral energy flux || Wave propagation mode
2	icur	REAL	0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data	Current interaction

1^st 2^nd 3^rd 4^th 5^th 6^th 7^th 8^th 9^th 10^th 11^th 12^th 13^th 14^th 15^th

iprp icur ibk irs kout ibnd iwet ibf iark iarkr akap bf ark arkr iwvbk 16^th 17^th 18^th 19^th 20^th 21^st 22^nd 23^rd 24^th nonln igrav irunup imud iwnd isolv ixmdf iproc iview

Among these 24 parameters in *.std, the first 6 parameters are always required in CMS-Wave and the remaining ones starting any parameter after the 6^th will be assigned to the default values if not provided in the *.std. The more specific use and options associated with each of these 1^st to 24^th parameters are given below.

iprp = 0 (wave propagation with wind input in *.eng) 1 (wave propagation only, neglect wind input in *.eng) -1 (fast mode) 2 (forced grid internal rotation) 3 (without lateral energy flux)

icur = 0 (no current input) 1 (with current input *.cur) 2 (with *.cur, use only the 1st set current data)

ibk = 0 (no wave break info output) 1 (output breaking indices *.brk) 2 (output energy dissipation rate *.brk)

irs = 0 (no wave radiation stress calc) 1 (output radiation stress *.rad) 2 (calculate/output setup/max-water-level + *.rad)

kout = number of special wave output location, output spectrum in *.obs and parameters in selhts.out

ibnd = 0 (no input a parent spectrum *.nst) 1 (read *.nst, averaging input spectrum) 2 (read *.nst, spatially variable spectrum input)

iwet = 0 (allow wet/dry, default) 1 (without wet/dry) -1 (allow wet/dry, output swell and local sea files) -2 (output combined steering wav files) -3 (output swell, local sea, and combined wav files)

ibf = 0 (no bottom friction calc) 1 (constant Darcy-Weisbach coef, c_f) 2 (read variable c_f file, *.fric) 3 (constant Mannings n) 4 (read variable Mannings n file, *.fric)

iark = 0 (without forward reflection) 1 (with forward reflection)

iarkr = 0 (without backward reflection) 1 (with backward reflection)

akap = 0 to 4 (diffraction intensity, 0 for zero diffraction, 4 for strong diffraction)

bf = constant bottom friction coef c_f or n (typical value is 0.005 for c_f and 0.025 for Mannings n)

ark = 0 to 1 (constant forward reflection coef, global specification, 0 for zero reflection, 1 for 100% or fully reflection)

arkr = 0 to 1 (constant backward reflection coef, global specification, 0 for zero reflection, 1 for 100% or fully reflection) iwvbk = 0 to 3 (option for the primary wave breaking formula: 0 for Goda-extended, 1 for Miche-extended, 2 for Battjes and Janssen, 3 for Chawla and Kirby)

nonln = 0 (none, default) 1 (nonlinear wave-wave interaction)

igrav = 0 (none, default) 1 (infra-gravity wave enter inlets)

irunup = 0 (none, default) 1 (automatic, runup relative to absolute datum) 2 (automatic, runup relative to updated MWL)

imud = 0 (mud.dat, default) 1 (none) ---- useful to users who may not want to include the mud effect when the mud.dat exists (typical maximum kinematic viscosity in mud.dat is 0.04 m²/sec)

iwnd = 0 (wind.dat, default) 1 (none) ---- useful in steering if users decide not to read the spatially varied wind field input wind.dat when the wind.dat file exists

isolv = 0 (GSR solver, default) 1 (ADI)

ixmdf = 0 (output ascii, default) 1 (output xmdf) 2 (input & output xmdf)

iproc = 0 (same as 1, default) n (n processors for isolv = 0) optimum n = (total row number) /300

iview = 0 (half-plane, default) 1 (full-plane) --- in the full plane, users can provide additional input wave spectrum file wave.spc (same format as the *.eng) along the opposite side boundary (an imaginary origin for wave.spc at the opposite corner; users can rotate the CMS-Wave grid by 180 deg in SMS to generate this wave.spc)

Figure 6.1 shows the CMS-Wave interface window for Model Control in SMS11.

• Full-plane – In this mode, CMS-Wave performs two half-plane runs in the same grid. The first

run is in the half-plane with the principle wave direction toward the shore. The second run is in the seaward half-plane. Upon the completion of the second run, two half-plane results are combined to one full-plane solution. Because the run time for the full-plane is approximately twice of the regular half-plane, users shall consider the full-plane mode only if the full-plane features like wave generation and propagation in a bay or around an island. An example is to run the Shark River wave case, 2009.sim, in the full plane (modify 2009.std).