User:Jcreer/SMS:ADCIRC Model Control in SMS 13.0

From XMS Wiki
Jump to navigationJump to search


The ADCIRC Model Control dialog is where important project parameters are chosen and defined. The six different tabs are listed below.

  • General Parameters – Contains such options as model type, cold/hot start, terms, maximum number of iterations, bottom stress/friction, etc.
  • Model Formulation
  • Timing – Defines simulation run time, start time, and time step options
  • Output
  • Wind – Enables specific wind types and their options
  • Nodal Attributes
  • Harmonics Analysis – Defines tidal constituents and harmonic analysis options
  • Wind – Enables specific wind types and their options

General Tab

ADCIRC Model Control dialog showing the General Parameters tab

The General parameters tab contains the following:

  • Simulation description –
    • Project title – Title for the project run. 24 character limit.
    • Run ID – ID for the project run. 24 character limit.
  • Run options
    • Hot start file (initial conditions) – This option allows choosing a file to initialize the values of water level and velocity for a simulation. The files must have been saved from a previous run using the same finite element grid. Default is "None: cold start". Other options are for specific file types containing the initial conditions for the model. Select the desired file type and the Select button will appear. Using the Select button will bring up a browser window to select the hot start file.
    • Run parallel ADCIRC
  • Bottom stress/friction – NOLIBF – Controls the type of bottom stress parameterization used in a 2D ADCIRC run.
    • "Linear bottom friction "
      • Friction coefficient – The "Linear" option should only be used for analytical cases when verifying the code. Linear friction causes overdamping in deep water. For "Constant Quadratic" the value of TAU0 should be set based on the principle depth. A value of TAU0 = 0.005 is suggested for deeper water (greater than 10 m depth). For shallow water, TAU0 = 0.02 is recommended. If the domain includes both shallow and deep water, consider the varying quadratic option.
    • "Quadratic bottom friction "
      • Friction coefficient – The "Linear" option should only be used for analytical cases when verifying the code. Linear friction causes overdamping in deep water. For "Constant Quadratic" the value of TAU0 should be set based on the principle depth. A value of TAU0 = 0.005 is suggested for deeper water (greater than 10 m depth). For shallow water, TAU0 = 0.02 is recommended. If the domain includes both shallow and deep water, consider the varying quadratic option.
    • "Hybrid nonlinear bottom friction" – In deep water, the friction coefficient is constant and a quadratic bottom friction law results. In shallow water the friction coefficient increases as the depth decreases.
      • Friction coefficient
      • Break depth
      • Approach factor
      • Increase factor
    • "Quadratic friction nodal attributes – Fric"
      • Select
    • "Mannings N nodal attributes – ManningsN"
      • Select
    • "Chezy C nodal attribute – ChezyFric"
      • Select
  • Coordinates
    • Coriolis option – Constant, Variable. The type of Coriolis forcing used does not impact run time for ADCIRC at all. Variable should be used when the coordinate system is geographic. In this case, ADCIRC will load a Coriolis coefficient for each node based on the latitude of the node. If wanting to run ADCIRC in rectilinear coordinates, specify a constant Coriolis option because the model does not support reprojection from rectilinear to geographic to determine individual latitudes for each node.
    • Calculate SLAM0 and SFEA0 on export
  • Flow boundaries
    • Minimum angle for tangential flow
  • Logging
    • Auto-correct (non-fatal erro override) – NFOVER
    • Logging level
    • Log message output (negative values output to file)

____

  • Initial Values – Cold Start, Hot Start 1, Hot Start2. This option allows choosing a file to initialize the values of water level and velocity for a simulation. The files must have been saved from a previous run using the same finite element grid.
  • Minimum Angle For Tangential Flow
  • Terms – Finite Amplitude Terms, Advective Terms, Time Derivative Terms
    • Wetting/Drying – Activating this option will cause an Options button to appear. Clicking this button will bring up the Wetting/Drying Parameters dialog.
  • Solver Type – Lumped, Direct Banded, Iterative JCG – The Iterative JCG is the default and recommended solver. The other options are in the interface for backward compatibility and analytic tests.
  • Absolute Convergence Criteria
  • Maximum Number of Iterations per Time Step
  • Generalized Properties – Wave Continuity, Lateral Viscosity

