CMS-Wave Model Parameters: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
 
(57 intermediate revisions by 3 users not shown)
Line 1: Line 1:
=Model Control=
[[File:CMSWave Model Parameters GeneralTab V13.2.12.png|thumb|right|400px|CMS-Wave Model Control window Version 13.2.12.]]
[[Image:Shark_Fig38.PNG|thumb|right|400px|CMS-Wave Model Control window.]]
= Parameters =
To setup the model parameters for CMS-Wave:
[[File:CMS Wave Model Control V. 13.2.12.png|500px|thumb|center]]
==CMS Wave plane mode ==
 
===Full half plane with reverse spectra===
===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 as opposed to the first run.  Upon the completion of the second run, two half-plane results are combined to one full-plane solution (Lin et al. 2012).  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.
===Half Plane===
 
== Source terms ==
CMS-Wave is a phase-averaged model for propagation of directional irre-gular waves over complicated bathymetry and nearshore where wave re-fraction, diffraction, reflection, shoaling, and breaking simultaneously act at inlets. Wave diffraction terms are included in the governing equations following the method of Mase et al. (2005). Four different depth-limiting wave breaking formulas can be selected as options including the interaction with a current. The wave-current interaction is calculated based on the dispersion relationship including wave blocking by an opposing current (Larson and Kraus 2002). Wave generation and whitecapping dissipation are based on the parameterization source term and calibration using field data (Lin and Lin 2004a and b, 2006a). Bottom friction loss is estimated from the classical drag law formula (Collins 1972).
 
