CMS-Wave Model Parameters: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
 
(55 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 =
[[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}}
 
where 
 
<math>\rho_s</math>  = real and positive power parameter [-]
 
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> .
 
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]]]
 
== 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 =
[[File:CMS Wave Model Control OutputControl. 13.2.12.png|500px|thumb|center]]
== Limit observation output ==
== Radiation stresses ==
== Sea/swell ==
==Breaking type==
 
= Options =
[[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]==
 
==Infragravity wave effect==
==[https://cirpwiki.info/wiki/CMS-Wave:Diffraction 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:
To setup the model parameters for CMS-Wave:


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


2. Users can also specify constant or varied forward and backward reflection
2. Users can also specify constant or varied forward and backward reflection
coefficients in Settings,
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


<br style="clear:both" />
<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]


= CMS-Wave Model Control File =
= CMS-Wave Model Control File formats=
The *.std has a maximum of 24 parameters - the first 15 parameters are more the basic ones as
The CMS-Wave Options File (*.std) can have one of 3 input formats. Click one of the options below for the format description.
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'''
{|
{| border="1"
! Number !! Variable !! Argument Type !! Options/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 <br/> 1 - read *.nst, averaging input spectrum <br/> 2 - read *.nst, spatially variable spectrum input || Nesting option.
|-
| 7 || iwet|| INTEGER || 0 - allow wet/dry, default <br/> 1 - without wet/dry <br/> -1 allow wet/dry, output swell and local sea files <br/> -2 - output combined steering wav files <br/> -3 - output swell, local sea, and combined wav files || Wetting and drying options.
|-
| 8 || ibf || INTEGER || 0 - no bottom friction calc <br/> 1 - constant Darcy-Weisbach coef, c_f <br/> 2 -read variable c_f file, *.fric <br/> 3 - constant Mannings n <br/> 4 - read variable Mannings n file, *.fric || Bottom friction option.
|-
| 9 || iark  || INTEGER ||  0 - without forward reflection <br/> 1 - with forward reflection || Forward reflection option.
|-
| 10 || iarkr || INTEGER || 0 - without backward reflection, 1 - with backward reflection || backward reflection option.
|-
| 11 || bf || REAL  || >=0 || constant bottom friction coef c_f or n
(typical value is 0.005 for c_f and 0.025 for Mannings n)
|-
| 12 || ark || REAL  || 0.0<=ark<=1.0 || Constant forward reflection coef, global specification (0 for zero reflection, 1 for full reflection).
|-
| 13 || arkr || REAL || 0.0<=arkr<=1.0 || Constant backward reflection coef, global specification (
0 for zero reflection, 1 for full reflection)
|-
| 14 || iwvbk || INTEGER || 0 - Goda-extended <br/> 1 - Miche-extended <br/> 2 - Battjes and Janssen <br/> 3 - Chawla and Kirby) || Option for the primary wave breaking formula.
|-
|-
| 15 || nonln || INTEGER ||  0 - none, default <br/> 1 - nonlinear wave-wave interaction || Nonlinear wave-wave interaction
|1. [[Wave_NoCards|No cards - order of values on one line is very important.]] || '''(Used with SMS 11.1 and previous)'''
|-
|-
| 16 || igrav || INTEGER || 0 - none, default <br/> 1 - infra-gravity wave enter inlets || Infragravity waves option.
|2. [[Wave_CardFormat1|Card-based - [value]  [!card name]]] || '''(Used with SMS 11.2 to 13.1)'''
|-
|-
| 17 || irunup || INTEGER || 0 - none, default <br/> 1 - automatic, runup relative to absolute datum <br/>
|3. [[Wave_CardFormat2|Card-based - [card name] [value(s)]]] || '''(Used with SMS 13.2 and after)'''
2 - automatic, runup relative to updated MWL || Runup option.
|-
| 18 || imud || INTEGER  || 0 - none default <br/> 1 - Mud dissipation on || Mud dissipation option. The kinematic viscosity is specified in mud.dat in units of m<sup>2</sup>/sec
|-
| 19 || iwnd || INTEGER || 0 - none, default <br/> 1 - Spatially variable wind on || Spatially variable wind field option. The winds are specified in wind.dat in units of m/s and in the reference frame of the CMS-Wave grid
|-
| 20 || isolv || INTEGER || 0 - GSR solver,  default <br/>  1 - ADI || Matrix solver for CMS-Wave.
|-
| 21 || ixmdf || INTEGER || 0 - output ascii, default <br/> 1 - output xmdf <br/> 2 - input & output xmdf || XMDF input and output options.
|-
| 22 || iproc || INGEGER  || >=0 || Number of threads for parallel computing. Optimum number is approximately equal to the total row number divided by 300. Only for isolv = 0.
|-
| 23 || iview || INTEGER || 0 - half-plane, default <br/> 1 - full-plane || Half-plane/full-plane option. 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
|-
| 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
[[CMS-Wave_File_Formats | Back to CMS-Wave File Formats]]
the remaining ones starting any parameter after the 6<sup>th</sup> 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<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)
 
 
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).

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