Model Formulation Tab

ADCIRC Model Control dialog showing the Model Formulation tab
  • Model Type – IM
    • Drop-down
      • "Barotropic 2DDI" – Uses new GWCE and momentum equation formulations.
      • "Barotropic 3D" – Uses new GWCE and velocity based momentum equations.
      • "Baroclinic 3D" – Uses new GWCE and velocity based momentum equations.
      • "Barotropic 2DDI using lumped GWCE (explicit)" – Used to run ADCIRC in lumped + explicit mode, thereby bypassing the iterative solver.
      • "Barotropic 3D using lumped GWCE (explicit)" – Used to run ADCIRC in lumped + explicit mode, thereby bypassing the iterative solver.
    • Form (IDEN) drop-down – Only appears when "Baroclinic 3D" is selected from the above drop-down:
      • "Diagnostic w/ salinity and temperature"
      • "Diagnostic w/ temperature"
      • "Diagnostic w/ salinity"
      • "Diagnostic w/ sigma T"
      • "Prognostic w/ salinity and temperature"
      • "Prognostic w/ temperature"
      • "Prognostic w/ salinity"
      • "Prognostic w/ sigma T"
    • Spatially varying average horizontal eddy viscosity in sea water wrt depth (EVC)
    • Horizontal viscosity (ESLM) – Spatially constant horizontal eddy viscosity for the momentum equations (units of length2/time).
  • Nonlinear terms
    • Finite amplitude terms (NOLIFA) – Parameter controlling the finite amplitude terms in ADCIRC.
      • None – Finite amplitude terms are not included in the model run, and wetting and drying of elements is disabled.
      • No wetting/drying – Finite amplitude terms are included in the model run and wetting and drying of elements is disabled.
      • With wetting/drying – Finite amplitude terms ARE included in the model run and wetting and drying of elements is enabled.
    • Minimum depth (H0) – Minimum water depth.
    • Minimum vel – option only appears when "With wetting/drying" is selected
    • Advective terms (NOLICA) – Parameter controlling the advective terms in ADCIRC.
    • Time derivative terms (NOLICAT) – Parameter controlling the time derivative portion of the advective terms that occurs in the GWCE form of the continuity equation in ADCIRC. Automatically on when "No wetting/drying" or "With wetting/drying" selected, can be turned on or off otherwise.
    • Advection state nodal attribute
      • If turned on, allows selection of a dataset via a Select button that brings up a Dataset dialog. The Dataset dialog has three buttons:
        • Create – Brings up the Dataset Toolbox dialog to allow creation of a new dataset.
        • Select – Brings up the Select Dataset dialog to allow selection of an existing dataset.
        • Unselect – Clears the selection displayed in the Name field in the Dataset dialog.
  • Generalized wave continuity equation – GWCE
    • Weighting factor (TAU0) drop-down:
      • "Pure wave equation: TAU0 = 0" – The GWCE is a pure wave equation.
      • "Pure primitive equation: TAU0 specified" – The GWCE behaves like a pure primitive continuity equation.
        • TAU0 – Allows manual entry of the TAU0 factor.
      • "Spatially discrete/time constant: TAU0 = -1" – The TAU0 is spatially varying but constant in time.
      • "Spatially variable/time constant: TAU0 = -2" – The TAU0 is spatially varying but constant in time.
      • "From nodal attribute: TAU0 = -3" – The TAU0 varies spatially and in time,
        • If turned on, allows selection of a dataset via a Select button that brings up a Dataset dialog. The Dataset dialog has three buttons:
          • Create – Brings up the Dataset Toolbox dialog to allow creation of a new dataset.
          • Select – Brings up the Select Dataset dialog to allow selection of an existing dataset.
          • Unselect – Clears the selection displayed in the Name field in the Dataset dialog.
      • "From local friction: TAU0 = -5" – The TAU0 varies spatially and in time, and is dependent on the local friction.
        • Tau0FullDomainMin and Tau0FullDomainMax – Specified values that the spatially and time varying TAU0 scheme must stay between.
    • Time weighting factors – These options are not editable when either of the explicit model type options are selected.
      • k+1 (A00) – Time weighting factor in the GWCE.
      • k (B00) – Time weighting factor in the GWCE.
      • k-1 (C00) – Time weighting factor in the GWCE.
  • Solver option (ITITER) drop-down – Specifies the solver that will be used for the GWCE.
    • "Lumped"
    • "JCG solver"
  • Solver message level (ISLDIA) drop-down – Determines which messages will be displayed.
    • "Fatal" – Fatal error messgs only from ITPACKV 2D.
    • "Warnings and minimum" – Warning messgs and minimum output from ITPACKV 2D.
    • "Algorithm progress" – Reasonable summary of algorithm progress from ITPACKV 2D.
    • "Parameters and info" – Parameter values and informative comments from ITPACKV 2D.
    • "Approx solution/iter" – Approximate solution after each iteration from ITPACKV 2D.
    • "Original system" – Original system from ITPACKV 2D.
  • Absolute convergence (CONVCR) – Absolute convergence criteria (should be no smaller than 500 times the machine precision).
  • Maximum number of iterations (ITMAX) – Maximum number of iterations each time step.