== [https://cirpwiki.info/wiki/CMS-Flow:Wave-current_Interaction Current interaction] ==
 
== Bottom friction ==
Bottom friction could be assigned as constant or by a dataset inside SMS. To see more details, see the [https://cirpwiki.info/wiki/CMS-Flow_Model_Parameters#Bottom_and_wall_friction Bottom and wall friction of the CMS Flow].
Usually, the same roughness method used on CMS Wave is used on the CMS Flow parameters bottom friction.
For a detail explanation of how Bottom Friction variable is used on the CMS Wave model, see [https://cirpwiki.info/wiki/Bottom_Friction Bottom Friction]
 
==Surge Fields==
==Wind Fields==
If the spatial wind field input is required, users shall prepare a wind.dat file or *.wind (in the same format as *.cur) to provide the x- and y-component wind data corresponding to the incident wave conditions in the model grid.
 
==Matrix Solver==
https://cirpwiki.info/wiki/CMSFlow_Matrix_Solver
 
= Boundary control =
[[File:CMS Wave Model Control Boundray ControlV. 13.2.12.png|500px|thumb|center]]
Water level and wind information are optional source as specified under
Wave Source in addition to the spectral input data.
==Source==
Spatially varied spectral input – This is simply the case as in a child grid that spatially varied wave spectra are permitted to assign at user specified locations along or near the seaward boundary of the child grid. To apply spatially varied spectra for wave input without a parent grid, users will need to prepare the wave input file with the format as described in the child grid run.
 
https://cirpwiki.info/wiki/CMS-Wave_Input_Spectra
 
==Interpolation==
*Inverse distance weighting
 
The inverse-distance interpolation also referred to as Shepard interpola-tion is given by (Shepard 1968)
 
{{Equation|<math>\phi(\overrightarrow{x}) = \sum_{i=1}^N w_i \phi_i</math>|17-1}}
 
where the interpolation weights are given by
 
{{Equation|<math>w_i = \frac{d_i ^{-\rho_s}}{ \sum_{j} d_j^{-\rho_s}}</math>|17-2}}


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


2. Users can also specify constant or varied forward and backward reflection
<math>\rho_s</math>  = real and positive power parameter [-]
coefficients in Settings,


3. Water level and wind information are optional source as specified under
d = distance between the known points <math>\overrightarrow{x_i}</math>  and the unknown interpolation points <math>\overrightarrow{x}</math>  equal to the Euclidean norm <math>d = ||\overrightarrow{x} - \overrightarrow{x}_i||</math> .
Wave Source in addition to the spectral input data,


4. File, Save As, Wave.sim (selecting the Save As Type as
In this interpolation, the weight of each point decreases with distance from the interpolated point. One advantage of the inverse-distance interpolation is the interpolation weights are independent of the interpolation function, and therefore only need to be calculated once and can be saved for computational efficiency.[[User_Guide_028|[1]]]


<br style="clear:both" />
== Computational spectral grid==
Spectral waves or wave parameters can be generated for the wave grid forcing, or wind direction and speeds can provide the necessary information for wind- wave generation. Full (directional) spectra can be imported into the SMS for the CMS-Wave, as well as simplified wave parameters (angle, wave height, and period, etc).


= CMS-Wave Model Control File =
==Sides==
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'''
==Case data==
{| border="1"
! Number !! Variable !! Argument Type !! Range !! Description
|-
|  1 || iprp ||  INTEGER || 0 - waves and wind input in *.eng <br /> 1 - waves only, neglect wind input in *.eng <br /> -1 - fast mode<br /> 2 - forced grid internal rotation <br /> 3 - without lateral energy flux || Wave propagation mode.
|-
| 2 || icur || INTEGER  || 0 - no current input <br /> 1 - with current input *.cur <br /> 2 -with *.cur, use only the 1st set current data || Current interaction
|-
| 3 || ibk  || INTEGER  || 0 - no wave breaking output <br /> 1 - output breaking indices <br /> 2 - output energy dissipation rate || Wave breaking output option
|-
| 4 || irs || INTEGER  || 0 - no wave radiation stress calculation or output <br /> 1 - calculate and output radiation stresses <br/> 2 - calculate and output radiation stresses plus setup/max-water-level|| Radiation stress and runup options.
|-
| 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
|}


===Wind direction angle conversion===
*Cartesian
*Meteorologic
*Oceanographic
*Shore normal


Among these 24 parameters in *.std, the first 6 parameters are always required in CMS-Wave and
===Populate from Spectra===
the remaining ones starting any parameter after the 6<sup>th</sup> will be assigned to the
===Set Reference Time===
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.
akap  = 0 to 4 (diffraction intensity, 0 for zero diffraction, 4 for strong diffraction)


bf    = constant bottom friction coef c_f or n
= Output control =
(typical value is 0.005 for c_f and 0.025 for Mannings n)
[[File:CMS Wave Model Control OutputControl. 13.2.12.png|500px|thumb|center]]
== Limit observation output ==
== Radiation stresses ==
== Sea/swell ==
==Breaking type==


ark  = 0 to 1 (constant forward reflection coef, global specification,
= Options =
0 for zero reflection, 1 for 100% or fully reflection)
[[File:CMS Wave Model Control Options 13.2.12.png|500px|thumb|center]]
To include (trigger) either of wave run-up, infra-gravity wave, nonlinear wave-wave interaction, binary (xmdf or *.h5) output, multiple processors, muddy bed, and spatial wind field input is just a one-click step in the SMS interface. Additional files are required for the muddy bed and spatial wind field input.
==[https://cirpwiki.info/wiki/CMS-Flow_NUmerical_Methods:_Wetting_and_Drying Allow wetting and drying]==


arkr  = 0 to 1 (constant backward reflection coef, global specification,
==Infragravity wave effect==
0 for zero reflection, 1 for 100% or fully reflection)
==[https://cirpwiki.info/wiki/CMS-Wave:Diffraction Diffraction intensity]==
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)
==Nonlinear wave effect==
==Run up==
==Fast-mode run==
==Roller effects==
The wave roller parametric formulation is commonly applied in the wave spectral model to modify breaking wave energy dissipation nearshore to mimic better the surf zone dynamics. The wave roller effect used in CMS-Wave is based on the roller model developed by Zhang et al. (2014).


igrav  = 0 (none, default) 1 (infra-gravity wave enter inlets)
==Forward reflection==
==Backward reflection==
==Muddy bed==
If the muddy bed calculation is required, users shall prepare a mud.dat file or *.mud (in the same format as *.dep) to list the spatial varying max-imum kinematic viscosity for the entire grid (recommended maximum kinematic viscosity for mud is 0.04 m2/sec).
==Wave breaking formula==
==Date format==


irunup = 0 (none, default)  1 (automatic, runup relative to absolute datum)
To setup the model parameters for CMS-Wave:
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
1. Go to CMS-Wave, Model Control, Options and turn on Allow wetting and drying
is 0.04 m<sup>2</sup>/sec)
and Bed friction


iwnd  = 0 (wind.dat, default)  1 (none)  ----  useful in steering if users decide not
2. Users can also specify constant or varied forward and backward reflection
to read the spatially varied wind field input wind.dat when the wind.dat file exists
coefficients in Settings,


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


ixmdf  = 0 (output ascii, default) 1 (output xmdf)  2 (input & output xmdf)
<br style="clear:both" />
To read more about CMS Wave defintions, see [https://cirpwiki.info/wiki/CMS-Wave_Model_Control CMS Wave Model Control Definitions]


iproc  = 0 (same as 1, default)  n (n processors for isolv = 0)
= CMS-Wave Model Control File formats=
optimum n = (total row number) /300
The CMS-Wave Options File (*.std) can have one of 3 input formats. Click one of the options below for the format description.


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)
|1. [[Wave_NoCards|No cards - order of values on one line is very important.]] || '''(Used with SMS 11.1 and previous)'''
|-
|2. [[Wave_CardFormat1|Card-based - [value]  [!card name]]] || '''(Used with SMS 11.2 to 13.1)'''
|-
|3. [[Wave_CardFormat2|Card-based - [card name] [value(s)]]] || '''(Used with SMS 13.2 and after)'''
|}


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
[[CMS-Wave_File_Formats | Back to CMS-Wave File Formats]]
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).

Latest revision as of 19:52, 20 December 2022

CMS-Wave Model Control window Version 13.2.12.

Parameters

CMS Wave Model Control V. 13.2.12.png

CMS Wave plane mode

Full half plane with reverse spectra

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 as opposed to the first run. Upon the completion of the second run, two half-plane results are combined to one full-plane solution (Lin et al. 2012). 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.

Half Plane

Source terms

CMS-Wave is a phase-averaged model for propagation of directional irre-gular waves over complicated bathymetry and nearshore where wave re-fraction, diffraction, reflection, shoaling, and breaking simultaneously act at inlets. Wave diffraction terms are included in the governing equations following the method of Mase et al. (2005). Four different depth-limiting wave breaking formulas can be selected as options including the interaction with a current. The wave-current interaction is calculated based on the dispersion relationship including wave blocking by an opposing current (Larson and Kraus 2002). Wave generation and whitecapping dissipation are based on the parameterization source term and calibration using field data (Lin and Lin 2004a and b, 2006a). Bottom friction loss is estimated from the classical drag law formula (Collins 1972).

Current interaction

Bottom friction

Bottom friction could be assigned as constant or by a dataset inside SMS. To see more details, see the Bottom and wall friction of the CMS Flow. Usually, the same roughness method used on CMS Wave is used on the CMS Flow parameters bottom friction. For a detail explanation of how Bottom Friction variable is used on the CMS Wave model, see Bottom Friction

Surge Fields

Wind Fields

If the spatial wind field input is required, users shall prepare a wind.dat file or *.wind (in the same format as *.cur) to provide the x- and y-component wind data corresponding to the incident wave conditions in the model grid.

Matrix Solver

https://cirpwiki.info/wiki/CMSFlow_Matrix_Solver

Boundary control

CMS Wave Model Control Boundray ControlV. 13.2.12.png

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

Source

Spatially varied spectral input – This is simply the case as in a child grid that spatially varied wave spectra are permitted to assign at user specified locations along or near the seaward boundary of the child grid. To apply spatially varied spectra for wave input without a parent grid, users will need to prepare the wave input file with the format as described in the child grid run.

https://cirpwiki.info/wiki/CMS-Wave_Input_Spectra

Interpolation

  • Inverse distance weighting

The inverse-distance interpolation also referred to as Shepard interpola-tion is given by (Shepard 1968)

  (17-1)

where the interpolation weights are given by

  (17-2)

where

= real and positive power parameter [-]

d = distance between the known points and the unknown interpolation points equal to the Euclidean norm .

In this interpolation, the weight of each point decreases with distance from the interpolated point. One advantage of the inverse-distance interpolation is the interpolation weights are independent of the interpolation function, and therefore only need to be calculated once and can be saved for computational efficiency.[1]

Computational spectral grid

Spectral waves or wave parameters can be generated for the wave grid forcing, or wind direction and speeds can provide the necessary information for wind- wave generation. Full (directional) spectra can be imported into the SMS for the CMS-Wave, as well as simplified wave parameters (angle, wave height, and period, etc).

Sides

Case data

Wind direction angle conversion

  • Cartesian
  • Meteorologic
  • Oceanographic
  • Shore normal

Populate from Spectra

Set Reference Time

Output control

CMS Wave Model Control OutputControl. 13.2.12.png

Limit observation output

Radiation stresses

Sea/swell

Breaking type

Options

CMS Wave Model Control Options 13.2.12.png

To include (trigger) either of wave run-up, infra-gravity wave, nonlinear wave-wave interaction, binary (xmdf or *.h5) output, multiple processors, muddy bed, and spatial wind field input is just a one-click step in the SMS interface. Additional files are required for the muddy bed and spatial wind field input.

Allow wetting and drying

Infragravity wave effect

Diffraction intensity

Nonlinear wave effect

Run up

Fast-mode run

Roller effects

The wave roller parametric formulation is commonly applied in the wave spectral model to modify breaking wave energy dissipation nearshore to mimic better the surf zone dynamics. The wave roller effect used in CMS-Wave is based on the roller model developed by Zhang et al. (2014).

Forward reflection

Backward reflection

Muddy bed

If the muddy bed calculation is required, users shall prepare a mud.dat file or *.mud (in the same format as *.dep) to list the spatial varying max-imum kinematic viscosity for the entire grid (recommended maximum kinematic viscosity for mud is 0.04 m2/sec).

Wave breaking formula

Date format

To setup the model parameters for CMS-Wave:

1. Go to CMS-Wave, Model Control, Options and turn on Allow wetting and drying and Bed friction

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



To read more about CMS Wave defintions, see CMS Wave Model Control Definitions

CMS-Wave Model Control File formats

The CMS-Wave Options File (*.std) can have one of 3 input formats. Click one of the options below for the format description.

1. No cards - order of values on one line is very important. (Used with SMS 11.1 and previous)
2. Card-based - [value] [!card name] (Used with SMS 11.2 to 13.1)
3. Card-based - [card name] [value(s)] (Used with SMS 13.2 and after)


Back to CMS-Wave File Formats

Retrieved from "https://cirpwiki.info/index.php?title=CMS-Wave_Model_Parameters&oldid=13920"