CMS-Wave Model Parameters: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
Line 25: Line 25:
! Number !! Variable !! Argument Type !! Range !! Description
! 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
|  1 || iprp ||  INTEGER || 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
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  
| 2 || icur || INTEGER  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 3 || ibk  || INTEGER  || 0 - no wave breaking output, 1 - output breaking indices, 2 - output energy dissipation rate || Wave breaking output option
|-
| 4 || irs || INTEGER  || 0 - no wave radiation stress calculation or output, 1 - calculate and output radiation stresses, 2 - calculate and output radiation stresses plus setup/max-water-level|| Current interaction
|-
| 5 || kout || INTEGER || >= 0 || Number of special wave output location, output spectrum in *.obs
and parameters in selhts.out
|-
| 6 || ibnd || INTEGER ||0 - no input a parent spectrum *.nst, 1 - read *.nst, averaging input spectrum,
2 - read *.nst, spatially variable spectrum input || Nesting option.
|-
| 7 || iwet|| INTEGER || 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 || Wetting and drying options.
|-
| 8 || ibf || INTEGER || 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 || Bottom friction option.
|-
| 9 || iark  || INTEGER ||  0 - without forward reflection, 1 - with forward reflection || Forward reflection option.
|-
| 10 || iarkr || INTEGER || 0 - without backward reflection, 1 - with backward reflection || backward reflection option.
|-
| 11 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 12 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 13 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 14 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 15 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 16 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 17 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 18 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 19 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 20 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 21 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 22 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 23 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 24 || icur || REAL  || 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data || Current interaction  
|}
|}


1<sup>st</sup>      2<sup>nd</sup>    3<sup>rd</sup>    4<sup>th</sup>    5<sup>th</sup>     
6<sup>th </sup>    7<sup>th</sup>      8<sup>th</sup>    9<sup>th</sup>   
10<sup>th</sup>    11<sup>th</sup>    12<sup>th</sup>  13<sup>th</sup> 
14<sup>th</sup>  15<sup>th</sup>
iprp  icur  ibk    irs    kout  ibnd  iwet  ibf    iark  iarkr  akap  bf    ark    arkr  iwvbk
16<sup>th</sup>      17<sup>th</sup>    18<sup>th</sup>      19<sup>th</sup>     
20<sup>th</sup>    21<sup>st</sup>    22<sup>nd</sup>      23<sup>rd</sup>   
24<sup>th</sup>
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  
Among these 24 parameters in *.std, the first 6 parameters are always required in CMS-Wave and  
Line 46: Line 82:
default values if not provided in the *.std.  The more specific use and options associated with  
default values if not provided in the *.std.  The more specific use and options associated with  
each of these 1<sup>st</sup> to 24<sup>th</sup> parameters are given below.
each of these 1<sup>st</sup> to 24<sup>th</sup> 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)
akap  = 0 to 4 (diffraction intensity, 0 for zero diffraction, 4 for strong diffraction)



Revision as of 17:53, 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 INTEGER 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 INTEGER 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
3 ibk INTEGER 0 - no wave breaking output, 1 - output breaking indices, 2 - output energy dissipation rate Wave breaking output option
4 irs INTEGER 0 - no wave radiation stress calculation or output, 1 - calculate and output radiation stresses, 2 - calculate and output radiation stresses plus setup/max-water-level Current interaction
5 kout INTEGER >= 0 Number of special wave output location, output spectrum in *.obs

and parameters in selhts.out

6 ibnd INTEGER 0 - no input a parent spectrum *.nst, 1 - read *.nst, averaging input spectrum,

2 - read *.nst, spatially variable spectrum input || Nesting option.

7 iwet INTEGER 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 Wetting and drying options.
8 ibf INTEGER 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 Bottom friction option.
9 iark INTEGER 0 - without forward reflection, 1 - with forward reflection Forward reflection option.
10 iarkr INTEGER 0 - without backward reflection, 1 - with backward reflection backward reflection option.
11 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
12 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
13 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
14 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
15 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
16 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
17 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
18 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
19 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
20 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
21 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
22 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
23 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction
24 icur REAL 0 - no current input, 1 - with current input *.cur, 2 -with *.cur, use only the 1st set current data Current interaction


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. 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).