Timing Tab

ADCIRC Model Control dialog showing the '"Timing tab

The Timing tab contains the following options:

  • Interpolation reference date – Allows selection of a start time. If using a hot start file, this time much exactly match the start time of the hot start file.
  • Timing
    • Time step (seconds) (DTDP) – ADCIRC time step (in seconds).
    • Length of run (RUNDAY) –
  • Ramp Options – NRAMP – Controls whether a ramp is applied to ADCIRC forcing functions.
    • Number of hyperpoic tangent spin up ramps drop-down:
      • "0" – No ramp function is used with forcing functions.
      • "1" – A hyperbolic tangent ramp function is specified and applied to forcing from surface elevation specified boundary conditions, nonzero flux boundary conditions, tidal potential, wind and atmospheric pressure and wave radiation stress.
      • "2" – A hyperbolic tangent ramp function is specified and applied to forcing from surface elevation specified boundary conditions, tidal potential, wind and atmospheric pressure and wave radiation stress. A separate hyperbolic tangent ramp function is specified and applied to the nonzero external flux boundary conditions.
      • "3" – A hyperbolic tangent ramp function is specified and applied to forcing from tidal potential, wind and atmospheric pressure and wave radiation stress. Two additional hyperbolic tangent ramp functions are specified and applied to the nonzero external flux boundary conditions and the nonzero internal flux boundary conditions.
      • "4" – A hyperbolic tangent ramp function is specified and applied to forcing from tidal potential, wind and atmospheric pressure and wave radiation stress. Two additional hyperbolic tangent ramp functions are specified and applied to the nonzero external flux boundary conditions, the nonzero internal flux boundary conditions and the surface elevation specified boundary conditions.
      • "5" – A hyperbolic tangent ramp function is specified and applied to forcing from wind and atmospheric pressure and wave radiation stress. Four additional hyperbolic tangent ramp functions are specified and applied to the nonzero external flux boundary conditions and the nonzero internal flux boundary conditions, the surface elevation specified boundary conditions and to the tidal potential.
      • "6" – A hyperbolic tangent ramp function is specified and applied to forcing from wave radiation stress. Five additional hyperbolic tangent ramp functions are specified and applied to the nonzero external flux boundary conditions, the nonzero internal flux boundary conditions, the surface elevation specified boundary conditions, to the tidal potential and the wind and atmospheric pressure.
      • "7" – A general hyperbolic tangent ramp function is specified, and six additional hyperbolic tangent ramp functions are specified and applied to the nonzero external flux boundary conditions, the nonzero internal flux boundary conditions, the surface elevation specified boundary conditions, to the tidal potential, the wind and atmospheric pressure, and the wave radiation stress.
      • "8" – A general hyperbolic tangent ramp function is specified, and seven additional hyperbolic tangent ramp functions are specified and applied to the nonzero external flux boundary conditions, the nonzero internal flux boundary conditions, the surface elevation specified boundary conditions, to the tidal potential, the wind and atmospheric pressure, the wave radiation stress, and the meteorological ramp delay.

