CMS-Wave Model Parameters: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
Line 16: Line 16:
<br style="clear:both" />
<br style="clear:both" />


= CMS-Wave Advanced Features =
= CMS-Wave Model Control File =
The most recent CMS-Wave code developed is Version 3.2.  Several new capabilities and
advanced features in this version include:
* Full-plane
* Automatic wave run-up calculation
* Infra-gravity wave
* Nonlinear wave-wave interaction
* Muddy bottom
* Binary file output
* Selection of multiple processors
* Permeable structure
* Spatially varied wind input
* Spatially varied spectral input
* Grid nesting
Users can use SMS11 or simply edit the existing model control file *.std to specify some of
these advanced features.
The *.std has a maximum of 24 parameters - the first 15 parameters are more the basic ones as  
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  
described in the CMS-Wave Technical Report (CHL-TR-08-13) while the remaining 9 parameters  

Revision as of 17:37, 16 March 2011

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

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

iprp icur ibk irs kout ibnd iwet ibf iark iarkr akap bf ark arkr iwvbk 16th 17th 18th 19th 20th 21st 22nd 23rd 24th 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 6th will be assigned to the default values if not provided in the *.std. The more specific use and options associated with each of these 1st to 24th 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 m2/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).