The following options appear when one or more spin up ramps are selected from the above drop-down:

  • Duration (days) (all of the options below are values in decimal days):
    • General Ramp (DRAMP) – Value used to compute the ramp function that ramps up ADCIRC forcings from zero. Available if NRAMP drop-down is set to "1" or higher.
    • External flux ramp (DRAMPExtFlux) – Value used to compute the ramp function that ramps up the nonzero external flux boundary condition. Available if NRAMP drop-down is set to "2" or higher.
    • Flux setting time (FluxSettingTime) – Time that it takes for the river flux boundary condition and the river bottom friction to equilibrate so the water surface elevation can find its steady state. Available if NRAMP drop-down is set to "2" or higher.
    • Internal flux ramp (DRAMPIntFlux) – Value used to compute the ramp function that ramps up the nonzero internal flux boundary condition. Available if NRAMP drop-down is set to "3" or higher.
    • Specified WSE ramp (DRAMPElev) – Value used to compute the ramp function that ramps up the elevation-specified boundary condition. Available if NRAMP drop-down is set to "4" or higher.
    • Tidal potential ramp (DRAMPTip) – Value used to compute the ramp function that ramps up the tidal potential. Available if NRAMP drop-down is set to "5" or higher.
    • Wind/Pressure ramp (DRAMPMete) – Value used to compute the ramp function that ramps up the wind and atmospheric pressure. Available if NRAMP drop-down is set to "6" or higher.
    • Wave stress ramp (DRAMPWRad) – Value used to compute the ramp function that ramps up the wave radiation stress. Available if NRAMP drop-down is set to "7" or higher.
    • Delay (DUnRampMete) – The meteorological ramp delay parameter. It delays the application of the meteorological ramp for the specified length of time, relative to the ADCIRC cold start time. Available if NRAMP drop-down is set to "8".

Output Tab

ADCIRC Model Control dialog showing the Output tab
  • Global output
    • Water surface elevation
      • Output format
      • Append on hot start
      • Start
      • End
      • Increment
    • Velocity
      • Ouput Format
      • Append on hot start
      • Start
      • End
      • Increment
    • Meteorological
      • Ouput Format
      • Append on hot start
      • Start
      • End
      • Increment
  • Recording stations
    • Water surface elevation
    • Velocity
    • Meteorological
  • Hot start input
    • Format
    • Increment

Wind Tab

ADCIRC Model Control dialog showing the Wind tab

ADCIRC can model and compute wind velocities and stresses. The ADCIRC model is able to use wind data from a variety of different file types. The SMS interface is not able to read and create all of the available wind data file types. There are three types that SMS is able to read and create. Even though SMS does not read and create these other file types, it does tell ADCIRC that there is data there, where it is, and how to access it for a simulation to run correctly.

One of the parameters found in the fort.22 file (Single File Meteorological Forcing Input File) will tell which of the wind data file types SMS can read and create. The first parameter specified in the fort.22 file is called the Node Wind Stress (NWS) value. This value represents what type of wind data is in the file. It can range from 0 to 100+. SMS can read and create fort.22 files that are of the NWS value of 1, 2, or 5. All other values of the NWS value are only shown and directed to the ADCIRC model by SMS so that a simulation will run correctly.

  • NWS = 0 – No Wind
  • NWS = 1 – Stress/Pressure – every node/every time step
  • NWS = 2 – Stress/Pressure – every node/every time interval
  • NWS = 3 – Velocity in US Navy feet Numeric format
  • NWS = 4 – Velocity/Pressure in PBL/JAG format
  • NWS = 5 – Velocity/Pressure – every node/time interval
  • NWS = 6 – Velocity/Pressure on grid
  • NWS = 8 – Symmetric cyclone storm from path
  • NWS = 10 – 10m Velocity/Pressure NWS AVN files
  • NWS = 11 – 10m Velocity/Pressure NWS ETA 29km files
  • NWS = 12 – Velocity/Pressure OWI files
  • NWS = 15 – HWind files
  • NWS = 16 – GFDL files
  • NWS = 19 – Asymmetric cyclonic storm from path
  • NWS = 20 – GAHM Format cyclonic storm path

Nodal Attributes Tab

  • Surface submergence state
  • Sea surface height above geoid
  • Wave refraction in swan
  • Bridge pilings friction parameters
  • Elemental slope limiter
  • Initial river elevation
ADCIRC Model Control dialog showing the Nodal Attributes tab

Harmonics Analysis Tab

ADCIRC Model Control dialog showing the Harmonics Analysis tab

Tidal constituents are variations in tides that are created by different frequencies of astronomical forcing. They arise due to the gravitational influences of the Moon and Sun on the Earth, the tilt of the rotational axis of the Earth, the elliptical shape of the Moon's orbit around the Earth, the shape of the Earth's orbit around the Sun, and other such factors. Many constituents have been defined and are classified based on their cycle lengths. Most of the tidal constituents used in the ADCIRC interface of SMS are either diurnal (one cycle per day) or semidiurnal (two cycles per day) in nature. Harmonic constituents are those variations that have periods of less than half a day. For a table of commonly used constituents see Principal Tidal Constituents.[1]

SMS uses the following databases to provide the tidal constituent data to an ADCIRC simulation:

  • LeProvost – The LeProvost database is a set of *.legi files which provide amplitude and phase information. The LeProvost database can be downloaded at: sms.aquaveo.com/leprovost.zip.
An image of the LeProvost tidal database domain coverage can be found here – LeProvost Domain.png[2]
  • ADCIRC – The ADCIRC database includes a grid file and a harmonics file that uses an *.exe extractor to put information into SMS. The ADCIRC database can be downloaded at: sms.aquaveo.com/adcirctides.zip. More information about the ADCIRC database can be obtained at the ADCIRC web page: ADCIRC Tides Databases. The ADCIRC database is only valid for projections in North America.

The path to each database on a computer must be specified in the Preferences dialog by going to Edit | Preferences then selecting the File Locations tab. The two tidal constituent databases (LeProvost and ADCIRC) are only accessible if the SMS project has a global (not local) display projection.

Advanced Controls

  • Changing bathymetry
    • Use hostrarted fort.141
    • Domain
    • Update increment
    • Increment duration
ADCIRC Model Control dialog showing the Advanced Controls tab

References

  1. ^ Cite error: Invalid <ref> tag; no text was provided for refs named ADCIRC article
  2. ^ Le Provost, C., Genco, M. L., and Lyard, F. (1995). “Modeling and predicting tides over the World Ocean,” Quantitative Skill Assessment for Coastal Ocean Models, Coastal and Estuaring Studies, Vol. 47, pp 175-201.

Related Links


v • e
SMS – Surface-water Modeling System
Modules:1D Grid  •   Cartesian Grid  •   Curvilinear Grid  •   GIS  •   Map  •   Mesh  •  Particle  •   Quadtree  •  Raster  •   Scatter   •   UGridSMS.png
General Models:3D Structure  •   FVCOM  •   Generic  •   PTM
Coastal Models:ADCIRC   •   BOUSS-2D  •   CGWAVE  •   CMS-Flow   •   CMS-Wave  •   GenCade  •   STWAVE  •   WAM
Riverine/Estuarine Models:AdH   •   HEC-RAS   •   HYDRO AS-2D   •   RMA2   •   RMA4   •   SRH-2D   •   TUFLOW   •   TUFLOW FV
Aquaveo   •   SMS Tutorials   •   SMS Workflows
Retrieved from "https://www.xmswiki.com/index.php?title=User:Jcreer/SMS:ADCIRC_Model_Control_in_SMS_13.0&oldid=137615"