https://wiki.adcirc.org/api.php?action=feedcontributions&user=Wpringle&feedformat=atomADCIRCWiki - User contributions [en]2024-03-29T14:40:03ZUser contributionsMediaWiki 1.38.1https://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=1094Idealized Channel Problem2021-04-08T03:01:02Z<p>Wpringle: </p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2021. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. Environmental Modelling and Software, 105045. https://doi.org/10.1016/j.envsoft.2021.105045</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. https://doi.org/10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. https://doi.org/10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Note that the short 6 hour length of the test is chosen only to limit simulation time for the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite] where the test case been found. Users may extend the simulation length to simulate more of the inundating phase of the incoming wave. <br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1000px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br />
[[File:Channel_Elev.gif|500px|thumb|Elevation time series for the idealized channel problem]] [[File:Channel_Vel.gif|500px|thumb|North-south velocity time series for the idealized channel problem]]<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files. [[Grid_Development_and_Editing#OceanMesh2D|OceanMesh2D]] functions can be used to automatically generate the sponge_generator_layer attribute ([https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/utilities/Calc_Sponge.m Calc_Sponge]) and the input files ([https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/utilities/Make_f5354.m Make_f5354]).<br />
*[[Fort.14_file_format|IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=NWS&diff=1090NWS2020-11-17T19:50:47Z<p>Wpringle: /* Parameter Summary */</p>
<hr />
<div>'''<code>NWS</code>''' is a parameter in the [[fort.15 file]] that selects the meteorological forcing input type. The value on the "[[fort.15_file_format#NWS|NWS line]]" of the fort.15 file also implicitly includes [[#Value Seen in fort.15 File|other parameters]] affecting wave coupling and ice inputs. Further, <code>NWS</code> affects not just the file type and handling of meteorological data, but also changes what the [[Fort.15_file_format#WTIMINC|meteorological parameter line]] (informally, the <code>[[WTIMINC]]</code> line) looks like in the fort.15 file. See the [[supplemental meteorological/wave/ice parameters]] page for information on the format of this line. <br />
<br />
ADCIRC supports a wide range of meteorological input formats, including moving/fixed gridded data in several file formats, tropical cyclone track and parameter data that can be turned into wind/pressure fields via one of several internal vortex models, and direct specification of wind speeds or stresses on nodes. As a result of the great flexibility and importance of this choice, several pages are devoted to the topic. In particular, see also the [[fort.22 file]] and [[wind stress]] pages. <br />
<br />
== Value Seen in fort.15 File ==<br />
In the fort.15 file, what we call "NWS" is actually a combination of several parameters. For example, given a 5-digit value on that line, <br />
-12305 ! TRICKY NWS IMPOSTER<br />
the first two digits (ten-thousands and thousands) tell us the format of ice data <code>[[NCICE]]=12</code>, the 3rd digit (hundreds) tells us the wave coupling mode <code>[[NRS]]=3</code>, and the last two digits (tens and ones) combined with the sign of the entire value tell us the meteorological forcing mode <code>NWS=-5</code>. If the value has only 3 digits then ADCIRC assumes no ice input <code>[[NCICE]]=0</code>, and if it's 2 digits then ADCIRC further assumes no wave coupling <code>[[NRS]]=0</code>. It is often presumed that when one refers to "NWS", one is referring to the 2-digit value, not what is in the fort.15 file, but be mindful of the ambiguity here. <br />
<br />
== Parameter Summary ==<br />
The following table is a summary of possible <code>NWS</code> values, their descriptions, and associated meteorological input files (required and optional).<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! <code>NWS</code> Value<br />
! Short-name<br />
! Description<br />
! Required Input Files<br />
! Optional Input Files<br />
|-<br />
| 1<br />
| Wind stress, every node, every timestep<br />
| Wind stress and atmospheric pressure are read in at all grid nodes at every model time step from the [[fort.22_file_format#NWS = 1 or 101|fort.22 file]]<br />
| [[fort.22_file_format#NWS = 1 or 101|fort.22]]<br />
|<br />
|-<br />
| 2<br />
| Wind stress, every node, every [[WTIMINC]]<br />
| Wind stress and atmospheric pressure are read in at all grid nodes at a time interval that does not equal the model time step from the [[fort.22_file_format#NWS = 2, -2, 102 or -102|fort.22 file]]. Interpolation in time is used to synchronize the wind and pressure information with the model time step. The wind time interval ([[WTIMINC]]) is specified in the [[fort.15_file_format|fort.15 file]].<br />
| [[fort.22_file_format#NWS = 2, -2, 102 or -102|fort.22]]<br />
|<br />
|-<br />
| 3<br />
| US Navy Fleet Numeric<br />
| Wind velocity is read in from a wind file from the [[fort.22_file_format#NWS = 3 or 103, Fleet Numeric Format|fort.22 file]] in US Navy Fleet Numeric format. This information is interpolated in space onto the ADCIRC grid and in time to synchronize the wind and pressure information with the model time step. Garret’s formula is used to compute wind stress from the wind velocity. Several parameters describing the Fleet Numeric wind file must be specified in the [[fort.15_file_format|fort.15 file]].<br />
| [[fort.22_file_format#NWS = 3 or 103, Fleet Numeric Format|fort.22]]<br />
|<br />
|-<br />
| 4<br />
| PBL/JAG<br />
| Wind velocity and atmospheric pressure are read in (PBL/JAG format) at selected ADCIRC grid nodes from the [[fort.22_file_format#NWS = 4, -4, 104 or -104 - PBL Hurricane Model format|fort.22]] file. Interpolation in time is used to synchronize the wind and pressure information with the model time step. Garret’s formula is used to compute wind stress from wind velocity.<br />
| [[fort.22_file_format#NWS = 4, -4, 104 or -104 - PBL Hurricane Model format|fort.22]]<br />
| <br />
|-<br />
| 5<br />
| Wind velocity, every node, every [[WTIMINC]]<br />
| Wind velocity and atmospheric pressure are read in at all grid nodes from the [[fort.22_file_format#NWS = 5, -5, 105, or -105|fort.22]] File. Interpolation in time is used to synchronize the wind and pressure information with the model time step. Garret’s formula is used to compute wind stress from wind velocity.<br />
| [[fort.22_file_format#NWS = 5, -5, 105, or -105|fort.22]]<br />
|<br />
|-<br />
| 6<br />
| wind velocity, rectangular grid, every [[WTIMINC]]<br />
| Meteorological data (U,V,P) is input on a rectangular grid (either in Longitude, Latitude or Cartesian coordinates, consistent with the grid coordinates) and interpolated in space onto the ADCIRC grid. Wind velocity (U,V @ 10 m above the water surface) must be input in units of m/s and surface atmospheric pressure (P) must be input in units of Pascals = Newtons/square meter. The meteorological grid MUST cover the entire ADCIRC mesh; that is, the ADCIRC mesh must be ENTIRELY within the meteorological grid or an error will result.<br />
| [[fort.22_file_format#NWS = 6 or 106|fort.22]]<br />
|<br />
|-<br />
| 7 {{ADC version|version=future|relation=eq}}<br />
| Wind stress, regular grid, every [[WTIMINC]]<br />
| Surface stress and pressure values are read in on a regular grid from the [[fort.22_file_format|fort.22 file]]. Currently, this feature is not supported for parallel runs because adcprep cannot decompose the files. See [https://github.com/adcirc/adcirc-cg/issues/215]. <br />
| [[fort.22_file_format#NWS = 7 or -7|fort.22]]<br />
|<br />
|-<br />
| 8<br />
| Symmetric Holland Vortex<br />
| Wind velocity and atmospheric pressure are calculated at every node on the fly by ADCIRC internally using the Dynamic Holland model.<br />
| [[Fort.22_file_format#NWS_.3D_8|fort.22]]<br />
| <br />
|-<br />
| 10<br />
| NCDC GFS<br />
| Wind velocity (10 m) and atmospheric pressure are read in from a sequence of National Weather Service (NWS) Aviation (AVN) model output files. Each AVN file is assumed to contain data on a Gaussian longitude, latitude grid at a single time.<br />
| [[fort.200_file_format|fort.200, fort.200+N, fort.200+2*N, fort.200+3*N,….,]] where N is the time interval (in hours) between successive meteorological data<br />
| <br />
|-<br />
| 11<br />
| National Weather Service Eta-29 file<br />
| Wind velocity (10 m) and atmospheric pressure are read in from a sequence of stripped down National Weather Service (NWS) ETA 29km model output files<br />
| [[fort.200_file_format|fort.200, fort.201, fort.202, fort.203,….,]]<br />
| <br />
|-<br />
|-<br />
| 12 / -12 <br />
| OWI ASCII, every [[WTIMINC]]<br />
| See [[NWS12]] for details. Wind velocities (U10, V10) and atmospheric sea level pressure (SLP) are provided in the OWI ASCII format on one to three rectangular (lat/lon) grid(s)<br />
| [[Fort.22_file_format#NWS_.3D_12|fort.22]], [[fort.221]], [[fort.222]]<br />
| [[fort.223]], [[fort.224]], [[fort.217]], [[fort.218]]<br />
|-<br />
| 13 {{ADC version|version=55|relation=ge}}<br />
| OWI NetCDF<br />
| See [[NWS13]] for details. Wind velocities (U10, V10) and atmospheric sea level pressure (SLP) fields are provided in the OWI NetCDF format as 1 or more meshgrid overlays stored in netCDF groups, supporting storm following grids on overlay 2 and on.<br />
| default is [[fort.22.nc]], see [[NWS13]]<br />
| <br />
|-<br />
| 14 / -14 {{ADC version|version=55|relation=ge}}<br />
| GRIB2/NetCDF Binary, every [[WTIMINC]]<br />
| Gridded data of wind velocities (U10, V10) and atmospheric sea level pressure (SLP) are provided in GRIB2 (e.g., GFS, CFSv2) or NetCDF (e.g., ERA5, WRF) binary files. Gridded data may be on a standard rectangular lat/lon grid, a [https://en.wikipedia.org/wiki/Gaussian_grid Gaussian grid], or a projected WRF-like grid. Requires that ADCIRC is compiled with DATETIME, NetCDF and if required, GRIB2 flags enabled (the static libraries must be compiled). Will find and read time-snaps based on the reference date, [[NCDATE]] located near or at the bottom of the [[fort.15_file_format|fort.15 file]] taking into account hot-start times etc. If the negative value is used, OWI ASCII (see NWS = 12) meteorology will overwrite the GRIB2/NetCDF meteorology data in the overlap region (except during the "skipping OWI time snap" phase). <br />
| [[fort.22x.grb2 file|fort.221.grb2, fort.222.grb2]] <br/>'''or'''<br/> [[Fort.22_file#NWS_.3D_.C2.B114_Gridded_GRIB2_or_NetCDF_Wind_and_Pressure|fort.22]], [[fort.221.nc]], [[fort.222.nc]]<br />
|<br />
|-<br />
| 15<br />
| HWIND<br />
| Uses data assimilated snapshots of the wind velocity fields of tropical cyclones that are produced by the NOAA Hurricane Research Division (HRD)<br />
| [[fort.22_file_format#NWS = 15|fort.22]]<br />
| Additional HWIND files specified in the [[fort.22_file_format#NWS = 15|fort.22]] file<br />
|-<br />
| 16<br />
| GFDL<br />
| GFDL model output files produced by the Geophysical Fluid Dynamics Laboratory at NOAA. Each ASCII GFDL model output file contains one or more nested grid dataset where the nested grids are allowed to change in time. Coarse grid data is not stored where finer nest data is given.<br />
| [[fort.22_file_format#NWS = 16|fort.22]]<br />
| Additional GFDL files specified in the [[fort.22_file_format#NWS = 16|fort.22]] file<br />
|-<br />
| 19<br />
| Dynamic Asymmetric Model<ref group="note" name="nws19bad">Use of this [[Typical_ADCIRC_Parameter_Selections#Discouraged_Parameter_Selections|is discouraged]].</ref><br />
| Wind velocity and atmospheric pressure are calculated at exact finite element mesh node locations and directly coupled to ADCIRC at every time step using the asymmetric hurricane vortex formulation based on the Holland gradient wind model. The input file is assumed to correspond to the ATCF Best Track/Objective Aid/Wind Radii Format. This option uses the radii at specific wind speeds (34, 50, 64, 100 knots) reported in the four quadrants (NE, SE, SW, NW) of the storm to calculate the radius of maximum winds as a function of the azimuthal angle. Garret’s formula is used to compute wind stress from the wind velocity. This option allows the user to set a value for Rmax and Holland B Parameter. Additionally the user can select the isotachs to be used for each of the 4 quadrants. The utility program aswip_1.0.3.F located in the /wind folder will generate the NWS=19 formatted file from a NWS=9 formatted fort.22 input file.<br />
| [[fort.22_file_format#NWS = 19|fort.22]]<br />
|<br />
|-<br />
| 20<br />
| [[Generalized Asymmetric Holland Model]]<br />
| The Generalized Asymmetric Holland Model (GAHM) provides a set of theoretical and practical improvements over previous parametric meteorological vortex models in ADCIRC. The track file format is similar to that of the older Dynamic Asymmetric Model (NWS = 19) but with 8 additional columns of data.<br />
| [[fort.22_file_format#NWS = 20, Generalized Asymmetric Holland Model (GAHM)|fort.22]]<br />
|<br />
|}<br />
<br />
==Extended NWS with Ice + Waves==<br />
The following presents a summary of the extended <code>NWS</code> values to included ice-coverage and/or wind wave-coupling<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! style="font-weight:bold;" |Meteorological Data Format<br />
! style="font-weight:bold;" |Met. Only<br />
! style="font-weight:bold;" |Met. plus Waves from fort.23<br />
! style="font-weight:bold;" |Met. plus Waves SWAN<br />
! style="font-weight:bold;" |Met. plus Waves STWAVE<br />
! style="font-weight:bold;" |Met. plus Ice Coverage, Waves off<br />
! style="font-weight:bold;" |Met. plus Ice Coverage OWI-like format plus Waves from fort.23<br />
! style="font-weight:bold;" |Met. plus Ice Coverage OWI-like format plus Waves from SWAN<br />
! style="font-weight:bold;" |Met. plus Ice Coverage OWI-like format plus Waves from STWAVE<br />
|-<br />
|none<br />
| 0<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
|-<br />
|wind stress, every node, every timestep<br />
| 1<br />
| 101<br />
| 301<br />
| 401<br />
|<br />
| 12101<br />
| 12301<br />
| 12401<br />
|-<br />
|wind stress, every node, every WTIMINC<br />
| 2<br />
| 102<br />
| 302<br />
| 402<br />
|<br />
| 12102<br />
| 12302<br />
| 12402<br />
|-<br />
|US Navy Fleet Numeric<br />
| 3<br />
| 103<br />
| 303<br />
| 403<br />
|<br />
| 12103<br />
| 12303<br />
| 12403<br />
|-<br />
|PBL/JAG<br />
| 4<br />
| 104<br />
| 304<br />
| 404<br />
|<br />
| 12104<br />
| 12304<br />
| 12404<br />
|-<br />
|wind velocity, every node, every WTIMINC<br />
| 5<br />
| 105<br />
| 305<br />
| 405<br />
|<br />
| 12105<br />
| 12305<br />
| 12405<br />
|-<br />
|wind velocity, rectangular grid, every WTIMINC<br />
| 6<br />
| 106<br />
| 306<br />
| 406<br />
|<br />
| 12106<br />
| 12306<br />
| 12406<br />
|-<br />
|wind stress, regular grid, every WTIMINC<br />
| 7<br />
| 107<br />
| 307<br />
| 407<br />
|<br />
| 12107<br />
| 12307<br />
| 12407<br />
|-<br />
|symmetrc vortex model<br />
| 8<br />
| 108<br />
| 308<br />
| 408<br />
|<br />
| 12108<br />
| 12308<br />
| 12408<br />
|-<br />
|asymmetric vortex model (no longer available)<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
|-<br />
|National Weather Service AVN<br />
| 10<br />
| 110<br />
| 310<br />
| 410<br />
| 10010 (ice in 4th column of AVN file)<br />
| 12110<br />
| 12310<br />
| 12410<br />
|-<br />
|National Weather Service ETA 29km<br />
| 11<br />
| 111<br />
| 311<br />
| 411<br />
|<br />
| 12111<br />
| 12311<br />
| 12411<br />
|-<br />
|Oceanweather Inc (OWI)<br />
| 12<br />
| 112<br />
| 312<br />
| 412<br />
|<br />
| 12112<br />
| 12312<br />
| 12412<br />
|-<br />
|Oceanweather Inc (OWI) NetCDF<br />
| 13<br />
| 113?<br />
| 313?<br />
| 413?<br />
|<br />
| 12113?<br />
| 12313?<br />
| 12413?<br />
|-<br />
|GRIB2/NetCDF <br />
| 14<br />
| 114<br />
| 314<br />
| 414<br />
| 14014 (GRIB2/NetCDF format ice)<br />
| 14114 (GRIB2/NetCDF format ice)<br />
| 14314 (GRIB2/NetCDF format ice)<br />
| 14414 (GRIB2/NetCDF format ice)<br />
|-<br />
|H*Wind<br />
| 15<br />
| 115<br />
| 315<br />
| 415<br />
|<br />
| 12115<br />
| 12315<br />
| 12415<br />
|-<br />
|Dynamic Asymmetric Holland Model<ref group="note" name="nws19bad"></ref><br />
| 19<br />
| 119<br />
| 319<br />
| 419<br />
|<br />
| 12119<br />
| 12319<br />
| 12419<br />
|-<br />
|[[Generalized Asymmetric Holland Model]]<br />
| 20<br />
| 120<br />
| 320<br />
| 420<br />
|<br />
| 12120<br />
| 12320<br />
| 12420<br />
| <br />
|}<br />
<br />
== Notes ==<br />
<references group="note" /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1089IM2020-08-21T19:30:43Z<p>Wpringle: /* Default IM Values */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| -<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| -<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=+}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
The most recent version (55+) also has an option that improves the (default) semi-implicit consistent GWCE mass matrix mode to compute the complete (total depth) gravity wave term (free surface gradient) implicitly; toggled by setting IMDigit-6 to 3. The default version (IMDigit-6=1), only computes the initial still water depth component of the free surface gradient implicitly, which might make it more susceptible to CFL violations in shallow depths and can encounter Matrix diagonality issues overland where the initial still water depth is negative.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1088IM2020-07-19T22:02:47Z<p>Wpringle: /* Six-digit IM Codes */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111114<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611114<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=+}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
The most recent version (55+) also has an option that improves the (default) semi-implicit consistent GWCE mass matrix mode to compute the complete (total depth) gravity wave term (free surface gradient) implicitly; toggled by setting IMDigit-6 to 3. The default version (IMDigit-6=1), only computes the initial still water depth component of the free surface gradient implicitly, which might make it more susceptible to CFL violations in shallow depths and can encounter Matrix diagonality issues overland where the initial still water depth is negative.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1087IM2020-07-12T20:34:28Z<p>Wpringle: /* Default IM Values */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111114<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611114<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=+}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| A value of 4-6 does the same as 1-3 (same order) but in baroclinic mode<br/><code>CBaroclinic=.TRUE.</code> {{ADC version|version=55|relation=+}}<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| See above <br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
The most recent version (55+) also has an option that improves the (default) semi-implicit consistent GWCE mass matrix mode to compute the complete (total depth) gravity wave term (free surface gradient) implicitly; toggled by setting IMDigit-6 to 3. The default version (IMDigit-6=1), only computes the initial still water depth component of the free surface gradient implicitly, which might make it more susceptible to CFL violations in shallow depths and can encounter Matrix diagonality issues overland where the initial still water depth is negative.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1086IM2020-07-11T01:13:59Z<p>Wpringle: /* Six-digit IM Codes */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111113<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611113<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=+}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| A value of 4-6 does the same as 1-3 (same order) but in baroclinic mode<br/><code>CBaroclinic=.TRUE.</code> {{ADC version|version=55|relation=+}}<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| See above <br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
The most recent version (55+) also has an option that improves the (default) semi-implicit consistent GWCE mass matrix mode to compute the complete (total depth) gravity wave term (free surface gradient) implicitly; toggled by setting IMDigit-6 to 3. The default version (IMDigit-6=1), only computes the initial still water depth component of the free surface gradient implicitly, which might make it more susceptible to CFL violations in shallow depths and can encounter Matrix diagonality issues overland where the initial still water depth is negative.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1085AliDispersionControl2020-07-10T21:46:43Z<p>Wpringle: /* Version */</p>
<hr />
<div>'''Ali's Dispersion Correction''' is a correction to the dispersive behavior of very long waves in a compressible ocean on an elastic Earth. It is intended to be used in lieu of the self-attraction and loading tide prescribed through the [[fort.24 file]]. It should be used with the sixth-digit of [[IM]] equal to 3 (fully implicit gravity wave term).<br />
<br />
== Version ==<br />
{{Version support box|version=55|relation=+|support=tp}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing.<br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used (the following floats are the default values):<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn the Mach number based correction off. <br />
* <code>Ad</code> the power law constant. <br />
* <code>Bd</code> the power law exponent.<br />
<br />
For the following equation:<br />
<math>Correction = 1 - \frac{Ma^2}{4} - {A_d}H^{B_d}</math><br />
<br />
where Ma is the Mach number (<math>Ma = \frac{\sqrt{gH}}{Cs}</math>) and H is the total water depth.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1084IM2020-07-10T21:43:58Z<p>Wpringle: /* Six-digit IM Codes */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111113<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611113<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=+}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| A value of 4-6 does the same as 1-3 (same order) but in baroclinic mode<br/><code>CBaroclinic=.TRUE.</code> {{Version support box|version=55|relation=+|support=tp}}<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| See above <br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
The most recent version (55+) also has an option that improves the (default) semi-implicit consistent GWCE mass matrix mode to compute the complete (total depth) gravity wave term (free surface gradient) implicitly; toggled by setting IMDigit-6 to 3. The default version (IMDigit-6=1), only computes the initial still water depth component of the free surface gradient implicitly, which might make it more susceptible to CFL violations in shallow depths and can encounter Matrix diagonality issues overland where the initial still water depth is negative.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1069AliDispersionControl2020-07-10T15:26:45Z<p>Wpringle: </p>
<hr />
<div>'''Ali's Dispersion Correction''' is a correction to the dispersive behavior of very long waves in a compressible ocean on an elastic Earth. It is intended to be used in lieu of the self-attraction and loading tide prescribed through the [[fort.24 file]]. It should be used with the sixth-digit of [[IM]] equal to 3 (fully implicit gravity wave term).<br />
<br />
== Version ==<br />
{{ADC version|version=55: Technical Preview|relation=+}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing. <br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used (the following floats are the default values):<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn the Mach number based correction off. <br />
* <code>Ad</code> the power law constant. <br />
* <code>Bd</code> the power law exponent.<br />
<br />
For the following equation:<br />
<math>Correction = 1 - \frac{Ma^2}{4} - {A_d}H^{B_d}</math><br />
<br />
where Ma is the Mach number (<math>Ma = \frac{\sqrt{gH}}{Cs}</math>) and H is the total water depth.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1068AliDispersionControl2020-07-10T15:24:45Z<p>Wpringle: /* Controlling Dispersive Behavior */</p>
<hr />
<div>'''Ali's Dispersion Correction''' is a correction to the dispersive behavior of very long waves in a compressible ocean on an elastic Earth. It is intended to be used in lieu of the self-attraction and loading tide prescribed through the [[fort.24 file]].<br />
<br />
== Version ==<br />
{{ADC version|version=55: Technical Preview|relation=+}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing. <br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used (the following floats are the default values):<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn the Mach number based correction off. <br />
* <code>Ad</code> the power law constant. <br />
* <code>Bd</code> the power law exponent.<br />
<br />
For the following equation:<br />
<math>Correction = 1 - \frac{Ma^2}{4} - {A_d}H^{B_d}</math><br />
<br />
where Ma is the Mach number (<math>Ma = \frac{\sqrt{gH}}{Cs}</math>) and H is the total water depth.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1067AliDispersionControl2020-07-10T15:24:17Z<p>Wpringle: /* Controlling Dispersive Behavior */</p>
<hr />
<div>'''Ali's Dispersion Correction''' is a correction to the dispersive behavior of very long waves in a compressible ocean on an elastic Earth. It is intended to be used in lieu of the self-attraction and loading tide prescribed through the [[fort.24 file]].<br />
<br />
== Version ==<br />
{{ADC version|version=55: Technical Preview|relation=+}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing. <br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used (the following floats are the default values):<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn the Mach number based correction off. <br />
* <code>Ad</code> the power law constant. <br />
* <code>Bd</code> the power law exponent.<br />
<br />
For the following equation:<br />
<math>Correction = 1 - Ma^2/4 - {A_d}H^{B_d}</math><br />
<br />
where Ma is the Mach number (<math>Ma = \frac{\sqrt{gH}}{Cs}</math>) and H is the total water depth.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1066AliDispersionControl2020-07-10T15:23:59Z<p>Wpringle: /* Controlling Dispersive Behavior */</p>
<hr />
<div>'''Ali's Dispersion Correction''' is a correction to the dispersive behavior of very long waves in a compressible ocean on an elastic Earth. It is intended to be used in lieu of the self-attraction and loading tide prescribed through the [[fort.24 file]].<br />
<br />
== Version ==<br />
{{ADC version|version=55: Technical Preview|relation=+}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing. <br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used (the following floats are the default values):<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn the Mach number based correction off. <br />
* <code>Ad</code> the power law constant. <br />
* <code>Bd</code> the power law exponent.<br />
<br />
For the following equation:<br />
<math>Correction = 1 - Ma^2/4 - {A_d}H^{B_d}</math><br />
<br />
where Ma is the Mach number (<math>Ma = \frac{\sqrt(gH)}{Cs}</math>) and H is the total water depth.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1065AliDispersionControl2020-07-10T15:18:28Z<p>Wpringle: </p>
<hr />
<div>'''Ali's Dispersion Correction''' is a correction to the dispersive behavior of very long waves in a compressible ocean on an elastic Earth. It is intended to be used in lieu of the self-attraction and loading tide prescribed through the [[fort.24 file]].<br />
<br />
== Version ==<br />
{{ADC version|version=55: Technical Preview|relation=+}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing. <br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used (the following floats are the default values):<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn the Mach number based correction off. <br />
* <code>Ad</code> the power law constant. <br />
* <code>Bd</code> the power law exponent.<br />
<br />
For the following equations<br />
<math>Correction = A_d*h^{B_d}</math></div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1064AliDispersionControl2020-07-10T15:07:53Z<p>Wpringle: /* Controlling Dispersive Behavior */</p>
<hr />
<div>'''Dynamic Water Level Correction''' is a process by which modeled water levels are dynamically adjusted by use of a forcing term. The correction can be applied as constant or varying in space and/or time. The correction is applied as a forcing term in the momentum equations whose mathematical form is equivalent to that of an atmospheric pressure term. This means that, for gradually-varying corrections, corrected water levels should closely follow the input correction, though these may deviate if a correction is applied very quickly or to an area that has a very weak connection to an open boundary through which water can flow. Further discussion is below in the [[#FAQ|FAQ]]. <br />
<br />
Overviews and examples of this capability have been provided in multiple presentations (Luettich et al. 2017<ref name="luettichPres2017">Luettich, R.L., T.G. Asher, B.O. Blanton, J.G. Fleming. Representing Low Frequency, Spatially Varying Water Level Anomalies in Storm Surge Computations. 2017 American Meteorological Society Annual Meeting. [https://ams.confex.com/ams/97Annual/webprogram/Paper316033.html Link to talk]</ref>, Asher et al. 2018<ref name="asherPres2018">Asher, T.G., R.L. Luettich, J.G. Fleming, B.O.Blanton. Assimilation of Observed Water Levels into Storm Surge Model Predictions. 2018 American Meteorological Society Annual Meeting. [https://ams.confex.com/ams/98Annual/webprogram/Paper334044.html Link to talk]</ref>) and a journal article (Asher et al. 2019<ref name="asher2019">Asher, T.G., Luettich Jr., R.A., Fleming, J.G., Blanton, B.O., 2019. Low frequency water level correction in storm surge models using data assimilation. Ocean Modelling 144, 101483. https://doi.org/10.1016/j.ocemod.2019.101483</ref>) with details and an application to Hurricane Matthew. Users looking for ways to generate water level correction surfaces can look to that same article and this digital publication/data repository<ref>Asher, T., 2019. Hurricane Matthew (2016) Storm Surge and Wave Simulations with Data Assimilation. https://doi.org/10.17603/2Z8H-7K90</ref>, which holds the code base used in the aforementioned paper. <br />
<br />
== Version ==<br />
{{ADC version|version=55: Technical Preview|relation=+}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing. <br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used (the following floats are the default values):<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn the Mach number based correction off. <br />
* <code>Ad</code> the power law constant. <br />
* <code>Bd</code> the power law exponent.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1063AliDispersionControl2020-07-10T15:04:58Z<p>Wpringle: /* Controlling Dispersive Behavior */</p>
<hr />
<div>'''Dynamic Water Level Correction''' is a process by which modeled water levels are dynamically adjusted by use of a forcing term. The correction can be applied as constant or varying in space and/or time. The correction is applied as a forcing term in the momentum equations whose mathematical form is equivalent to that of an atmospheric pressure term. This means that, for gradually-varying corrections, corrected water levels should closely follow the input correction, though these may deviate if a correction is applied very quickly or to an area that has a very weak connection to an open boundary through which water can flow. Further discussion is below in the [[#FAQ|FAQ]]. <br />
<br />
Overviews and examples of this capability have been provided in multiple presentations (Luettich et al. 2017<ref name="luettichPres2017">Luettich, R.L., T.G. Asher, B.O. Blanton, J.G. Fleming. Representing Low Frequency, Spatially Varying Water Level Anomalies in Storm Surge Computations. 2017 American Meteorological Society Annual Meeting. [https://ams.confex.com/ams/97Annual/webprogram/Paper316033.html Link to talk]</ref>, Asher et al. 2018<ref name="asherPres2018">Asher, T.G., R.L. Luettich, J.G. Fleming, B.O.Blanton. Assimilation of Observed Water Levels into Storm Surge Model Predictions. 2018 American Meteorological Society Annual Meeting. [https://ams.confex.com/ams/98Annual/webprogram/Paper334044.html Link to talk]</ref>) and a journal article (Asher et al. 2019<ref name="asher2019">Asher, T.G., Luettich Jr., R.A., Fleming, J.G., Blanton, B.O., 2019. Low frequency water level correction in storm surge models using data assimilation. Ocean Modelling 144, 101483. https://doi.org/10.1016/j.ocemod.2019.101483</ref>) with details and an application to Hurricane Matthew. Users looking for ways to generate water level correction surfaces can look to that same article and this digital publication/data repository<ref>Asher, T., 2019. Hurricane Matthew (2016) Storm Surge and Wave Simulations with Data Assimilation. https://doi.org/10.17603/2Z8H-7K90</ref>, which holds the code base used in the aforementioned paper. <br />
<br />
== Version ==<br />
{{ADC version|version=55: Technical Preview|relation=+}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing. <br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used:<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn the Mach number based correction off. <br />
* <code>Ad</code> the power law constant. <br />
* <code>Bd</code> the power law exponent.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1062AliDispersionControl2020-07-10T15:04:02Z<p>Wpringle: /* Controlling Dispersive Behavior */</p>
<hr />
<div>'''Dynamic Water Level Correction''' is a process by which modeled water levels are dynamically adjusted by use of a forcing term. The correction can be applied as constant or varying in space and/or time. The correction is applied as a forcing term in the momentum equations whose mathematical form is equivalent to that of an atmospheric pressure term. This means that, for gradually-varying corrections, corrected water levels should closely follow the input correction, though these may deviate if a correction is applied very quickly or to an area that has a very weak connection to an open boundary through which water can flow. Further discussion is below in the [[#FAQ|FAQ]]. <br />
<br />
Overviews and examples of this capability have been provided in multiple presentations (Luettich et al. 2017<ref name="luettichPres2017">Luettich, R.L., T.G. Asher, B.O. Blanton, J.G. Fleming. Representing Low Frequency, Spatially Varying Water Level Anomalies in Storm Surge Computations. 2017 American Meteorological Society Annual Meeting. [https://ams.confex.com/ams/97Annual/webprogram/Paper316033.html Link to talk]</ref>, Asher et al. 2018<ref name="asherPres2018">Asher, T.G., R.L. Luettich, J.G. Fleming, B.O.Blanton. Assimilation of Observed Water Levels into Storm Surge Model Predictions. 2018 American Meteorological Society Annual Meeting. [https://ams.confex.com/ams/98Annual/webprogram/Paper334044.html Link to talk]</ref>) and a journal article (Asher et al. 2019<ref name="asher2019">Asher, T.G., Luettich Jr., R.A., Fleming, J.G., Blanton, B.O., 2019. Low frequency water level correction in storm surge models using data assimilation. Ocean Modelling 144, 101483. https://doi.org/10.1016/j.ocemod.2019.101483</ref>) with details and an application to Hurricane Matthew. Users looking for ways to generate water level correction surfaces can look to that same article and this digital publication/data repository<ref>Asher, T., 2019. Hurricane Matthew (2016) Storm Surge and Wave Simulations with Data Assimilation. https://doi.org/10.17603/2Z8H-7K90</ref>, which holds the code base used in the aforementioned paper. <br />
<br />
== Version ==<br />
{{ADC version|version=55: Technical Preview|relation=+}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing. <br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used:<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn this component off. <br />
* <code>Ad</code> the power law constant. <br />
* <code>Bd</code> the power law exponent.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=AliDispersionControl&diff=1061AliDispersionControl2020-07-10T15:02:04Z<p>Wpringle: Created page with "'''Dynamic Water Level Correction''' is a process by which modeled water levels are dynamically adjusted by use of a forcing term. The correction can be applied as constant o..."</p>
<hr />
<div>'''Dynamic Water Level Correction''' is a process by which modeled water levels are dynamically adjusted by use of a forcing term. The correction can be applied as constant or varying in space and/or time. The correction is applied as a forcing term in the momentum equations whose mathematical form is equivalent to that of an atmospheric pressure term. This means that, for gradually-varying corrections, corrected water levels should closely follow the input correction, though these may deviate if a correction is applied very quickly or to an area that has a very weak connection to an open boundary through which water can flow. Further discussion is below in the [[#FAQ|FAQ]]. <br />
<br />
Overviews and examples of this capability have been provided in multiple presentations (Luettich et al. 2017<ref name="luettichPres2017">Luettich, R.L., T.G. Asher, B.O. Blanton, J.G. Fleming. Representing Low Frequency, Spatially Varying Water Level Anomalies in Storm Surge Computations. 2017 American Meteorological Society Annual Meeting. [https://ams.confex.com/ams/97Annual/webprogram/Paper316033.html Link to talk]</ref>, Asher et al. 2018<ref name="asherPres2018">Asher, T.G., R.L. Luettich, J.G. Fleming, B.O.Blanton. Assimilation of Observed Water Levels into Storm Surge Model Predictions. 2018 American Meteorological Society Annual Meeting. [https://ams.confex.com/ams/98Annual/webprogram/Paper334044.html Link to talk]</ref>) and a journal article (Asher et al. 2019<ref name="asher2019">Asher, T.G., Luettich Jr., R.A., Fleming, J.G., Blanton, B.O., 2019. Low frequency water level correction in storm surge models using data assimilation. Ocean Modelling 144, 101483. https://doi.org/10.1016/j.ocemod.2019.101483</ref>) with details and an application to Hurricane Matthew. Users looking for ways to generate water level correction surfaces can look to that same article and this digital publication/data repository<ref>Asher, T., 2019. Hurricane Matthew (2016) Storm Surge and Wave Simulations with Data Assimilation. https://doi.org/10.17603/2Z8H-7K90</ref>, which holds the code base used in the aforementioned paper. <br />
<br />
== Version ==<br />
{{ADC version|version=55: Technical Preview|relation=+}} <br />
This is considered a technical preview in version 55. Theoretical work is still ongoing. <br />
<br />
==Controlling Dispersive Behavior==<br />
The feature is triggered by the presence of the &AliDispersionControl namelist at the bottom of the [[fort.15 file]]. Here is an example of how this line is used:<br/><br />
<code>&AliDispersionControl CAliDisp=T, Cs=1500.0, Ad = 0.0050189, Bd = 0.23394/</code><br />
<br />
* <code>CAliDisp</code> logical flag to turn Ali's dispersion correction on (F=false by default). <br />
* <code>Cs</code> is the speed of sound in water [m/s] for the Mach number based dispersive correction. Set Cs to a negative value to turn this component off. <br />
* <code>Ad</code> coefficient. <br />
* <code>Bd</code> coefficient.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=Fort.15_file_format&diff=1058Fort.15 file format2020-07-10T14:32:49Z<p>Wpringle: /* Namelists */</p>
<hr />
<div>The basic file structure of the [[fort.15 file]] is shown below. Each line of input data is represented by a line containing the input variable name(s). Inputs in the [[fort.15 file]] must be entered in the exact order shown on this page. '''Blank lines and headings are only to enhance readability.''' Loops indicate multiple lines of input. <br/><br />
==Main Controls==<br />
===Metadata and Logging===<br />
<!-- Comments<br />
<code>[[RUNDES]]</code><br/><br />
<code>[[RUNID]]</code><br/><br />
<code>[[NFOVER]]</code><br/><br />
<code>[[NABOUT]]</code><br/><br />
<code>[[NSCREEN]]</code><br/><br />
--><br />
<br />
{| class="wikitable" border="1" style="text-align: center"<br />
|-<br />
! Parameter<br />
! Type<br />
! Required?<br />
! Description<br />
! Values<br />
|-<br />
| <code>RUNDES</code><br />
| <math>\leq</math>32 character string<br />
| Always<br />
| Run description<br />
| Any alpha-numeric <br />
|- style="background:#efefef;"<br />
| <code>RUNID</code><br />
| <math>\leq</math>24 character string<br />
| Always<br />
| Run identification<br />
| Any alpha-numeric <br />
|- style="background:#efefef;"<br />
| <code>[[NFOVER]]</code><br />
| integer<br />
| Always<br />
| Non-fatal error override option <br />
| 0 or 1<br />
|- style="background:#efefef;"<br />
| <code>[[NABOUT]]</code><br />
| integer<br />
| Always<br />
| Logging level <br />
| -1, 0, 1, 2, or 3<br />
|- style="background:#efefef;"<br />
| <code>[[NSCREEN]]</code><br />
| integer<br />
| Always<br />
| Logging output destination<br />
| -1, 0, or 1<br />
|}<br />
<br />
===Numerics & Physics===<br />
'''<code>[[IHOT]]'''</code> - whether to read a hotstart file<br/><br />
'''<code>[[ICS]]</code>''' - coordinate projection to run in<br/><br />
'''<code>[[IM]]</code>''' - model run mode<br/><br />
'''<code>[[IDEN]]</code>''' - density forcing mode, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[IM]] = 20, 21, 30, 31</code>, or if <code>IM</code>'s last 6-digit entry is > 4 (e.g., <code>51311<b>5</b></code>).<br/><br />
'''<code>[[NOLIBF]]</code>''' - bottom stress parameterization mode<br/><br />
'''<code>[[NOLIFA]]</code>''' - finite amplitude term mode<br/><br />
'''<code>[[NOLICA]]</code>''' - advection term mode<br/><br />
'''<code>[[NOLICAT]]</code>''' - advection term mode<br/><br />
'''<code>[[NWP]]</code>''' - number of [[nodal attribute]]s<br/><br />
''for j=1 to <code>[[NWP]]</code>''<br/><br />
: '''<code>[[AttrName(j)]]</code>''' - nodal attributes to use, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NWP]] > 0</code><br/><br />
''end j loop''<br/><br />
'''<code>[[NCOR]]</code>''' - Coriolis control<br/><br />
'''<code>[[NTIP]]</code>''' - tidal potential forcing control<br/><br />
<span id="NWS"/>'''<code>[[NWS#Value_Seen_in_fort.15_File|NWS]]</code>''' - meteorological, wave, and ice forcing control<br/><br />
'''<code>[[NRAMP]]</code>''' - forcing ramping control<br/><br />
'''<code>[[G]]</code>''' - acceleration due to gravity<br/><br />
'''<code>[[TAU0]]</code>''' - affects numerical diffusion/stability of governing equations<br/><br />
'''<code>[[Tau0FullDomainMin]] [[Tau0FullDomainMax]]</code>''' - limits on <code>[[TAU0]]</code>, ''<span style="background:blanchedalmond">include this line only if:</span>''<code>TAU0 = -5.0</code>.<br/><br />
'''<code>[[DTDP]]</code>''' - model time step (seconds) and predictor-corrector control<br/><br />
'''<code>[[STATIM]]</code>''' - shifts numeric value of starting simulation time (days)<br/><br />
<span id="REFTIM"/>'''<code>[[REFTIM]]</code>''' - shifts reference time (days) for tidal harmonic analysis<br/><br />
<span id="WTIMINC"/><span id="RSTIMINC"/><span id="CICE_TIMINC"/>'''[[Supplemental_Meteorological/Wave/Ice_Parameters|Meteorological controls including <code>WTIMINC, RSTIMINC</code>]]''' - ''<span style="background:blanchedalmond">include this line:</span>'' for most cases of <code>NWS ≠ 0</code>, see linked page for details.<br/><br />
'''<code>[[RNDAY]]</code>''' - end time of simulation (days)<br/><br />
'''[[Ramping|Ramping controls including <code>DRAMP, FluxSettlingTime</code>]]''' - ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>NRAMP > 0</code>, see linked page for details<br/><br />
'''<code>[[A00 B00 C00]]</code>''' - time weighting factors in GWCE<br/><br />
'''<code>[[H0]]</code>''' - minimum depth, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIFA]] = 0</code> or <code>1</code>''<br/><br />
'''<code>[[H0]]</code> <code>INTEGER</code> <code>INTEGER</code> <code>[[VELMIN]]</code>''' - alternate minimum depth controls, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIFA]] = 2</code> or <code>3</code>''<br/><br />
'''<code>[[SLAM0]] [[SFEA0]]</code>''' - longitude and latitude for center of CPP coordinate projection<br/><br />
'''<code>[[TAU]]</code>''' - linear bottom friction coefficient, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIBF]] = 0</code>''<br/><br />
'''<code>[[CF]]</code>''' - quadratic bottom friction coefficient or limit, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIBF]] = 1</code>''<br/><br />
'''<code>[[CF]]</code> <code>[[HBREAK]]</code> <code>[[FTHETA]]</code> <code>[[FGAMMA]]</code>''' - alternate quadratic bottom friction controls, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIBF]] = 2</code>''<br/><br />
'''<code>[[ESLM]]</code>''' - horizontal eddy viscosity controls, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[IM]] = 0</code>, <code>1</code>, or <code>2</code>''<br/><br />
'''<code>[[ESLM]]</code> <code>[[ESLC]]</code>''' - alternate horizontal eddy viscosity controls, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[IM]] = 10</code>''<br/><br />
'''<code>[[CORI]]</code>''' - constant Coriolis coefficient, ''<span style="background:mistyRose">only used if</span>'' <code>NCOR=0</code><br />
<br />
===Periodic (Tidal) Body Forcing===<br />
[[NTIF]]<br/><br />
for k=1 to [[NTIF]]<br />
: [[TIPOTAG(k)]]<br />
: [[TPK(k)]], [[AMIGT(k)]], [[ETRF(k)]], [[FFT(k)]], [[FACET(k)]]<br />
end k loop<br />
<br />
===Periodic (Tidal) Boundary Elevations===<br />
[[NBFR]]<br/><br />
for k=1 to [[NBFR]]<br/><br />
: [[BOUNTAG(k)]]<br/><br />
: [[AMIG(k)]], [[FF(k)]], [[FACE(k)]]<br/><br />
end k loop<br/><br />
for k=1 to [[NBFR]]<br/><br />
: [[ALPHA(k)]]<br/><br />
: for j=1 to [[NETA]]<br/><br />
:: [[EMO(k,j), EFA(k,j)]]<br/><br />
: end j loop<br/><br />
end k loop<br />
<br />
===Periodic (Tidal) Boundary Velocities===<br />
[[ANGINN]]<br/><br />
[[NFFR]] - include this line only if [[IBTYPE]] = 2, 12, 22, 32 or 52 in the Grid and Boundary Information File<br/><br />
for k=1 to [[NFFR]]<br/><br />
: [[FBOUNTAG(k)]]<br/><br />
: [[FAMIGT(k),FFF(k),FFACE(k)]]<br/><br />
end k loop<br/><br />
for k=1 to [[NFFR]]<br/><br />
: [[ALPHA(k)]]<br/><br />
: for j=1 to [[NVEL]]<br/><br />
:: [[QNAM(k,j), QNPH(k,j)]] - use this line if [[IBTYPE]] = 2, 12, 22 in the Grid and Boundary Information File<br/><br />
:: [[QNAM(k,j), QNPH(k,j)]], [[ENAM(k,j), ENPH(k,j)]] - use this line if [[IBTYPE]] = 32 in the Grid and Boundary Information File<br/><br />
:end j loop<br/><br />
end k loop<br />
<br />
===Model Output===<br />
[[NOUTE]], [[TOUTSE]], [[TOUTFE]], [[NSPOOLE]]<br/><br />
[[NSTAE]]<br/><br />
for k=1 to [[NSTAE]]<br/><br />
: [[XEL(k), YEL(k)]] - use these lines if [[NSTAE]] is positive. If negative, stations are listed in the [[elev_stat.151]] file<br/><br />
end k loop<br/><br />
[[NOUTV]], [[TOUTSV]], [[TOUTFV]], [[NSPOOLV]]<br/><br />
[[NSTAV]]<br/><br />
for k=1 to [[NSTAV]]<br/><br />
: [[XEV(k), YEV(k)]] - use these lines if [[NSTAV]] is positive. If negative, stations are listed in the [[vel_stat.151]] file<br/><br />
end k loop<br/><br />
[[NOUTC]], [[TOUTSC]], [[TOUTFC]], [[NSPOOLC]] - include this line only if IM =10<br/><br />
[[NSTAC]] - include this line only if IM =10<br/><br />
for k=1 to [[NSTAC]]<br/><br />
: [[XEC(k), YEC(k)]]<br/><br />
end k loop<br/><br />
[[NOUTM]], [[TOUTSM]], [[TOUTFM]], [[NSPOOLM]] - include this line only if NWS is not equal to zero.<br/><br />
[[NSTAM]] - include this line only if NWS is not equal to zero.<br/><br />
for k=1 to [[NSTAM]]<br/><br />
: [[XEM(k), YEM(k)]] - use these lines if [[NSTAM]] is positive. If negative, stations are listed in the [[met_stat.151]] file<br/><br />
end k loop<br/><br />
[[NOUTGE]], [[TOUTSGE]], [[TOUTFGE]], [[NSPOOLGE]]<br/><br />
[[NOUTGV]], [[TOUTSGV]], [[TOUTFGV]], [[NSPOOLGV]]<br/><br />
[[NOUTGC]], [[TOUTSGC]], [[TOUTFGC]], [[NSPOOLGC]] - include this line only if IM =10<br/><br />
[[NOUTGW]], [[TOUTSGW]], [[TOUTFGW]], [[NSPOOLGW]] - include this line only if NWS is not equal to zero.<br />
<br />
====Harmonic Analysis====<br />
[[NFREQ]]<br/><br />
for k=1 to [[NFREQ]]<br/><br />
: [[NAMEFR(k)]]<br/><br />
: [[HAFREQ(k), HAFF(k), HAFACE(k)]]<br/><br />
end k loop<br/><br />
[[THAS]], [[THAF]], [[NHAINC]], [[FMV]]<br/><br />
[[NHASE]], [[NHASV]], [[NHAGE]], [[NHAGV]]<br />
<br />
====Hotstart Output and Numeric Controls====<br />
[[NHSTAR]], [[NHSINC]]<br/><br />
[[ITITER]], [[ISLDIA]], [[CONVCR]], [[ITMAX]]<br/><br />
<br />
''For a 2DDI ADCIRC run that does not use netCDF nor namelists, the file ends here. For those controls, see further below in the [[#NetCDF Controls|NetCDF Controls]] and [[#Namelists|Namelists]] sections.''<br />
<br />
==3D Model Run==<br />
[[IDEN]]<br/><br />
[[ISLIP]], [[KP]]<br/><br />
[[Z0S,Z0B]]<br/><br />
[[ALP1,ALP2,ALP3]]<br/><br />
[[IGC]], [[NFEN]]<br/><br />
for k=1 to [[NFEN]] (include this loop only if [[IGC]] = 0, k=1 at bottom, k= [[NFEN]] at surface)<br/><br />
: [[SIGMA(k)]]<br/><br />
end k loop<br/><br />
[[IEVC]], [[EVMIN]], [[EVCON]]<br/><br />
for k=1 to [[NFEN]] (include this loop only if [[IEVC]] = 0, k=1 at bottom, k= [[NFEN]] at surface)<br/><br />
: [[EVTOT(k)]]<br/><br />
end k loop<br/><br />
[[THETA1, THETA2]](include this line only if [[IEVC]] = 50 or 51)<br/><br />
[[I3DSD,TO3DSDS,TO3DSDF,NSPO3DSD]]<br/><br />
[[NSTA3DD]]<br/><br />
for k=1 to [[NSTA3DD]]<br/><br />
: [[X3DS(k), Y3DS(k)]]<br/><br />
end k loop<br/><br />
[[I3DSV,TO3DSVS,TO3DSVF,NSPO3DSV]]<br/><br />
[[NSTA3DV]]<br/><br />
for k=1 to [[NSTA3DV]]<br/><br />
: [[X3DS(k), Y3DS(k)]]<br/><br />
end k loop<br/><br />
[[I3DST,TO3DSTS,TO3DSTF,NSPO3DST]]<br/><br />
[[NSTA3DT]]<br/><br />
for k=1 to [[NSTA3DT]]<br/><br />
: [[X3DS(k), Y3DS(k)]]<br/><br />
end k loop<br/><br />
[[I3DGD]],[[TO3DGDS]],[[TO3DGDF]],[[NSPO3DGD]]<br/><br />
[[I3DGV]],[[TO3DGVS]],[[TO3DGVF]],[[NSPO3DGV]]<br/><br />
[[I3DGT]],[[TO3DGTS]],[[TO3DGTF]],[[NSPO3DGT]]<br/><br />
The following line will be read in if [[IM]] is 21 or 31.<br/><br />
[[RES_BC_FLAG]], [[BCFLAG_LNM]], [[BCFLAG_TEMP]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] is negative.<br/><br />
[[RBCTIMEINC]]<br/><br />
[[BCSTATIM]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] = 2.<br/><br />
[[RBCTIMEINC]], [[SBCTIMEINC]]<br/><br />
[[BCSTATIM]], [[SBCSTATIM]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] = 3.<br/><br />
[[RBCTIMEINC]], [[TBCTIMEINC]]<br/><br />
[[BCSTATIM]], [[TBCSTATIM]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] = 4.<br/><br />
[[RBCTIMEINC]], [[SBCTIMEINC]], [[TBCTIMEINC]]<br/><br />
[[BCSTATIM]], [[SBCSTATIM]], [[TBCSTATIM]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] = 3 or 4 and [[BCFLAG_TEMP]] is not equal to 0.<br/><br />
[[TTBCTIMEINC]], [[TTBCSTATIM]]<br/><br />
[[TTBCTIMEINC]]<br/><br />
The following two lines will be read in only if [[IM]] is 21 or 31.<br/><br />
[[SPONGEDIST]]<br/><br />
[[EQNSTATE]]<br/><br />
The following lines will be read in only if [[IDEN]] is > 0.<br/><br />
[[NLSD, NVSD]]<br/><br />
[[NLTD, NVTD]]<br/><br />
[[ALP4]]<br/><br />
The following line will be read in only if [[IDEN]] = 3 or 4.<br/><br />
[[NTF]]<br />
<br />
==NetCDF Controls==<br />
The following lines will be read in only if the NetCDF output or hotstart format is chosen<br/><br />
NCPROJ<br/><br />
NCINST<br/><br />
NCSOUR<br/><br />
NCHIST<br/><br />
NCREF<br/><br />
NCCOM<br/><br />
NCHOST<br/><br />
NCCONV<br/><br />
NCCONT<br/><br />
NCDATE<br />
<br />
==Namelists==<br />
The following Fortran namelist lines are optional, but if they appear, they must appear at the very end of the fort.15 file.<br/><br />
<code>&metControl WindDragLimit=floatValue, DragLawString='stringValue', rhoAir=floatValue, outputWindDrag=logicalValue /</code><br/><br />
<code>&timeBathyControl NDDT=integerValue, BTIMINC=floatValue, BCHGTIMINC=floatValue /</code><br/><br />
<code>&waveCoupling WindWaveMultiplier=floatValue /</code><br/><br />
<code>&SWANOutputControl SWAN_OutputHS=logicalValue, SWAN_OutputDIR=logicalValue, SWAN_OutputTM01=logicalValue, SWAN_OutputTPS=logicalValue, SWAN_OutputWIND=logicalValue, SWAN_OutputTM02=logicalValue, SWAN_OutputTMM10=logicalValue /</code><br/><br />
<code>&subdomainModeling subdomainOn=logicalValue/</code><br/><br />
<code>&wetDryControl outputNodeCode=logicalValue, outputNOFF=logicalValue, noffActive=logicalValue /</code><br/><br />
<code>&inundationOutputControl inundationOutput=logicalValue0, inunThresh =floatValue /</code><br/><br />
<code>&TVWControl use_TVW=logicalValue, TVW_file='stringValue', nout_TVW =integerValue, touts_TVW =floatValue, toutf_TVW=floatValue, nspool_TVW =integerValue /</code><br/><br />
<code>&WarnElevControl WarnElev=floatValue, ErrorElev=floatValue, WarnElevDump=logicalValue, WarnElevDumpLimit=integerValue /</code><br/><br />
<code>[[Dynamic_water_level_correction#Controlling_Water_Level_Correction|&dynamicWaterLevelCorrectionControl]] dynamicWaterLevelCorrectionFileName='stringValue' dynamicWaterLevelCorrectionMultiplier=floatValue, dynamicWaterLevelCorrectionRampStart=floatValue, dynamicWaterLevelCorrectionRampEnd=floatValue, dynamicWaterLevelCorrectionRampReferenceTime='stringValue', dynamicWaterLevelCorrectionSkipSnaps=integerValue /</code><br/><br />
<code>&[[AliDispersionControl]] CAliDisp=logicalValue, Cs=floatValue, Ad=floatValue, Bd=floatValue /</code><br/><br />
[[category:input files]]</div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1057IM2020-07-10T14:31:50Z<p>Wpringle: /* Six-digit IM Codes */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111113<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611113<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=+}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| A value of 4-6 does the same as 1-3 (same order) but in baroclinic mode<br/><code>CBaroclinic=.TRUE.</code> {{ADC version|version=55: Technical Preview|relation=+}}<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| See above <br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
The most recent version (55+) also has an option that improves the (default) semi-implicit consistent GWCE mass matrix mode to compute the complete (total depth) gravity wave term (free surface gradient) implicitly; toggled by setting IMDigit-6 to 3. The default version (IMDigit-6=1), only computes the initial still water depth component of the free surface gradient implicitly, which might make it more susceptible to CFL violations in shallow depths and can encounter Matrix diagonality issues overland where the initial still water depth is negative.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Fort.15_file_format&diff=1056Fort.15 file format2020-07-10T14:26:44Z<p>Wpringle: /* Namelists */</p>
<hr />
<div>The basic file structure of the [[fort.15 file]] is shown below. Each line of input data is represented by a line containing the input variable name(s). Inputs in the [[fort.15 file]] must be entered in the exact order shown on this page. '''Blank lines and headings are only to enhance readability.''' Loops indicate multiple lines of input. <br/><br />
==Main Controls==<br />
===Metadata and Logging===<br />
<!-- Comments<br />
<code>[[RUNDES]]</code><br/><br />
<code>[[RUNID]]</code><br/><br />
<code>[[NFOVER]]</code><br/><br />
<code>[[NABOUT]]</code><br/><br />
<code>[[NSCREEN]]</code><br/><br />
--><br />
<br />
{| class="wikitable" border="1" style="text-align: center"<br />
|-<br />
! Parameter<br />
! Type<br />
! Required?<br />
! Description<br />
! Values<br />
|-<br />
| <code>RUNDES</code><br />
| <math>\leq</math>32 character string<br />
| Always<br />
| Run description<br />
| Any alpha-numeric <br />
|- style="background:#efefef;"<br />
| <code>RUNID</code><br />
| <math>\leq</math>24 character string<br />
| Always<br />
| Run identification<br />
| Any alpha-numeric <br />
|- style="background:#efefef;"<br />
| <code>[[NFOVER]]</code><br />
| integer<br />
| Always<br />
| Non-fatal error override option <br />
| 0 or 1<br />
|- style="background:#efefef;"<br />
| <code>[[NABOUT]]</code><br />
| integer<br />
| Always<br />
| Logging level <br />
| -1, 0, 1, 2, or 3<br />
|- style="background:#efefef;"<br />
| <code>[[NSCREEN]]</code><br />
| integer<br />
| Always<br />
| Logging output destination<br />
| -1, 0, or 1<br />
|}<br />
<br />
===Numerics & Physics===<br />
'''<code>[[IHOT]]'''</code> - whether to read a hotstart file<br/><br />
'''<code>[[ICS]]</code>''' - coordinate projection to run in<br/><br />
'''<code>[[IM]]</code>''' - model run mode<br/><br />
'''<code>[[IDEN]]</code>''' - density forcing mode, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[IM]] = 20, 21, 30, 31</code>, or if <code>IM</code>'s last 6-digit entry is > 4 (e.g., <code>51311<b>5</b></code>).<br/><br />
'''<code>[[NOLIBF]]</code>''' - bottom stress parameterization mode<br/><br />
'''<code>[[NOLIFA]]</code>''' - finite amplitude term mode<br/><br />
'''<code>[[NOLICA]]</code>''' - advection term mode<br/><br />
'''<code>[[NOLICAT]]</code>''' - advection term mode<br/><br />
'''<code>[[NWP]]</code>''' - number of [[nodal attribute]]s<br/><br />
''for j=1 to <code>[[NWP]]</code>''<br/><br />
: '''<code>[[AttrName(j)]]</code>''' - nodal attributes to use, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NWP]] > 0</code><br/><br />
''end j loop''<br/><br />
'''<code>[[NCOR]]</code>''' - Coriolis control<br/><br />
'''<code>[[NTIP]]</code>''' - tidal potential forcing control<br/><br />
<span id="NWS"/>'''<code>[[NWS#Value_Seen_in_fort.15_File|NWS]]</code>''' - meteorological, wave, and ice forcing control<br/><br />
'''<code>[[NRAMP]]</code>''' - forcing ramping control<br/><br />
'''<code>[[G]]</code>''' - acceleration due to gravity<br/><br />
'''<code>[[TAU0]]</code>''' - affects numerical diffusion/stability of governing equations<br/><br />
'''<code>[[Tau0FullDomainMin]] [[Tau0FullDomainMax]]</code>''' - limits on <code>[[TAU0]]</code>, ''<span style="background:blanchedalmond">include this line only if:</span>''<code>TAU0 = -5.0</code>.<br/><br />
'''<code>[[DTDP]]</code>''' - model time step (seconds) and predictor-corrector control<br/><br />
'''<code>[[STATIM]]</code>''' - shifts numeric value of starting simulation time (days)<br/><br />
<span id="REFTIM"/>'''<code>[[REFTIM]]</code>''' - shifts reference time (days) for tidal harmonic analysis<br/><br />
<span id="WTIMINC"/><span id="RSTIMINC"/><span id="CICE_TIMINC"/>'''[[Supplemental_Meteorological/Wave/Ice_Parameters|Meteorological controls including <code>WTIMINC, RSTIMINC</code>]]''' - ''<span style="background:blanchedalmond">include this line:</span>'' for most cases of <code>NWS ≠ 0</code>, see linked page for details.<br/><br />
'''<code>[[RNDAY]]</code>''' - end time of simulation (days)<br/><br />
'''[[Ramping|Ramping controls including <code>DRAMP, FluxSettlingTime</code>]]''' - ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>NRAMP > 0</code>, see linked page for details<br/><br />
'''<code>[[A00 B00 C00]]</code>''' - time weighting factors in GWCE<br/><br />
'''<code>[[H0]]</code>''' - minimum depth, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIFA]] = 0</code> or <code>1</code>''<br/><br />
'''<code>[[H0]]</code> <code>INTEGER</code> <code>INTEGER</code> <code>[[VELMIN]]</code>''' - alternate minimum depth controls, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIFA]] = 2</code> or <code>3</code>''<br/><br />
'''<code>[[SLAM0]] [[SFEA0]]</code>''' - longitude and latitude for center of CPP coordinate projection<br/><br />
'''<code>[[TAU]]</code>''' - linear bottom friction coefficient, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIBF]] = 0</code>''<br/><br />
'''<code>[[CF]]</code>''' - quadratic bottom friction coefficient or limit, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIBF]] = 1</code>''<br/><br />
'''<code>[[CF]]</code> <code>[[HBREAK]]</code> <code>[[FTHETA]]</code> <code>[[FGAMMA]]</code>''' - alternate quadratic bottom friction controls, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[NOLIBF]] = 2</code>''<br/><br />
'''<code>[[ESLM]]</code>''' - horizontal eddy viscosity controls, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[IM]] = 0</code>, <code>1</code>, or <code>2</code>''<br/><br />
'''<code>[[ESLM]]</code> <code>[[ESLC]]</code>''' - alternate horizontal eddy viscosity controls, ''<span style="background:blanchedalmond">include this line only if:</span>'' <code>[[IM]] = 10</code>''<br/><br />
'''<code>[[CORI]]</code>''' - constant Coriolis coefficient, ''<span style="background:mistyRose">only used if</span>'' <code>NCOR=0</code><br />
<br />
===Periodic (Tidal) Body Forcing===<br />
[[NTIF]]<br/><br />
for k=1 to [[NTIF]]<br />
: [[TIPOTAG(k)]]<br />
: [[TPK(k)]], [[AMIGT(k)]], [[ETRF(k)]], [[FFT(k)]], [[FACET(k)]]<br />
end k loop<br />
<br />
===Periodic (Tidal) Boundary Elevations===<br />
[[NBFR]]<br/><br />
for k=1 to [[NBFR]]<br/><br />
: [[BOUNTAG(k)]]<br/><br />
: [[AMIG(k)]], [[FF(k)]], [[FACE(k)]]<br/><br />
end k loop<br/><br />
for k=1 to [[NBFR]]<br/><br />
: [[ALPHA(k)]]<br/><br />
: for j=1 to [[NETA]]<br/><br />
:: [[EMO(k,j), EFA(k,j)]]<br/><br />
: end j loop<br/><br />
end k loop<br />
<br />
===Periodic (Tidal) Boundary Velocities===<br />
[[ANGINN]]<br/><br />
[[NFFR]] - include this line only if [[IBTYPE]] = 2, 12, 22, 32 or 52 in the Grid and Boundary Information File<br/><br />
for k=1 to [[NFFR]]<br/><br />
: [[FBOUNTAG(k)]]<br/><br />
: [[FAMIGT(k),FFF(k),FFACE(k)]]<br/><br />
end k loop<br/><br />
for k=1 to [[NFFR]]<br/><br />
: [[ALPHA(k)]]<br/><br />
: for j=1 to [[NVEL]]<br/><br />
:: [[QNAM(k,j), QNPH(k,j)]] - use this line if [[IBTYPE]] = 2, 12, 22 in the Grid and Boundary Information File<br/><br />
:: [[QNAM(k,j), QNPH(k,j)]], [[ENAM(k,j), ENPH(k,j)]] - use this line if [[IBTYPE]] = 32 in the Grid and Boundary Information File<br/><br />
:end j loop<br/><br />
end k loop<br />
<br />
===Model Output===<br />
[[NOUTE]], [[TOUTSE]], [[TOUTFE]], [[NSPOOLE]]<br/><br />
[[NSTAE]]<br/><br />
for k=1 to [[NSTAE]]<br/><br />
: [[XEL(k), YEL(k)]] - use these lines if [[NSTAE]] is positive. If negative, stations are listed in the [[elev_stat.151]] file<br/><br />
end k loop<br/><br />
[[NOUTV]], [[TOUTSV]], [[TOUTFV]], [[NSPOOLV]]<br/><br />
[[NSTAV]]<br/><br />
for k=1 to [[NSTAV]]<br/><br />
: [[XEV(k), YEV(k)]] - use these lines if [[NSTAV]] is positive. If negative, stations are listed in the [[vel_stat.151]] file<br/><br />
end k loop<br/><br />
[[NOUTC]], [[TOUTSC]], [[TOUTFC]], [[NSPOOLC]] - include this line only if IM =10<br/><br />
[[NSTAC]] - include this line only if IM =10<br/><br />
for k=1 to [[NSTAC]]<br/><br />
: [[XEC(k), YEC(k)]]<br/><br />
end k loop<br/><br />
[[NOUTM]], [[TOUTSM]], [[TOUTFM]], [[NSPOOLM]] - include this line only if NWS is not equal to zero.<br/><br />
[[NSTAM]] - include this line only if NWS is not equal to zero.<br/><br />
for k=1 to [[NSTAM]]<br/><br />
: [[XEM(k), YEM(k)]] - use these lines if [[NSTAM]] is positive. If negative, stations are listed in the [[met_stat.151]] file<br/><br />
end k loop<br/><br />
[[NOUTGE]], [[TOUTSGE]], [[TOUTFGE]], [[NSPOOLGE]]<br/><br />
[[NOUTGV]], [[TOUTSGV]], [[TOUTFGV]], [[NSPOOLGV]]<br/><br />
[[NOUTGC]], [[TOUTSGC]], [[TOUTFGC]], [[NSPOOLGC]] - include this line only if IM =10<br/><br />
[[NOUTGW]], [[TOUTSGW]], [[TOUTFGW]], [[NSPOOLGW]] - include this line only if NWS is not equal to zero.<br />
<br />
====Harmonic Analysis====<br />
[[NFREQ]]<br/><br />
for k=1 to [[NFREQ]]<br/><br />
: [[NAMEFR(k)]]<br/><br />
: [[HAFREQ(k), HAFF(k), HAFACE(k)]]<br/><br />
end k loop<br/><br />
[[THAS]], [[THAF]], [[NHAINC]], [[FMV]]<br/><br />
[[NHASE]], [[NHASV]], [[NHAGE]], [[NHAGV]]<br />
<br />
====Hotstart Output and Numeric Controls====<br />
[[NHSTAR]], [[NHSINC]]<br/><br />
[[ITITER]], [[ISLDIA]], [[CONVCR]], [[ITMAX]]<br/><br />
<br />
''For a 2DDI ADCIRC run that does not use netCDF nor namelists, the file ends here. For those controls, see further below in the [[#NetCDF Controls|NetCDF Controls]] and [[#Namelists|Namelists]] sections.''<br />
<br />
==3D Model Run==<br />
[[IDEN]]<br/><br />
[[ISLIP]], [[KP]]<br/><br />
[[Z0S,Z0B]]<br/><br />
[[ALP1,ALP2,ALP3]]<br/><br />
[[IGC]], [[NFEN]]<br/><br />
for k=1 to [[NFEN]] (include this loop only if [[IGC]] = 0, k=1 at bottom, k= [[NFEN]] at surface)<br/><br />
: [[SIGMA(k)]]<br/><br />
end k loop<br/><br />
[[IEVC]], [[EVMIN]], [[EVCON]]<br/><br />
for k=1 to [[NFEN]] (include this loop only if [[IEVC]] = 0, k=1 at bottom, k= [[NFEN]] at surface)<br/><br />
: [[EVTOT(k)]]<br/><br />
end k loop<br/><br />
[[THETA1, THETA2]](include this line only if [[IEVC]] = 50 or 51)<br/><br />
[[I3DSD,TO3DSDS,TO3DSDF,NSPO3DSD]]<br/><br />
[[NSTA3DD]]<br/><br />
for k=1 to [[NSTA3DD]]<br/><br />
: [[X3DS(k), Y3DS(k)]]<br/><br />
end k loop<br/><br />
[[I3DSV,TO3DSVS,TO3DSVF,NSPO3DSV]]<br/><br />
[[NSTA3DV]]<br/><br />
for k=1 to [[NSTA3DV]]<br/><br />
: [[X3DS(k), Y3DS(k)]]<br/><br />
end k loop<br/><br />
[[I3DST,TO3DSTS,TO3DSTF,NSPO3DST]]<br/><br />
[[NSTA3DT]]<br/><br />
for k=1 to [[NSTA3DT]]<br/><br />
: [[X3DS(k), Y3DS(k)]]<br/><br />
end k loop<br/><br />
[[I3DGD]],[[TO3DGDS]],[[TO3DGDF]],[[NSPO3DGD]]<br/><br />
[[I3DGV]],[[TO3DGVS]],[[TO3DGVF]],[[NSPO3DGV]]<br/><br />
[[I3DGT]],[[TO3DGTS]],[[TO3DGTF]],[[NSPO3DGT]]<br/><br />
The following line will be read in if [[IM]] is 21 or 31.<br/><br />
[[RES_BC_FLAG]], [[BCFLAG_LNM]], [[BCFLAG_TEMP]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] is negative.<br/><br />
[[RBCTIMEINC]]<br/><br />
[[BCSTATIM]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] = 2.<br/><br />
[[RBCTIMEINC]], [[SBCTIMEINC]]<br/><br />
[[BCSTATIM]], [[SBCSTATIM]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] = 3.<br/><br />
[[RBCTIMEINC]], [[TBCTIMEINC]]<br/><br />
[[BCSTATIM]], [[TBCSTATIM]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] = 4.<br/><br />
[[RBCTIMEINC]], [[SBCTIMEINC]], [[TBCTIMEINC]]<br/><br />
[[BCSTATIM]], [[SBCSTATIM]], [[TBCSTATIM]]<br/><br />
The following two lines will be read in if [[RES_BC_FLAG]] = 3 or 4 and [[BCFLAG_TEMP]] is not equal to 0.<br/><br />
[[TTBCTIMEINC]], [[TTBCSTATIM]]<br/><br />
[[TTBCTIMEINC]]<br/><br />
The following two lines will be read in only if [[IM]] is 21 or 31.<br/><br />
[[SPONGEDIST]]<br/><br />
[[EQNSTATE]]<br/><br />
The following lines will be read in only if [[IDEN]] is > 0.<br/><br />
[[NLSD, NVSD]]<br/><br />
[[NLTD, NVTD]]<br/><br />
[[ALP4]]<br/><br />
The following line will be read in only if [[IDEN]] = 3 or 4.<br/><br />
[[NTF]]<br />
<br />
==NetCDF Controls==<br />
The following lines will be read in only if the NetCDF output or hotstart format is chosen<br/><br />
NCPROJ<br/><br />
NCINST<br/><br />
NCSOUR<br/><br />
NCHIST<br/><br />
NCREF<br/><br />
NCCOM<br/><br />
NCHOST<br/><br />
NCCONV<br/><br />
NCCONT<br/><br />
NCDATE<br />
<br />
==Namelists==<br />
The following Fortran namelist lines are optional, but if they appear, they must appear at the very end of the fort.15 file.<br/><br />
<code>&metControl WindDragLimit=floatValue, DragLawString='stringValue', rhoAir=floatValue, outputWindDrag=logicalValue /</code><br/><br />
<code>&timeBathyControl NDDT=integerValue, BTIMINC=floatValue, BCHGTIMINC=floatValue /</code><br/><br />
<code>&waveCoupling WindWaveMultiplier=floatValue /</code><br/><br />
<code>&SWANOutputControl SWAN_OutputHS=logicalValue, SWAN_OutputDIR=logicalValue, SWAN_OutputTM01=logicalValue, SWAN_OutputTPS=logicalValue, SWAN_OutputWIND=logicalValue, SWAN_OutputTM02=logicalValue, SWAN_OutputTMM10=logicalValue /</code><br/><br />
<code>&subdomainModeling subdomainOn=logicalValue/</code><br/><br />
<code>&wetDryControl outputNodeCode=logicalValue, outputNOFF=logicalValue, noffActive=logicalValue /</code><br/><br />
<code>&inundationOutputControl inundationOutput=logicalValue0, inunThresh =floatValue /</code><br/><br />
<code>&TVWControl use_TVW=logicalValue, TVW_file='stringValue', nout_TVW =integerValue, touts_TVW =floatValue, toutf_TVW=floatValue, nspool_TVW =integerValue /</code><br/><br />
<code>&WarnElevControl WarnElev=floatValue, ErrorElev=floatValue, WarnElevDump=logicalValue, WarnElevDumpLimit=integerValue /</code><br/><br />
<code>[[Dynamic_water_level_correction#Controlling_Water_Level_Correction|&dynamicWaterLevelCorrectionControl]] dynamicWaterLevelCorrectionFileName='stringValue' dynamicWaterLevelCorrectionMultiplier=floatValue, dynamicWaterLevelCorrectionRampStart=floatValue, dynamicWaterLevelCorrectionRampEnd=floatValue, dynamicWaterLevelCorrectionRampReferenceTime='stringValue', dynamicWaterLevelCorrectionSkipSnaps=integerValue /</code><br/><br />
<code>&AliDispersionControl CAliDisp=logicalValue, Cs=floatValue, Ad=floatValue, Bd=floatValue /</code><br/><br />
[[category:input files]]</div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1055IM2020-07-10T14:24:20Z<p>Wpringle: /* Six-digit IM Codes */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111113<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611113<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| A value of 4-6 does the same as 1-3 (same order) but in baroclinic mode<br/><code>CBaroclinic=.TRUE.</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| See above <br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
The most recent version (55+) also has an option that improves the (default) semi-implicit consistent GWCE mass matrix mode to compute the complete (total depth) gravity wave term (free surface gradient) implicitly; toggled by setting IMDigit-6 to 3. The default version (IMDigit-6=1), only computes the initial still water depth component of the free surface gradient implicitly, which might make it more susceptible to CFL violations in shallow depths and can encounter Matrix diagonality issues overland where the initial still water depth is negative.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1054IM2020-07-10T14:23:27Z<p>Wpringle: /* Six-digit IM Codes */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111113<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611113<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| A value of 4-6 does the same as 1-3 (same order) but in baroclinic mode<br/><code>CBaroclinic=.TRUE.</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| See above <br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
The most recent version (55+) also has an option that improves the (default) semi-implicit consistent GWCE mass matrix mode to compute the complete (total depth) gravity wave term (free surface gradient) implicitly; toggled by setting IMDig6 to 3. The default version (IMDig=1), only computes the initial still water depth component of the free surface gradient implicitly, which might make it more susceptible to CFL violations in shallow depths and can encounter Matrix diagonality issues overland where the initial still water depth is negative.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1053IM2020-07-10T14:22:59Z<p>Wpringle: /* Six-digit IM Codes */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111113<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611113<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| A value of 4-6 does the same as 1-3 (same order) but in baroclinic mode<br/><code>CBaroclinic=.TRUE.</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| See above <br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
The most recent version (55+) also has an option that improves the (default) semi-implicit consistent GWCE mass matrix mode to include compute the complete (total depth) gravity wave term (free surface gradient) implicitly; toggled by setting IMDig6 to 3. The default version (IMDig=1), only computes the initial still water depth component of the free surface gradient implicitly, which might make it more susceptible to CFL violations in shallow depths and can encounter Matrix diagonality issues overland where the initial still water depth is negative.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1052IM2020-07-10T14:10:11Z<p>Wpringle: </p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>; users should be aware that the <code>[[A00, B00, C00]]</code> coefficients must be specified differently in these two cases. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111113<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611113<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| A value of 4-6 does the same as 1-3 (same order) but in baroclinic mode<br/><code>CBaroclinic=.TRUE.</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| See above <br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=IM&diff=1051IM2020-07-10T14:07:31Z<p>Wpringle: /* Six-digit IM Codes */</p>
<hr />
<div>'''<code>IM</code>''' is an important parameter in the [[fort.15 file]] that defines numerical model formulation and dimension. Among other things, <code>IM</code> specifies whether ADCIRC is solved in two-dimensional depth-integrated (2DDI) or in three-dimensions (3D), solution of the governing equations is semi-implicit or explicit in time, and whether the model formulation is barotropic or baroclinic. Popular values for 2D barotropic ADCIRC include <code>IM=0</code> and <code>IM=111112</code>, though the latter also requires modifying <code>[[A00, B00, C00]]</code>. <br />
<br />
== Default IM Values ==<br />
Default simulation option combinations can be specified through single or double digit values, some of which are shortcuts to the six-digit codes described in the next heading. <br />
The available <code>IM</code> values are specified in the table below and in the following section on 6-digit values:<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! IM Value<br />
! Six-digit Equivalent<br />
! Description<br />
|-<br />
| 0<br />
| 111111<br />
| Barotropic 2DDI <br />
|-<br />
| 1<br />
| 611111<br />
| Barotropic 3D velocity-based momentum<br />
|-<br />
| 2<br />
| -<br />
| Barotropic 3D stress-based momentum<br />
|-<br />
| 10<br />
| -<br />
| Barotropic 2DDI with passive scalar transport<br />
|-<br />
| 11<br />
| -<br />
| Barotropic 3D velocity-based momentum with passive scalar transport<br />
|-<br />
| 20<br />
| 111113<br />
| Baroclinic 2DDI<br />
|-<br />
| 21<br />
| 611113<br />
| Baroclinic 3D velocity-based momentum<br />
|-<br />
| 30<br />
| -<br />
| Baroclinic 2DDI with passive scalar transport<br />
|-<br />
| 31<br />
| -<br />
| Baroclinic 3D velocity-based momentum with passive scalar transport<br />
|}<br />
<br />
Note that all default <code>IM</code> values employ the semi-implicit consistent GWCE mass matrix solver. It has less numerical error and tends to be more stable than the explicit mass-lumping approach at the expense of computational time and memory.<br />
<br />
== Six-digit IM Codes ==<br />
For fine-grained control of various options six-digit codes for <code>IM</code> can be specified. Each digit represents a specific option regarding the dimension and the formulation of certain terms or integration methods in the GWCE or momentum equations. <br />
The available options for each digit are specified below, with the first digit being the left-most. The internal flags that are set are listed to help users dig through the code. <br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! Value<br />
! Digit 1: 2DDI/3D, Lateral Stress in GWCE<ref name=Kendra1>K.M. Dresback, R.L. Kolar, R.A. Luettich, Jr. (2005). On the Form of the Momentum Equation and Lateral Stress Closure Law in Shallow Water Modeling, in: Estuar. Coast. Model., American Society of Civil Engineers, Reston, VA, 399–418. doi:10.1061/40876(209)23</ref><br />
! Digit 2: Advection in GWCE<ref name=Kendra2>K.M. Dresback, R.L. Kolar, J.C. Dietrich (2005). On the Form of the Momentum Equation for Shallow Water Models Based on the Generalized Wave Continuity Equation: Conservative vs. Non-Conservative. Advances in Water Resources, 28(4), 345-358. doi:10.1016/j.advwatres.2004.11.011</ref><br />
! Digit 3: Lateral Stress in Momentum<ref name=Kendra1></ref><br />
! Digit 4: Advection in Momentum<ref name=Kendra2></ref><br />
! Digit 5: Area Integration in Momentum<br />
! Digit 6: GWCE Mass Matrix, Barotropic/Baroclinic<br />
|-<br />
| 1 (default)<br />
| 2DDI, Kolar-Gray flux-based<br/><code>CGWCE_LS_KGQ=.TRUE.</code><br />
| Non conservative<br/><code>CGWCE_Advec_NC=.TRUE.</code><br />
| Integration by parts, velocity-based<br/><code>CME_LS_IBPV=.TRUE.</code><br />
| Non conservative<br/><code>CME_New_NC=.TRUE.</code><br />
| Corrected <br/><code>CME_AreaInt_Corr=.TRUE.</code><br />
| Consistent (implicit for linear part of gravity wave term), barotropic<br/><code>ILump=0</code><br />
|-<br />
| 2<br />
| 2DDI, 2-part flux-based<br/><code>CGWCE_LS_2PartQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CGWCE_Advec_C1=.TRUE.</code><br />
| Integration by parts, flux-based<br/><code>CME_LS_IBPQ=.TRUE.</code><br />
| Conservative form 1<br/><code>CME_New_C1=.TRUE.</code><br />
| Original <br/><code>CME_AreaInt_Orig=.TRUE.</code><br />
| Lumped (explicit), barotropic<br/><code>CGWCE_Lump=.TRUE.</code>, <code>ILump=1</code><br />
|-<br />
| 3<br />
| 2DDI, 2-part velocity-based<br/><code>CGWCE_LS_2PartV=.TRUE.</code><br />
| Conservative form 2<br/><code>CGWCE_Advec_C2=.TRUE.</code><br />
| Integration by parts, velocity-based symmetrical<br/><code>CME_LS_IBPSV=.TRUE.</code><br />
| Conservative form 2<br/><code>CME_New_C2=.TRUE.</code><br />
| -<br />
| Consistent (implicit for full gravity wave term), barotropic<br/><code>CGWCE_HDP=.TRUE.</code>, <code>IFNL_HDP=1</code>, <code>ILump=0</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 4<br />
| 2DDI, 2-part flux-based symmetrical<br/><code>CGWCE_LS_2PartSQ=.TRUE.</code><br />
| -<br />
| Integration by parts, flux-based symmetrical<br/><code>CME_LS_IBPSQ=.TRUE.</code><br />
| -<br />
| -<br />
| A value of 4-6 does the same as 1-3 (same order) but in baroclinic mode<br/><code>CBaroclinic=.TRUE.</code> {{ADC version|version=55|relation=eq}}<br />
|-<br />
| 5<br />
| 2DDI, 2-part velocity-based symmetrical<br/><code>CGWCE_LS_2PartSV=.TRUE.</code><br />
| -<br />
| 2 Part, velocity-based (''not implemented'')<br/><code>CME_LS_2PartV=.TRUE.</code><br />
| -<br />
| -<br />
| -<br />
|-<br />
| 6<br />
| 3D, Kolar-Gray flux-based<br/><code>C2DDI=.FALSE.</code>, <code>CGWCE_LS_KGQ=.TRUE.</code>, <code>C3D=.TRUE.</code>, <code>C3DVS=.TRUE.</code>, <code>ILump=0</code><br />
| -<br />
| 2 Part, flux-based (''not implemented'')<br/><code>CME_LS_2PartQ=.TRUE.</code><br />
| -<br />
| -<br />
| See above <br />
|}<br />
<br />
A common code combination is <code>IM=111112</code>, which is identical to the default <code>111111</code> (same as <code>IM=0</code>), but simulates in explicit mass-lumping mode. Note that <code>[[A00, B00, C00]]</code> must be set to <code>0.0 1.0 0.0</code> when in this mode. Lumped explicit mode is a useful alternative to the (default) semi-implicit consistent GWCE mass matrix mode, because the latter requires a matrix solve that increases computational time and memory. By comparison, the explicit mass-lumping mode is about twice as fast and scales to fewer grid nodes per computational core.<ref>S. Tanaka, S. Bunya, J.J. Westerink, C. Dawson, R.A. Luettich, Scalability of an Unstructured Grid Continuous Galerkin Based Hurricane Storm Surge Model, J. Sci. Comput. 46 (2011) 329–358. doi:10.1007/s10915-010-9402-1</ref> Moreover, for model setups that are sufficiently resolved in space and time, differences in the solution between approaches should be small. Though, many users have reported somewhat lower stability in lumped explicit mode.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Global_Astronomical_M2_Tide&diff=1036Global Astronomical M2 Tide2020-07-06T21:37:34Z<p>Wpringle: /* Options/Features Tested */</p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of the astronomical M2 tidal constituent on the spherical Earth under equilibrium tidal forcing with the inclusion of the self-attraction and loading tide. The results of interest are the M2 tidal constituent amplitudes and phases of elevations and velocities from the least-squares harmonic analysis. The test finishes in about 5 minutes in serial ADCIRC for a full month of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_global-tide-2d GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh is a coarse representation of the spherical Earth with minimum resolution of approximately 50 km, comprised of 27,330 vertices and 50,859 triangular elements. <br />
<br />
== Options/Features Tested ==<br />
*<code>[[ICS]]</code> = -22: Uses the Mercator projection with a coordinate rotation to remove the pole singularity (need to provide a [[fort.rotm]]). <br />
*<code>[[IM]]</code> = 513113: Uses the fully implicit scheme for the gravity wave term (computational time step is 12 minutes). <br />
*<code>[[NTIP]]</code> = 2: equilibrium tide + self-attraction and loading tide forcing (read from a [[fort.24 file]]).<br />
*<code>[[A00, B00, C00]]</code> = 0.5, 0.5, 0: Ensures that the fully implicit scheme is stable with a large time step.<br />
*<code>[[ESLM]]</code> = -0.2: enables the Smagorinsky turbulence closure with a coefficient of 0.2.<br />
*<code>[[NHAGE]]</code> = 5: outputs the harmonic constituent elevations into a netCDF4 [[fort.53 file]]. <br />
*<code>[[NHAGV]]</code> = 5: outputs the harmonic constituent velocities into a netCDF4 [[fort.54 file]]. <br />
*[[fort.13_file#Internal_Tide_Energy_Conversion|internal_tide_friction]]: spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.<br />
*[[fort.13_file#Quadratic_Friction_coefficient|quadratic_friction_coefficient]]: spatially varying quadratic bottom friction [[fort.13 file]] attribute.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=Fort.200_file_format&diff=1027Fort.200 file format2020-06-29T23:13:53Z<p>Wpringle: /* File Structure */</p>
<hr />
<div>The basic file structure is shown below. Each line of input data is represented by a line containing the input variable name(s). Blank lines are only to enhance readability. Loops indicate multiple lines of input. Conditional input is indicated by an if clause. <br />
<br />
=NWS = 10 or 110 or 10010 =<br />
<br />
==File Structure==<br />
for k=1, [[LONB*LATB]]<br />
<br />
: [[PG(k)]], [[UG(k)]], [[VG(k)]] , ''if NWS=10010:'' [[IG(k)]] <br />
<br />
end k loop<br />
<br />
'''Key:'''<br /><br />
[[PG(k)]]: surface pressure (Pa)<br /><br />
[[UG(k)]]: U10 wind speed (m/s)<br /><br />
[[VG(k)]]: V10 wind speed (m/s)<br /><br />
{{ADC version|version=55|relation=ge}} [[IG(k)]]: ice concentration (0-1 decimal fraction)<br />
<br />
==Notes==<br />
NWS = 10, US National Climatic Data Center (NCDC) Global Forecast System (GFS) model<br />
<br />
Meteorological data is input to a T1534, Gaussian grid (NWLON*NWLAT is hardwired to 3072*1536 = 4,718,592) and interpolated in space onto the ADCIRC mesh. Data are ordered in these files the same way they are generated by the degribbing program.<br />
<br />
The met data are contained in a sequence of files with names as follows: fort.200, fort.200+N, fort.200+2*N, fort.200+3*N,…., where N is the time interval (in hours) between successive meteorological data (e.g., at 6 hr intervals this would be fort.200, fort.206, fort.212, etc.). N is determined as N = WTIMINC/3600 where WTIMINC is the meteorological data time interval (in seconds) and is specified in the Model Parameter and Periodic Boundary Condition File. WTIMINC must be evenly divisible by 3600.<br />
<br />
If the model is cold started, it is assumed that the winds are at rest at the beginning of the model run and no data file (fort.200) is read corresponding to TIME=STATIM. Rather, the first required data file (fort.200+N) corresponds to TIME=STATIM+WTIMINC. If the model is hot started, an initial file (fort.200) is required corresponding to TIME=HOT START TIME. In either case, additional sets of meteorological data must be provided every WTIMINC, in appropriately named files. Meteorological data is interpolated in time to the ADCIRC time step.<br />
<br />
Wind velocity (@ 10 m above the water surface) must be input in units of m/s. Surface atmospheric pressure must be input in units of Pascals = Newtons/square meter.<br />
<br />
The following relations are used to compute wind stress from the input wind velocity.<br />
<br />
WIND_SPEED = magnitude of WIND_VEL<br />
<br />
DRAG_COEFF = 0.001*(0.75+0.067*WIND_SPEED)<br />
<br />
If (DRAG_COEFF.gt.0.003) DRAG_COEFF=0.003<br />
<br />
WIND_STRESS = DRAG_COEFF*0.001293*WIND_VEL*WIND_SPEED<br />
<br />
The following relationship is used in ADCIRC to convert to pressure in meters of water from pressure in Pascal:<br />
<br />
PRESSURE{m H2O}=PRESSURE{Pascal}/(GRAVITY*DENSITY of H2O)<br />
<br />
=NWS = 11 or 111=<br />
<br />
==File Structure==<br />
for k=1, 8<br />
<br />
: IYEAR, IMONTH, IDAY, IHOUR<br />
<br />
: for j=1, [[LONB*LATB]]<br />
<br />
:: [[UE(j)]], [[VE(j)]], [[PE(j)]]<br />
<br />
: end j loop<br />
<br />
end k loop<br />
<br />
==Notes==<br />
NWS = 11, 111 – National Weather Service Eta-29 file.<br />
<br />
Meteorological data is input to the standard Eta-29 “E” grid, (non-cartesian grid, NWLON*NWLAT is hardwired to 181*271 = 49,051), and interpolated in space onto the ADCIRC grid. The ADCIRC grid must be in lon, lat coordinates.<br />
<br />
The met data are contained in a sequence of files with names: fort.200, fort.201, fort.202, fort.203,…., Each file contains 8 data sets spaced 3 hrs apart in time (WTIMINC is the met. data time interval = 3600*3 = 10800 sec) for a total of 1 day of data per file. The first and last data sets in each file are assumed to correspond to 03:00 hrs and 24:00 hrs, respectively, on the day of the file. If the model is cold started, it is assumed that the winds are at rest at the beginning of the model run and no data file (fort.200) is read corresponding to TIME=STATIM. Rather, the initial data set in the first required data file (fort.201) corresponds to TIME=STATIM+WTIMINC. If the model is hot started, an initial file (fort.200) is required in which the last data set corresponds to TIME=HOT START TIME. In either case, 8 sets of met. data must be provided every<br />
WTIMINC in each file and appropriately named files must be provided for each day of the model run. Met data is interpolated in time to the ADCIRC time step.<br />
<br />
The file type is binary: ACCESS=’sequential’, FORM=’unformatted’<br />
<br />
Wind velocity (@ 10 m above the water surface) must be input in units of m/s oriented along the E grid coordinate directions. Surface atmospheric pressure must be input in units of millibars.<br />
<br />
The following relations are used to compute wind stress from the input wind velocity.<br />
<br />
WIND_SPEED = magnitude of WIND_VEL<br />
<br />
DRAG_COEFF = 0.001*(0.75+0.067*WIND_SPEED)<br />
<br />
If (DRAG_COEFF.gt.0.003) DRAG_COEFF=0.003<br />
<br />
WIND_STRESS = DRAG_COEFF*0.001293*WIND_VEL*WIND_SPEED<br />
<br />
The following relationship is used in ADCIRC to convert to pressure in meters of water from pressure in millibars.<br />
<br />
PRESSURE{m H2O}=PRESSURE{millibars}/(100*GRAVITY*DENSITY of H2O)<br />
<br />
The integer values IYEAR, IMONTH, IDAY, IHOUR are read in and ignored by ADCIRC.<br />
<br />
=General Note=<br />
<br />
Meteorological data must be provided for the entire model run, otherwise the run will crash.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=Fort.200_file_format&diff=1026Fort.200 file format2020-06-29T23:12:57Z<p>Wpringle: /* NWS = 10 or 110 */</p>
<hr />
<div>The basic file structure is shown below. Each line of input data is represented by a line containing the input variable name(s). Blank lines are only to enhance readability. Loops indicate multiple lines of input. Conditional input is indicated by an if clause. <br />
<br />
=NWS = 10 or 110 or 10010 =<br />
<br />
==File Structure==<br />
for k=1, [[LONB*LATB]]<br />
<br />
: [[PG(k)]], [[UG(k)]], [[VG(k)]] , ''if NWS=10010:'' [[IG(k)]] <br />
<br />
end k loop<br />
<br />
Key:<br /><br />
[[PG(k)]]: surface pressure (Pa)<br /><br />
[[UG(k)]]: U10 wind speed (m/s)<br /><br />
[[VG(k)]]: V10 wind speed (m/s)<br /><br />
{{ADC version|version=55|relation=ge}} [[IG(k)]]: ice concentration (0-1 decimal fraction) <br />
<br />
==Notes==<br />
NWS = 10, US National Climatic Data Center (NCDC) Global Forecast System (GFS) model<br />
<br />
Meteorological data is input to a T1534, Gaussian grid (NWLON*NWLAT is hardwired to 3072*1536 = 4,718,592) and interpolated in space onto the ADCIRC mesh. Data are ordered in these files the same way they are generated by the degribbing program.<br />
<br />
The met data are contained in a sequence of files with names as follows: fort.200, fort.200+N, fort.200+2*N, fort.200+3*N,…., where N is the time interval (in hours) between successive meteorological data (e.g., at 6 hr intervals this would be fort.200, fort.206, fort.212, etc.). N is determined as N = WTIMINC/3600 where WTIMINC is the meteorological data time interval (in seconds) and is specified in the Model Parameter and Periodic Boundary Condition File. WTIMINC must be evenly divisible by 3600.<br />
<br />
If the model is cold started, it is assumed that the winds are at rest at the beginning of the model run and no data file (fort.200) is read corresponding to TIME=STATIM. Rather, the first required data file (fort.200+N) corresponds to TIME=STATIM+WTIMINC. If the model is hot started, an initial file (fort.200) is required corresponding to TIME=HOT START TIME. In either case, additional sets of meteorological data must be provided every WTIMINC, in appropriately named files. Meteorological data is interpolated in time to the ADCIRC time step.<br />
<br />
Wind velocity (@ 10 m above the water surface) must be input in units of m/s. Surface atmospheric pressure must be input in units of Pascals = Newtons/square meter.<br />
<br />
The following relations are used to compute wind stress from the input wind velocity.<br />
<br />
WIND_SPEED = magnitude of WIND_VEL<br />
<br />
DRAG_COEFF = 0.001*(0.75+0.067*WIND_SPEED)<br />
<br />
If (DRAG_COEFF.gt.0.003) DRAG_COEFF=0.003<br />
<br />
WIND_STRESS = DRAG_COEFF*0.001293*WIND_VEL*WIND_SPEED<br />
<br />
The following relationship is used in ADCIRC to convert to pressure in meters of water from pressure in Pascal:<br />
<br />
PRESSURE{m H2O}=PRESSURE{Pascal}/(GRAVITY*DENSITY of H2O)<br />
<br />
=NWS = 11 or 111=<br />
<br />
==File Structure==<br />
for k=1, 8<br />
<br />
: IYEAR, IMONTH, IDAY, IHOUR<br />
<br />
: for j=1, [[LONB*LATB]]<br />
<br />
:: [[UE(j)]], [[VE(j)]], [[PE(j)]]<br />
<br />
: end j loop<br />
<br />
end k loop<br />
<br />
==Notes==<br />
NWS = 11, 111 – National Weather Service Eta-29 file.<br />
<br />
Meteorological data is input to the standard Eta-29 “E” grid, (non-cartesian grid, NWLON*NWLAT is hardwired to 181*271 = 49,051), and interpolated in space onto the ADCIRC grid. The ADCIRC grid must be in lon, lat coordinates.<br />
<br />
The met data are contained in a sequence of files with names: fort.200, fort.201, fort.202, fort.203,…., Each file contains 8 data sets spaced 3 hrs apart in time (WTIMINC is the met. data time interval = 3600*3 = 10800 sec) for a total of 1 day of data per file. The first and last data sets in each file are assumed to correspond to 03:00 hrs and 24:00 hrs, respectively, on the day of the file. If the model is cold started, it is assumed that the winds are at rest at the beginning of the model run and no data file (fort.200) is read corresponding to TIME=STATIM. Rather, the initial data set in the first required data file (fort.201) corresponds to TIME=STATIM+WTIMINC. If the model is hot started, an initial file (fort.200) is required in which the last data set corresponds to TIME=HOT START TIME. In either case, 8 sets of met. data must be provided every<br />
WTIMINC in each file and appropriately named files must be provided for each day of the model run. Met data is interpolated in time to the ADCIRC time step.<br />
<br />
The file type is binary: ACCESS=’sequential’, FORM=’unformatted’<br />
<br />
Wind velocity (@ 10 m above the water surface) must be input in units of m/s oriented along the E grid coordinate directions. Surface atmospheric pressure must be input in units of millibars.<br />
<br />
The following relations are used to compute wind stress from the input wind velocity.<br />
<br />
WIND_SPEED = magnitude of WIND_VEL<br />
<br />
DRAG_COEFF = 0.001*(0.75+0.067*WIND_SPEED)<br />
<br />
If (DRAG_COEFF.gt.0.003) DRAG_COEFF=0.003<br />
<br />
WIND_STRESS = DRAG_COEFF*0.001293*WIND_VEL*WIND_SPEED<br />
<br />
The following relationship is used in ADCIRC to convert to pressure in meters of water from pressure in millibars.<br />
<br />
PRESSURE{m H2O}=PRESSURE{millibars}/(100*GRAVITY*DENSITY of H2O)<br />
<br />
The integer values IYEAR, IMONTH, IDAY, IHOUR are read in and ignored by ADCIRC.<br />
<br />
=General Note=<br />
<br />
Meteorological data must be provided for the entire model run, otherwise the run will crash.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=Fort.200_file_format&diff=1025Fort.200 file format2020-06-29T23:07:10Z<p>Wpringle: /* File Structure */</p>
<hr />
<div>The basic file structure is shown below. Each line of input data is represented by a line containing the input variable name(s). Blank lines are only to enhance readability. Loops indicate multiple lines of input. Conditional input is indicated by an if clause. <br />
<br />
=NWS = 10 or 110=<br />
<br />
==File Structure==<br />
for k=1, [[LONB*LATB]]<br />
<br />
: [[PG(k)]], [[UG(k)]], [[VG(k)]] , ''optional:'' [[IG(k)]] <br />
<br />
end k loop<br />
<br />
==Notes==<br />
NWS = 10, US National Climatic Data Center (NCDC) Global Forecast System (GFS) model<br />
<br />
Meteorological data is input to a T1534, Gaussian grid (NWLON*NWLAT is hardwired to 3072*1536 = 4,718,592) and interpolated in space onto the ADCIRC mesh. Data are ordered in these files the same way they are generated by the degribbing program.<br />
<br />
The met data are contained in a sequence of files with names as follows: fort.200, fort.200+N, fort.200+2*N, fort.200+3*N,…., where N is the time interval (in hours) between successive meteorological data (e.g., at 6 hr intervals this would be fort.200, fort.206, fort.212, etc.). N is determined as N = WTIMINC/3600 where WTIMINC is the meteorological data time interval (in seconds) and is specified in the Model Parameter and Periodic Boundary Condition File. WTIMINC must be evenly divisible by 3600.<br />
<br />
If the model is cold started, it is assumed that the winds are at rest at the beginning of the model run and no data file (fort.200) is read corresponding to TIME=STATIM. Rather, the first required data file (fort.200+N) corresponds to TIME=STATIM+WTIMINC. If the model is hot started, an initial file (fort.200) is required corresponding to TIME=HOT START TIME. In either case, additional sets of meteorological data must be provided every WTIMINC, in appropriately named files. Meteorological data is interpolated in time to the ADCIRC time step.<br />
<br />
Wind velocity (@ 10 m above the water surface) must be input in units of m/s. Surface atmospheric pressure must be input in units of Pascals = Newtons/square meter.<br />
<br />
The following relations are used to compute wind stress from the input wind velocity.<br />
<br />
WIND_SPEED = magnitude of WIND_VEL<br />
<br />
DRAG_COEFF = 0.001*(0.75+0.067*WIND_SPEED)<br />
<br />
If (DRAG_COEFF.gt.0.003) DRAG_COEFF=0.003<br />
<br />
WIND_STRESS = DRAG_COEFF*0.001293*WIND_VEL*WIND_SPEED<br />
<br />
The following relationship is used in ADCIRC to convert to pressure in meters of water from pressure in Pascal:<br />
<br />
PRESSURE{m H2O}=PRESSURE{Pascal}/(GRAVITY*DENSITY of H2O)<br />
<br />
=NWS = 11 or 111=<br />
<br />
==File Structure==<br />
for k=1, 8<br />
<br />
: IYEAR, IMONTH, IDAY, IHOUR<br />
<br />
: for j=1, [[LONB*LATB]]<br />
<br />
:: [[UE(j)]], [[VE(j)]], [[PE(j)]]<br />
<br />
: end j loop<br />
<br />
end k loop<br />
<br />
==Notes==<br />
NWS = 11, 111 – National Weather Service Eta-29 file.<br />
<br />
Meteorological data is input to the standard Eta-29 “E” grid, (non-cartesian grid, NWLON*NWLAT is hardwired to 181*271 = 49,051), and interpolated in space onto the ADCIRC grid. The ADCIRC grid must be in lon, lat coordinates.<br />
<br />
The met data are contained in a sequence of files with names: fort.200, fort.201, fort.202, fort.203,…., Each file contains 8 data sets spaced 3 hrs apart in time (WTIMINC is the met. data time interval = 3600*3 = 10800 sec) for a total of 1 day of data per file. The first and last data sets in each file are assumed to correspond to 03:00 hrs and 24:00 hrs, respectively, on the day of the file. If the model is cold started, it is assumed that the winds are at rest at the beginning of the model run and no data file (fort.200) is read corresponding to TIME=STATIM. Rather, the initial data set in the first required data file (fort.201) corresponds to TIME=STATIM+WTIMINC. If the model is hot started, an initial file (fort.200) is required in which the last data set corresponds to TIME=HOT START TIME. In either case, 8 sets of met. data must be provided every<br />
WTIMINC in each file and appropriately named files must be provided for each day of the model run. Met data is interpolated in time to the ADCIRC time step.<br />
<br />
The file type is binary: ACCESS=’sequential’, FORM=’unformatted’<br />
<br />
Wind velocity (@ 10 m above the water surface) must be input in units of m/s oriented along the E grid coordinate directions. Surface atmospheric pressure must be input in units of millibars.<br />
<br />
The following relations are used to compute wind stress from the input wind velocity.<br />
<br />
WIND_SPEED = magnitude of WIND_VEL<br />
<br />
DRAG_COEFF = 0.001*(0.75+0.067*WIND_SPEED)<br />
<br />
If (DRAG_COEFF.gt.0.003) DRAG_COEFF=0.003<br />
<br />
WIND_STRESS = DRAG_COEFF*0.001293*WIND_VEL*WIND_SPEED<br />
<br />
The following relationship is used in ADCIRC to convert to pressure in meters of water from pressure in millibars.<br />
<br />
PRESSURE{m H2O}=PRESSURE{millibars}/(100*GRAVITY*DENSITY of H2O)<br />
<br />
The integer values IYEAR, IMONTH, IDAY, IHOUR are read in and ignored by ADCIRC.<br />
<br />
=General Note=<br />
<br />
Meteorological data must be provided for the entire model run, otherwise the run will crash.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=NWS&diff=1024NWS2020-06-29T23:04:35Z<p>Wpringle: /* Extended NWS with Ice + Waves */</p>
<hr />
<div>'''<code>NWS</code>''' is a parameter in the [[fort.15 file]] that selects the meteorological forcing input type. The value on the "[[fort.15_file_format#NWS|NWS line]]" of the fort.15 file also implicitly includes [[#Value Seen in fort.15 File|other parameters]] affecting wave coupling and ice inputs. Further, <code>NWS</code> affects not just the file type and handling of meteorological data, but also changes what the [[Fort.15_file_format#WTIMINC|meteorological parameter line]] (informally, the <code>[[WTIMINC]]</code> line) looks like in the fort.15 file. See the [[supplemental meteorological/wave/ice parameters]] page for information on the format of this line. <br />
<br />
ADCIRC supports a wide range of meteorological input formats, including moving/fixed gridded data in several file formats, tropical cyclone track and parameter data that can be turned into wind/pressure fields via one of several internal vortex models, and direct specification of wind speeds or stresses on nodes. As a result of the great flexibility and importance of this choice, several pages are devoted to the topic. In particular, see also the [[fort.22 file]] and [[wind stress]] pages. <br />
<br />
== Value Seen in fort.15 File ==<br />
In the fort.15 file, what we call "NWS" is actually a combination of several parameters. For example, given a 5-digit value on that line, <br />
-12305 ! TRICKY NWS IMPOSTER<br />
the first two digits (ten-thousands and thousands) tell us the format of ice data <code>[[NCICE]]=12</code>, the 3rd digit (hundreds) tells us the wave coupling mode <code>[[NRS]]=3</code>, and the last two digits (tens and ones) combined with the sign of the entire value tell us the meteorological forcing mode <code>NWS=-5</code>. If the value has only 3 digits then ADCIRC assumes no ice input <code>[[NCICE]]=0</code>, and if it's 2 digits then ADCIRC further assumes no wave coupling <code>[[NRS]]=0</code>. It is often presumed that when one refers to "NWS", one is referring to the 2-digit value, not what is in the fort.15 file, but be mindful of the ambiguity here. <br />
<br />
== Parameter Summary ==<br />
The following table is a summary of possible <code>NWS</code> values, their descriptions, and associated meteorological input files (required and optional).<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! <code>NWS</code> Value<br />
! Short-name<br />
! Description<br />
! Required Input Files<br />
! Optional Input Files<br />
|-<br />
| 1<br />
| Wind stress, every node, every timestep<br />
| Wind stress and atmospheric pressure are read in at all grid nodes at every model time step from the [[fort.22_file_format#NWS = 1 or 101|fort.22 file]]<br />
| [[fort.22_file_format#NWS = 1 or 101|fort.22]]<br />
|<br />
|-<br />
| 2<br />
| Wind stress, every node, every [[WTIMINC]]<br />
| Wind stress and atmospheric pressure are read in at all grid nodes at a time interval that does not equal the model time step from the [[fort.22_file_format#NWS = 2, -2, 102 or -102|fort.22 file]]. Interpolation in time is used to synchronize the wind and pressure information with the model time step. The wind time interval ([[WTIMINC]]) is specified in the [[fort.15_file_format|fort.15 file]].<br />
| [[fort.22_file_format#NWS = 2, -2, 102 or -102|fort.22]]<br />
|<br />
|-<br />
| 3<br />
| US Navy Fleet Numeric<br />
| Wind velocity is read in from a wind file from the [[fort.22_file_format#NWS = 3 or 103, Fleet Numeric Format|fort.22 file]] in US Navy Fleet Numeric format. This information is interpolated in space onto the ADCIRC grid and in time to synchronize the wind and pressure information with the model time step. Garret’s formula is used to compute wind stress from the wind velocity. Several parameters describing the Fleet Numeric wind file must be specified in the [[fort.15_file_format|fort.15 file]].<br />
| [[fort.22_file_format#NWS = 3 or 103, Fleet Numeric Format|fort.22]]<br />
|<br />
|-<br />
| 4<br />
| PBL/JAG<br />
| Wind velocity and atmospheric pressure are read in (PBL/JAG format) at selected ADCIRC grid nodes from the [[fort.22_file_format#NWS = 4, -4, 104 or -104 - PBL Hurricane Model format|fort.22]] file. Interpolation in time is used to synchronize the wind and pressure information with the model time step. Garret’s formula is used to compute wind stress from wind velocity.<br />
| [[fort.22_file_format#NWS = 4, -4, 104 or -104 - PBL Hurricane Model format|fort.22]]<br />
| <br />
|-<br />
| 5<br />
| Wind velocity, every node, every [[WTIMINC]]<br />
| Wind velocity and atmospheric pressure are read in at all grid nodes from the [[fort.22_file_format#NWS = 5, -5, 105, or -105|fort.22]] File. Interpolation in time is used to synchronize the wind and pressure information with the model time step. Garret’s formula is used to compute wind stress from wind velocity.<br />
| [[fort.22_file_format#NWS = 5, -5, 105, or -105|fort.22]]<br />
|<br />
|-<br />
| 6<br />
| wind velocity, rectangular grid, every [[WTIMINC]]<br />
| Meteorological data (U,V,P) is input on a rectangular grid (either in Longitude, Latitude or Cartesian coordinates, consistent with the grid coordinates) and interpolated in space onto the ADCIRC grid. Wind velocity (U,V @ 10 m above the water surface) must be input in units of m/s and surface atmospheric pressure (P) must be input in units of Pascals = Newtons/square meter. The meteorological grid MUST cover the entire ADCIRC mesh; that is, the ADCIRC mesh must be ENTIRELY within the meteorological grid or an error will result.<br />
| [[fort.22_file_format#NWS = 6 or 106|fort.22]]<br />
|<br />
|-<br />
| 7 {{ADC version|version=future|relation=eq}}<br />
| Wind stress, regular grid, every [[WTIMINC]]<br />
| Surface stress and pressure values are read in on a regular grid from the [[fort.22_file_format|fort.22 file]]. Currently, this feature is not supported for parallel runs because adcprep cannot decompose the files. See [https://github.com/adcirc/adcirc-cg/issues/215]. <br />
| [[fort.22_file_format#NWS = 7 or -7|fort.22]]<br />
|<br />
|-<br />
| 8<br />
| Symmetric Holland Vortex<br />
| Wind velocity and atmospheric pressure are calculated at every node on the fly by ADCIRC internally using the Dynamic Holland model.<br />
| [[Fort.22_file_format#NWS_.3D_8|fort.22]]<br />
| <br />
|-<br />
| 10<br />
| NCDC GFS<br />
| Wind velocity (10 m) and atmospheric pressure are read in from a sequence of National Weather Service (NWS) Aviation (AVN) model output files. Each AVN file is assumed to contain data on a Gaussian longitude, latitude grid at a single time.<br />
| [[fort.200_file_format|fort.200, fort.200+N, fort.200+2*N, fort.200+3*N,….,]] where N is the time interval (in hours) between successive meteorological data<br />
| <br />
|-<br />
| 11<br />
| National Weather Service Eta-29 file<br />
| Wind velocity (10 m) and atmospheric pressure are read in from a sequence of stripped down National Weather Service (NWS) ETA 29km model output files<br />
| [[fort.200_file_format|fort.200, fort.201, fort.202, fort.203,….,]]<br />
| <br />
|-<br />
|-<br />
| 12 / -12 <br />
| OWI ASCII, every [[WTIMINC]]<br />
| See [[NWS12]] for details. Wind velocities (U10, V10) and atmospheric sea level pressure (SLP) are provided in the OWI ASCII format on one to three rectangular (lat/lon) grid(s)<br />
| [[Fort.22_file_format#NWS_.3D_12|fort.22]], [[fort.221]], [[fort.222]]<br />
| [[fort.223]], [[fort.224]], [[fort.217]], [[fort.218]]<br />
|-<br />
| 13 {{ADC version|version=55|relation=ge}}<br />
| OWI NetCDF<br />
| See [[NWS13]] for details. Wind velocities (U10, V10) and atmospheric sea level pressure (SLP) fields are provided in the OWI NetCDF format as 1 or more meshgrid overlays stored in netCDF groups, supporting storm following grids on overlay 2 and on.<br />
| default is [[fort.22.nc]], see [[NWS13]]<br />
| <br />
|-<br />
| 14 / -14 {{ADC version|version=55|relation=ge}}<br />
| GRIB2/NetCDF Binary, every [[WTIMINC]]<br />
| Gridded data of wind velocities (U10, V10) and atmospheric sea level pressure (SLP) are provided in GRIB2 (e.g., GFS, CFSv2) or NetCDF (e.g., ERA5, WRF) binary files. Gridded data may be on a standard rectangular lat/lon grid, a [https://en.wikipedia.org/wiki/Gaussian_grid Gaussian grid], or a projected WRF-like grid. Requires that ADCIRC is compiled with DATETIME, NetCDF and if required, GRIB2 flags enabled (the static libraries must be compiled). Will find and read time-snaps based on the reference date, [[NCDATE]] located near or at the bottom of the [[fort.15_file_format|fort.15 file]] taking into account hot-start times etc. If the negative value is used, OWI ASCII (see NWS = 12) meteorology will overwrite the GRIB2/NetCDF meteorology data in the overlap region (except during the "skipping OWI time snap" phase). <br />
| [[fort.22x.grb2 file|fort.221.grb2, fort.222.grb2]] <br/>'''or'''<br/> [[Fort.22_file_format#NWS_.3D_14|fort.22]], [[fort.221.nc]], [[fort.222.nc]]<br />
|<br />
|-<br />
| 15<br />
| HWIND<br />
| Uses data assimilated snapshots of the wind velocity fields of tropical cyclones that are produced by the NOAA Hurricane Research Division (HRD)<br />
| [[fort.22_file_format#NWS = 15|fort.22]]<br />
| Additional HWIND files specified in the [[fort.22_file_format#NWS = 15|fort.22]] file<br />
|-<br />
| 16<br />
| GFDL<br />
| GFDL model output files produced by the Geophysical Fluid Dynamics Laboratory at NOAA. Each ASCII GFDL model output file contains one or more nested grid dataset where the nested grids are allowed to change in time. Coarse grid data is not stored where finer nest data is given.<br />
| [[fort.22_file_format#NWS = 16|fort.22]]<br />
| Additional GFDL files specified in the [[fort.22_file_format#NWS = 16|fort.22]] file<br />
|-<br />
| 19<br />
| Dynamic Asymmetric Model<ref group="note" name="nws19bad">Use of this [[Typical_ADCIRC_Parameter_Selections#Discouraged_Parameter_Selections|is discouraged]].</ref><br />
| Wind velocity and atmospheric pressure are calculated at exact finite element mesh node locations and directly coupled to ADCIRC at every time step using the asymmetric hurricane vortex formulation based on the Holland gradient wind model. The input file is assumed to correspond to the ATCF Best Track/Objective Aid/Wind Radii Format. This option uses the radii at specific wind speeds (34, 50, 64, 100 knots) reported in the four quadrants (NE, SE, SW, NW) of the storm to calculate the radius of maximum winds as a function of the azimuthal angle. Garret’s formula is used to compute wind stress from the wind velocity. This option allows the user to set a value for Rmax and Holland B Parameter. Additionally the user can select the isotachs to be used for each of the 4 quadrants. The utility program aswip_1.0.3.F located in the /wind folder will generate the NWS=19 formatted file from a NWS=9 formatted fort.22 input file.<br />
| [[fort.22_file_format#NWS = 19|fort.22]]<br />
|<br />
|-<br />
| 20<br />
| [[Generalized Asymmetric Holland Model]]<br />
| The Generalized Asymmetric Holland Model (GAHM) provides a set of theoretical and practical improvements over previous parametric meteorological vortex models in ADCIRC. The track file format is similar to that of the older Dynamic Asymmetric Model (NWS = 19) but with 8 additional columns of data.<br />
| [[fort.22_file_format#NWS = 20, Generalized Asymmetric Holland Model (GAHM)|fort.22]]<br />
|<br />
|}<br />
<br />
==Extended NWS with Ice + Waves==<br />
The following presents a summary of the extended <code>NWS</code> values to included ice-coverage and/or wind wave-coupling<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! style="font-weight:bold;" |Meteorological Data Format<br />
! style="font-weight:bold;" |Met. Only<br />
! style="font-weight:bold;" |Met. plus Waves from fort.23<br />
! style="font-weight:bold;" |Met. plus Waves SWAN<br />
! style="font-weight:bold;" |Met. plus Waves STWAVE<br />
! style="font-weight:bold;" |Met. plus Ice Coverage, Waves off<br />
! style="font-weight:bold;" |Met. plus Ice Coverage OWI-like format plus Waves from fort.23<br />
! style="font-weight:bold;" |Met. plus Ice Coverage OWI-like format plus Waves from SWAN<br />
! style="font-weight:bold;" |Met. plus Ice Coverage OWI-like format plus Waves from STWAVE<br />
|-<br />
|none<br />
| 0<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
|-<br />
|wind stress, every node, every timestep<br />
| 1<br />
| 101<br />
| 301<br />
| 401<br />
|<br />
| 12101<br />
| 12301<br />
| 12401<br />
|-<br />
|wind stress, every node, every WTIMINC<br />
| 2<br />
| 102<br />
| 302<br />
| 402<br />
|<br />
| 12102<br />
| 12302<br />
| 12402<br />
|-<br />
|US Navy Fleet Numeric<br />
| 3<br />
| 103<br />
| 303<br />
| 403<br />
|<br />
| 12103<br />
| 12303<br />
| 12403<br />
|-<br />
|PBL/JAG<br />
| 4<br />
| 104<br />
| 304<br />
| 404<br />
|<br />
| 12104<br />
| 12304<br />
| 12404<br />
|-<br />
|wind velocity, every node, every WTIMINC<br />
| 5<br />
| 105<br />
| 305<br />
| 405<br />
|<br />
| 12105<br />
| 12305<br />
| 12405<br />
|-<br />
|wind velocity, rectangular grid, every WTIMINC<br />
| 6<br />
| 106<br />
| 306<br />
| 406<br />
|<br />
| 12106<br />
| 12306<br />
| 12406<br />
|-<br />
|wind stress, regular grid, every WTIMINC<br />
| 7<br />
| 107<br />
| 307<br />
| 407<br />
|<br />
| 12107<br />
| 12307<br />
| 12407<br />
|-<br />
|symmetrc vortex model<br />
| 8<br />
| 108<br />
| 308<br />
| 408<br />
|<br />
| 12108<br />
| 12308<br />
| 12408<br />
|-<br />
|asymmetric vortex model (no longer available)<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
|-<br />
|National Weather Service AVN<br />
| 10<br />
| 110<br />
| 310<br />
| 410<br />
| 10010 (ice in 4th column of AVN file)<br />
| 12110<br />
| 12310<br />
| 12410<br />
|-<br />
|National Weather Service ETA 29km<br />
| 11<br />
| 111<br />
| 311<br />
| 411<br />
|<br />
| 12111<br />
| 12311<br />
| 12411<br />
|-<br />
|Oceanweather Inc (OWI)<br />
| 12<br />
| 112<br />
| 312<br />
| 412<br />
|<br />
| 12112<br />
| 12312<br />
| 12412<br />
|-<br />
|Oceanweather Inc (OWI) NetCDF<br />
| 13<br />
| 113?<br />
| 313?<br />
| 413?<br />
|<br />
| 12113?<br />
| 12313?<br />
| 12413?<br />
|-<br />
|GRIB2/NetCDF <br />
| 14<br />
| 114<br />
| 314<br />
| 414<br />
| 14014 (GRIB2/NetCDF format ice)<br />
| 14114 (GRIB2/NetCDF format ice)<br />
| 14314 (GRIB2/NetCDF format ice)<br />
| 14414 (GRIB2/NetCDF format ice)<br />
|-<br />
|H*Wind<br />
| 15<br />
| 115<br />
| 315<br />
| 415<br />
|<br />
| 12115<br />
| 12315<br />
| 12415<br />
|-<br />
|Dynamic Asymmetric Holland Model<ref group="note" name="nws19bad"></ref><br />
| 19<br />
| 119<br />
| 319<br />
| 419<br />
|<br />
| 12119<br />
| 12319<br />
| 12419<br />
|-<br />
|[[Generalized Asymmetric Holland Model]]<br />
| 20<br />
| 120<br />
| 320<br />
| 420<br />
|<br />
| 12120<br />
| 12320<br />
| 12420<br />
| <br />
|}<br />
<br />
== Notes ==<br />
<references group="note" /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=NWS&diff=1023NWS2020-06-26T21:46:47Z<p>Wpringle: /* Extended NWS with Ice + Waves */</p>
<hr />
<div>'''<code>NWS</code>''' is a parameter in the [[fort.15 file]] that selects the meteorological forcing input type. The value on the "[[fort.15_file_format#NWS|NWS line]]" of the fort.15 file also implicitly includes [[#Value Seen in fort.15 File|other parameters]] affecting wave coupling and ice inputs. Further, <code>NWS</code> affects not just the file type and handling of meteorological data, but also changes what the [[Fort.15_file_format#WTIMINC|meteorological parameter line]] (informally, the <code>[[WTIMINC]]</code> line) looks like in the fort.15 file. See the [[supplemental meteorological/wave/ice parameters]] page for information on the format of this line. <br />
<br />
ADCIRC supports a wide range of meteorological input formats, including moving/fixed gridded data in several file formats, tropical cyclone track and parameter data that can be turned into wind/pressure fields via one of several internal vortex models, and direct specification of wind speeds or stresses on nodes. As a result of the great flexibility and importance of this choice, several pages are devoted to the topic. In particular, see also the [[fort.22 file]] and [[wind stress]] pages. <br />
<br />
== Value Seen in fort.15 File ==<br />
In the fort.15 file, what we call "NWS" is actually a combination of several parameters. For example, given a 5-digit value on that line, <br />
-12305 ! TRICKY NWS IMPOSTER<br />
the first two digits (ten-thousands and thousands) tell us the format of ice data <code>[[NCICE]]=12</code>, the 3rd digit (hundreds) tells us the wave coupling mode <code>[[NRS]]=3</code>, and the last two digits (tens and ones) combined with the sign of the entire value tell us the meteorological forcing mode <code>NWS=-5</code>. If the value has only 3 digits then ADCIRC assumes no ice input <code>[[NCICE]]=0</code>, and if it's 2 digits then ADCIRC further assumes no wave coupling <code>[[NRS]]=0</code>. It is often presumed that when one refers to "NWS", one is referring to the 2-digit value, not what is in the fort.15 file, but be mindful of the ambiguity here. <br />
<br />
== Parameter Summary ==<br />
The following table is a summary of possible <code>NWS</code> values, their descriptions, and associated meteorological input files (required and optional).<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! <code>NWS</code> Value<br />
! Short-name<br />
! Description<br />
! Required Input Files<br />
! Optional Input Files<br />
|-<br />
| 1<br />
| Wind stress, every node, every timestep<br />
| Wind stress and atmospheric pressure are read in at all grid nodes at every model time step from the [[fort.22_file_format#NWS = 1 or 101|fort.22 file]]<br />
| [[fort.22_file_format#NWS = 1 or 101|fort.22]]<br />
|<br />
|-<br />
| 2<br />
| Wind stress, every node, every [[WTIMINC]]<br />
| Wind stress and atmospheric pressure are read in at all grid nodes at a time interval that does not equal the model time step from the [[fort.22_file_format#NWS = 2, -2, 102 or -102|fort.22 file]]. Interpolation in time is used to synchronize the wind and pressure information with the model time step. The wind time interval ([[WTIMINC]]) is specified in the [[fort.15_file_format|fort.15 file]].<br />
| [[fort.22_file_format#NWS = 2, -2, 102 or -102|fort.22]]<br />
|<br />
|-<br />
| 3<br />
| US Navy Fleet Numeric<br />
| Wind velocity is read in from a wind file from the [[fort.22_file_format#NWS = 3 or 103, Fleet Numeric Format|fort.22 file]] in US Navy Fleet Numeric format. This information is interpolated in space onto the ADCIRC grid and in time to synchronize the wind and pressure information with the model time step. Garret’s formula is used to compute wind stress from the wind velocity. Several parameters describing the Fleet Numeric wind file must be specified in the [[fort.15_file_format|fort.15 file]].<br />
| [[fort.22_file_format#NWS = 3 or 103, Fleet Numeric Format|fort.22]]<br />
|<br />
|-<br />
| 4<br />
| PBL/JAG<br />
| Wind velocity and atmospheric pressure are read in (PBL/JAG format) at selected ADCIRC grid nodes from the [[fort.22_file_format#NWS = 4, -4, 104 or -104 - PBL Hurricane Model format|fort.22]] file. Interpolation in time is used to synchronize the wind and pressure information with the model time step. Garret’s formula is used to compute wind stress from wind velocity.<br />
| [[fort.22_file_format#NWS = 4, -4, 104 or -104 - PBL Hurricane Model format|fort.22]]<br />
| <br />
|-<br />
| 5<br />
| Wind velocity, every node, every [[WTIMINC]]<br />
| Wind velocity and atmospheric pressure are read in at all grid nodes from the [[fort.22_file_format#NWS = 5, -5, 105, or -105|fort.22]] File. Interpolation in time is used to synchronize the wind and pressure information with the model time step. Garret’s formula is used to compute wind stress from wind velocity.<br />
| [[fort.22_file_format#NWS = 5, -5, 105, or -105|fort.22]]<br />
|<br />
|-<br />
| 6<br />
| wind velocity, rectangular grid, every [[WTIMINC]]<br />
| Meteorological data (U,V,P) is input on a rectangular grid (either in Longitude, Latitude or Cartesian coordinates, consistent with the grid coordinates) and interpolated in space onto the ADCIRC grid. Wind velocity (U,V @ 10 m above the water surface) must be input in units of m/s and surface atmospheric pressure (P) must be input in units of Pascals = Newtons/square meter. The meteorological grid MUST cover the entire ADCIRC mesh; that is, the ADCIRC mesh must be ENTIRELY within the meteorological grid or an error will result.<br />
| [[fort.22_file_format#NWS = 6 or 106|fort.22]]<br />
|<br />
|-<br />
| 7 {{ADC version|version=future|relation=eq}}<br />
| Wind stress, regular grid, every [[WTIMINC]]<br />
| Surface stress and pressure values are read in on a regular grid from the [[fort.22_file_format|fort.22 file]]. Currently, this feature is not supported for parallel runs because adcprep cannot decompose the files. See [https://github.com/adcirc/adcirc-cg/issues/215]. <br />
| [[fort.22_file_format#NWS = 7 or -7|fort.22]]<br />
|<br />
|-<br />
| 8<br />
| Symmetric Holland Vortex<br />
| Wind velocity and atmospheric pressure are calculated at every node on the fly by ADCIRC internally using the Dynamic Holland model.<br />
| [[Fort.22_file_format#NWS_.3D_8|fort.22]]<br />
| <br />
|-<br />
| 10<br />
| NCDC GFS<br />
| Wind velocity (10 m) and atmospheric pressure are read in from a sequence of National Weather Service (NWS) Aviation (AVN) model output files. Each AVN file is assumed to contain data on a Gaussian longitude, latitude grid at a single time.<br />
| [[fort.200_file_format|fort.200, fort.200+N, fort.200+2*N, fort.200+3*N,….,]] where N is the time interval (in hours) between successive meteorological data<br />
| <br />
|-<br />
| 11<br />
| National Weather Service Eta-29 file<br />
| Wind velocity (10 m) and atmospheric pressure are read in from a sequence of stripped down National Weather Service (NWS) ETA 29km model output files<br />
| [[fort.200_file_format|fort.200, fort.201, fort.202, fort.203,….,]]<br />
| <br />
|-<br />
|-<br />
| 12 / -12 <br />
| OWI ASCII, every [[WTIMINC]]<br />
| See [[NWS12]] for details. Wind velocities (U10, V10) and atmospheric sea level pressure (SLP) are provided in the OWI ASCII format on one to three rectangular (lat/lon) grid(s)<br />
| [[Fort.22_file_format#NWS_.3D_12|fort.22]], [[fort.221]], [[fort.222]]<br />
| [[fort.223]], [[fort.224]], [[fort.217]], [[fort.218]]<br />
|-<br />
| 13 {{ADC version|version=55|relation=ge}}<br />
| OWI NetCDF<br />
| See [[NWS13]] for details. Wind velocities (U10, V10) and atmospheric sea level pressure (SLP) fields are provided in the OWI NetCDF format as 1 or more meshgrid overlays stored in netCDF groups, supporting storm following grids on overlay 2 and on.<br />
| default is [[fort.22.nc]], see [[NWS13]]<br />
| <br />
|-<br />
| 14 / -14 {{ADC version|version=55|relation=ge}}<br />
| GRIB2/NetCDF Binary, every [[WTIMINC]]<br />
| Gridded data of wind velocities (U10, V10) and atmospheric sea level pressure (SLP) are provided in GRIB2 (e.g., GFS, CFSv2) or NetCDF (e.g., ERA5, WRF) binary files. Gridded data may be on a standard rectangular lat/lon grid, a [https://en.wikipedia.org/wiki/Gaussian_grid Gaussian grid], or a projected WRF-like grid. Requires that ADCIRC is compiled with DATETIME, NetCDF and if required, GRIB2 flags enabled (the static libraries must be compiled). Will find and read time-snaps based on the reference date, [[NCDATE]] located near or at the bottom of the [[fort.15_file_format|fort.15 file]] taking into account hot-start times etc. If the negative value is used, OWI ASCII (see NWS = 12) meteorology will overwrite the GRIB2/NetCDF meteorology data in the overlap region (except during the "skipping OWI time snap" phase). <br />
| [[fort.22x.grb2 file|fort.221.grb2, fort.222.grb2]] <br/>'''or'''<br/> [[Fort.22_file_format#NWS_.3D_14|fort.22]], [[fort.221.nc]], [[fort.222.nc]]<br />
|<br />
|-<br />
| 15<br />
| HWIND<br />
| Uses data assimilated snapshots of the wind velocity fields of tropical cyclones that are produced by the NOAA Hurricane Research Division (HRD)<br />
| [[fort.22_file_format#NWS = 15|fort.22]]<br />
| Additional HWIND files specified in the [[fort.22_file_format#NWS = 15|fort.22]] file<br />
|-<br />
| 16<br />
| GFDL<br />
| GFDL model output files produced by the Geophysical Fluid Dynamics Laboratory at NOAA. Each ASCII GFDL model output file contains one or more nested grid dataset where the nested grids are allowed to change in time. Coarse grid data is not stored where finer nest data is given.<br />
| [[fort.22_file_format#NWS = 16|fort.22]]<br />
| Additional GFDL files specified in the [[fort.22_file_format#NWS = 16|fort.22]] file<br />
|-<br />
| 19<br />
| Dynamic Asymmetric Model<ref group="note" name="nws19bad">Use of this [[Typical_ADCIRC_Parameter_Selections#Discouraged_Parameter_Selections|is discouraged]].</ref><br />
| Wind velocity and atmospheric pressure are calculated at exact finite element mesh node locations and directly coupled to ADCIRC at every time step using the asymmetric hurricane vortex formulation based on the Holland gradient wind model. The input file is assumed to correspond to the ATCF Best Track/Objective Aid/Wind Radii Format. This option uses the radii at specific wind speeds (34, 50, 64, 100 knots) reported in the four quadrants (NE, SE, SW, NW) of the storm to calculate the radius of maximum winds as a function of the azimuthal angle. Garret’s formula is used to compute wind stress from the wind velocity. This option allows the user to set a value for Rmax and Holland B Parameter. Additionally the user can select the isotachs to be used for each of the 4 quadrants. The utility program aswip_1.0.3.F located in the /wind folder will generate the NWS=19 formatted file from a NWS=9 formatted fort.22 input file.<br />
| [[fort.22_file_format#NWS = 19|fort.22]]<br />
|<br />
|-<br />
| 20<br />
| [[Generalized Asymmetric Holland Model]]<br />
| The Generalized Asymmetric Holland Model (GAHM) provides a set of theoretical and practical improvements over previous parametric meteorological vortex models in ADCIRC. The track file format is similar to that of the older Dynamic Asymmetric Model (NWS = 19) but with 8 additional columns of data.<br />
| [[fort.22_file_format#NWS = 20, Generalized Asymmetric Holland Model (GAHM)|fort.22]]<br />
|<br />
|}<br />
<br />
==Extended NWS with Ice + Waves==<br />
The following presents a summary of the extended <code>NWS</code> values to included ice-coverage and/or wind wave-coupling<br />
<br />
{| class="wikitable" border="1"<br />
|-<br />
! style="font-weight:bold;" |Meteorological Data Format<br />
! style="font-weight:bold;" |Met. Only<br />
! style="font-weight:bold;" |Met. plus Waves from fort.23<br />
! style="font-weight:bold;" |Met. plus Waves SWAN<br />
! style="font-weight:bold;" |Met. plus Waves STWAVE<br />
! style="font-weight:bold;" |Met. plus Ice Coverage, Waves off<br />
! style="font-weight:bold;" |Met. plus Ice Coverage OWI-like format plus Waves from fort.23<br />
! style="font-weight:bold;" |Met. plus Ice Coverage OWI-like format plus Waves from SWAN<br />
! style="font-weight:bold;" |Met. plus Ice Coverage OWI-like format plus Waves from STWAVE<br />
|-<br />
|none<br />
| 0<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
|-<br />
|wind stress, every node, every timestep<br />
| 1<br />
| 101<br />
| 301<br />
| 401<br />
|<br />
| 12101<br />
| 12301<br />
| 12401<br />
|-<br />
|wind stress, every node, every WTIMINC<br />
| 2<br />
| 102<br />
| 302<br />
| 402<br />
|<br />
| 12102<br />
| 12302<br />
| 12402<br />
|-<br />
|US Navy Fleet Numeric<br />
| 3<br />
| 103<br />
| 303<br />
| 403<br />
|<br />
| 12103<br />
| 12303<br />
| 12403<br />
|-<br />
|PBL/JAG<br />
| 4<br />
| 104<br />
| 304<br />
| 404<br />
|<br />
| 12104<br />
| 12304<br />
| 12404<br />
|-<br />
|wind velocity, every node, every WTIMINC<br />
| 5<br />
| 105<br />
| 305<br />
| 405<br />
|<br />
| 12105<br />
| 12305<br />
| 12405<br />
|-<br />
|wind velocity, rectangular grid, every WTIMINC<br />
| 6<br />
| 106<br />
| 306<br />
| 406<br />
|<br />
| 12106<br />
| 12306<br />
| 12406<br />
|-<br />
|wind stress, regular grid, every WTIMINC<br />
| 7<br />
| 107<br />
| 307<br />
| 407<br />
|<br />
| 12107<br />
| 12307<br />
| 12407<br />
|-<br />
|symmetrc vortex model<br />
| 8<br />
| 108<br />
| 308<br />
| 408<br />
|<br />
| 12108<br />
| 12308<br />
| 12408<br />
|-<br />
|asymmetric vortex model (no longer available)<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
| n/a<br />
|-<br />
|National Weather Service AVN<br />
| 10<br />
| 110<br />
| 310<br />
| 410<br />
| 14010 (GRIB2/NetCDF format ice)<br />
| 12110<br />
| 12310<br />
| 12410<br />
|-<br />
|National Weather Service ETA 29km<br />
| 11<br />
| 111<br />
| 311<br />
| 411<br />
|<br />
| 12111<br />
| 12311<br />
| 12411<br />
|-<br />
|Oceanweather Inc (OWI)<br />
| 12<br />
| 112<br />
| 312<br />
| 412<br />
|<br />
| 12112<br />
| 12312<br />
| 12412<br />
|-<br />
|Oceanweather Inc (OWI) NetCDF<br />
| 13<br />
| 113?<br />
| 313?<br />
| 413?<br />
|<br />
| 12113?<br />
| 12313?<br />
| 12413?<br />
|-<br />
|GRIB2/NetCDF <br />
| 14<br />
| 114<br />
| 314<br />
| 414<br />
| 14014 (GRIB2/NetCDF format ice)<br />
| 14114 (GRIB2/NetCDF format ice)<br />
| 14314 (GRIB2/NetCDF format ice)<br />
| 14414 (GRIB2/NetCDF format ice)<br />
|-<br />
|H*Wind<br />
| 15<br />
| 115<br />
| 315<br />
| 415<br />
|<br />
| 12115<br />
| 12315<br />
| 12415<br />
|-<br />
|Dynamic Asymmetric Holland Model<ref group="note" name="nws19bad"></ref><br />
| 19<br />
| 119<br />
| 319<br />
| 419<br />
|<br />
| 12119<br />
| 12319<br />
| 12419<br />
|-<br />
|[[Generalized Asymmetric Holland Model]]<br />
| 20<br />
| 120<br />
| 320<br />
| 420<br />
|<br />
| 12120<br />
| 12320<br />
| 12420<br />
| <br />
|}<br />
<br />
== Notes ==<br />
<references group="note" /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=File:Channel_Elev.gif&diff=884File:Channel Elev.gif2020-06-09T02:10:50Z<p>Wpringle: Wpringle uploaded a new version of File:Channel Elev.gif</p>
<hr />
<div>== Summary ==<br />
Elevation time series for the idealized channel problem</div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=883Idealized Channel Problem2020-06-09T02:08:50Z<p>Wpringle: </p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Note that the short 6 hour length of the test is chosen only to limit simulation time for the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite] where the test case been found. Users may extend the simulation length to simulate more of the inundating phase of the incoming wave. <br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1000px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br />
[[File:Channel_Elev.gif|500px|thumb|Elevation time series for the idealized channel problem]] [[File:Channel_Vel.gif|500px|thumb|North-south velocity time series for the idealized channel problem]]<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files. [[Grid_Development_and_Editing#OceanMesh2D|OceanMesh2D]] functions can be used to automatically generate the sponge_generator_layer attribute ([https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/utilities/Calc_Sponge.m Calc_Sponge]) and the input files ([https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/utilities/Make_f5354.m Make_f5354]).<br />
*[[Fort.14_file_format|IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=882Idealized Channel Problem2020-06-09T02:07:51Z<p>Wpringle: /* Options/Features Tested */</p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Note that the short 6 hour length of the test is chosen only to limit simulation time for the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite] where the test case been found. Users may extend the simulation length to simulate the inundating phase of the incoming wave. <br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1000px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br />
[[File:Channel_Elev.gif|500px|thumb|Elevation time series for the idealized channel problem]] [[File:Channel_Vel.gif|500px|thumb|North-south velocity time series for the idealized channel problem]]<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files. [[Grid_Development_and_Editing#OceanMesh2D|OceanMesh2D]] functions can be used to automatically generate the sponge_generator_layer attribute ([https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/utilities/Calc_Sponge.m Calc_Sponge]) and the input files ([https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/utilities/Make_f5354.m Make_f5354]).<br />
*[[Fort.14_file_format|IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=881Idealized Channel Problem2020-06-09T02:06:07Z<p>Wpringle: /* Options/Features Tested */</p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Note that the short 6 hour length of the test is chosen only to limit simulation time for the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite] where the test case been found. Users may extend the simulation length to simulate the inundating phase of the incoming wave. <br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1000px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br />
[[File:Channel_Elev.gif|500px|thumb|Elevation time series for the idealized channel problem]] [[File:Channel_Vel.gif|500px|thumb|North-south velocity time series for the idealized channel problem]]<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files. [[OceanMesh2D]] functions exist for generating the sponge_generator_layer attribute ([https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/utilities/Calc_Sponge.m Calc_Sponge]), and the input files ([https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/utilities/Make_f5354.m Make_f5354]).<br />
*[[Fort.14_file_format|IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=880Idealized Channel Problem2020-06-09T02:01:58Z<p>Wpringle: </p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Note that the short 6 hour length of the test is chosen only to limit simulation time for the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite] where the test case been found. Users may extend the simulation length to simulate the inundating phase of the incoming wave. <br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1000px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br />
[[File:Channel_Elev.gif|500px|thumb|Elevation time series for the idealized channel problem]] [[File:Channel_Vel.gif|500px|thumb|North-south velocity time series for the idealized channel problem]]<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
*[[Fort.14_file_format|IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=879Idealized Channel Problem2020-06-09T02:01:35Z<p>Wpringle: </p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Note that the short 6 hour length of the test is chosen to only limit simulation time for the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite] where the test case been found. Users may extend the simulation length to simulate the inundating phase of the incoming wave. <br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1000px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br />
[[File:Channel_Elev.gif|500px|thumb|Elevation time series for the idealized channel problem]] [[File:Channel_Vel.gif|500px|thumb|North-south velocity time series for the idealized channel problem]]<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
*[[Fort.14_file_format|IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=878Idealized Channel Problem2020-06-09T01:58:02Z<p>Wpringle: </p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1000px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br />
[[File:Channel_Elev.gif|500px|thumb|Elevation time series for the idealized channel problem]] [[File:Channel_Vel.gif|500px|thumb|North-south velocity time series for the idealized channel problem]]<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
*[[Fort.14_file_format|IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=File:Channel_Vel.gif&diff=877File:Channel Vel.gif2020-06-09T01:51:54Z<p>Wpringle: North-south velocity time series for the idealized channel problem</p>
<hr />
<div>== Summary ==<br />
North-south velocity time series for the idealized channel problem</div>Wpringlehttps://wiki.adcirc.org/index.php?title=File:Channel_Elev.gif&diff=876File:Channel Elev.gif2020-06-09T01:46:40Z<p>Wpringle: Elevation time series for the idealized channel problem</p>
<hr />
<div>== Summary ==<br />
Elevation time series for the idealized channel problem</div>Wpringlehttps://wiki.adcirc.org/index.php?title=Alaskan_Winter_Storm_with_Ice&diff=875Alaskan Winter Storm with Ice2020-06-08T23:43:52Z<p>Wpringle: </p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of the storm tides in a regional Alaska domain under astronomical and atmospheric forcing in November 2011 during a strong winter storm in the presence of sea ice (affecting the surface wind drag)<ref name=Brian>Joyce, B.R., Pringle, W.J., Wirasaet, D., Westerink, J.J., Van der Westhuysen, A.J., Grumbine, R., Feyen, J., 2019. High resolution modeling of western Alaskan tides and storm surge under varying sea ice conditions. Ocean Model. 141, 101421. doi:10.1016/j.ocemod.2019.101421</ref>. The results of interest are the global elevations, velocities and meteorology. The test finishes in about 5 minutes in serial ADCIRC for two weeks of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_alaska_ice-2d GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh was generated using the OceanMesh2D Alaska [https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/Examples/Example_8_AK.m Example_8_AK.m]. The domain encompasses the Gulf of Alaska, Bering Sea, and Chukchi Sea with a minimum resolution of 5 km, comprised of 15,876 vertices and 27,757 triangular elements.<br />
<br />
== Options/Features Tested ==<br />
*<code>[[ICS]]</code> = 20: Equal-Area cylindrical projection. <br />
*<code>[[IM]]</code> = 513111: Uses the implicit scheme for the linear component of the gravity wave term (computational time step is 4 minutes). <br />
*<code>[[NTIP]]</code> = 2: Equilibrium tide + self-attraction and loading tide (read from a [[fort.24 file]]) forcing for 8 tidal constituents.<br />
*<code>[[NWS]]</code> = 14014: Reads from GRIB2 files that specify the global atmospheric forcing and sea-ice concentration (6-hourly CFSv2 reanalysis data). Sea-ice concentration affects the wind drag coefficient<ref name=Brian></ref>. <br />
*<code>[[WTIMINC]]</code> = 21600, 21600: First value gives the temporal interval of the GRIB2 met data (6 hours), second value gives the temporal interval of the GRIB2 ice data (6 hours) - these should always be the same.<br />
*<code>[[A00, B00, C00]]</code> = 0.4, 0.4, 0.2: Ensures that the implicit scheme is stable with a fairly large time step.<br />
*<code>[[ESLM]]</code> = -0.2: Enables the Smagorinsky turbulence closure with a coefficient of 0.2.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[fort.13_file#Internal_Tide_Energy_Conversion|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.<br />
*[[Fort.15_file_format#Namelists|&WarnElevControl namelist]]: Set "WarnElev", the warning elevation level, to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).<br />
*[[Fort.15_file_format#Namelists|&metControl namelist]]: Set "rhoAir", to 1.29193 (density of air at 0 deg C for 1013 mbar); set "WindDragLimit" equal to 0.0025; set "invertedBarometerOnElevationBoundary" to true (in Alaska extremely large-scale low pressure systems persist and cross over the open boundaries, so it is important to have the inverted barometer condition along the elevation specified boundary); set "outputWindDrag" to true.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Fort.14_file_format&diff=874Fort.14 file format2020-06-08T23:38:55Z<p>Wpringle: </p>
<hr />
<div>The basic file structure is shown below. Each line of input data is represented by a line containing the input variable name(s). Blank lines are only to enhance readability. Loops indicate multiple lines of input. <br />
<br />
[[AGRID]]<br />
<br />
[[NE]], [[NP]]<br />
<br />
for k=1 to [[NP]]<br />
: [[JN]], [[X(JN)]], [[Y(JN)]], [[DP(JN)]]<br />
end k loop<br />
<br />
for k=1 to [[NE]]<br />
: [[JE]], [[NHY]], [[NM(JE,1),NM(JE,2), NM(JE,3)]]<br />
end k loop<br />
<br />
[[NOPE]]<br />
<br />
[[NETA]]<br />
<br />
for k=1 to [[NOPE]]<br />
: [[NVDLL(k)]], [[IBTYPEE(k)]]<br />
: for j=1 to [[NVDLL(k)]]<br />
:: [[NBDV(k,j)]]<br />
:end j loop<br />
end k loop<br />
<br />
[[NBOU]]<br />
<br />
[[NVEL]]<br />
<br />
for k=1 to [[NBOU]]<br />
: [[NVELL(k)]], IBTYPE(k)<br />
: for j=1,NVELL(k)<br />
:: [[NBVV(k,j)]] include this line only if [[IBTYPE(k)]] = 0, 1, 2, 10, 11, 12, 20, 21, 22, 30<br />
:: [[NBVV(k,j)]], [[IBCONN(k,j)]] include this line only if [[IBTYPE(k)]] = 94<br />
:: [[NBVV(k,j)]], [[BARLANHT(k,j)]], [[BARLANCFSP(k,j)]] include this line only if IBTYPE(k) = 3, 13, 23<br />
:: [[NBVV(k,j)]], [[IBCONN(k,j)]], [[BARINHT(k,j)]], [[BARINCFSB(k,j)]], [[BARINCFSP(k,j)]] include this line only if IBTYPE(k) = 4, 24<br />
:: [[NBVV(k,j)]], [[IBCONN(k,j)]], [[BARINHT(k,j)]], [[BARINCFSB(k,j)]], [[BARINCFSP(k,j)]], [[PIPEHT(k,j)]], [[PIPECOEF(k,j)]], [[PIPEDIAM(k,j)]], include this line only if IBTYPE(k) = 5, 25<br />
:end j loop<br />
end k loop</div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=873Idealized Channel Problem2020-06-08T23:34:10Z<p>Wpringle: /* Options/Features Tested */</p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1750px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
*[[Fort.14_file_format|IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=872Idealized Channel Problem2020-06-08T23:33:45Z<p>Wpringle: /* Options/Features Tested */</p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1750px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
*[[Fort.14_file_format IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=871Idealized Channel Problem2020-06-08T23:32:14Z<p>Wpringle: /* Options/Features Tested */</p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1750px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
[[IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=870Idealized Channel Problem2020-06-08T23:28:12Z<p>Wpringle: </p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1750px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br /><br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=869Idealized Channel Problem2020-06-08T23:26:39Z<p>Wpringle: /* Mesh */</p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1750px|thumb|Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. Center: Mesh topo-bathy. Right: The sponge strength coefficients.]]<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=868Idealized Channel Problem2020-06-08T23:23:40Z<p>Wpringle: /* Mesh */</p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
[[File:IdealChannel.png|1500px]]<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
<br />
== References ==<br />
<references /></div>Wpringlehttps://wiki.adcirc.org/index.php?title=File:IdealChannel.png&diff=867File:IdealChannel.png2020-06-08T23:20:56Z<p>Wpringle: Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations.
Center: Mesh topo-bathy.
Right: The sponge streng...</p>
<hr />
<div>== Summary ==<br />
Left: Mesh triangulation and resolution. Blue line shows the elevation specified boundary condition location, green and yellow lines on the sides show the periodic lateral boundary condition locations. <br />
Center: Mesh topo-bathy.<br />
Right: The sponge strength coefficients.</div>Wpringlehttps://wiki.adcirc.org/index.php?title=Idealized_Channel_Problem&diff=866Idealized Channel Problem2020-06-08T23:08:03Z<p>Wpringle: /* Options/Features Tested */</p>
<hr />
<div>This example tests ADCIRC version 55 (and beyond). It tests the simulation of a diurnal tide on a sloping beach with a channel along its centerline (adapted from<ref name=Keith>Roberts, K.J., Dietrich, J.C., Wirasaet, D., Pringle, W.J., Westerink, J.J., 2020. Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding. In Preparation, pp.37.</ref>). It tests lateral periodic boundary conditions and the absorption-generation sponge layer<ref name=Pringle>Pringle, W.J., Wirasaet, D., Suhardjo, A., Meixner, J., Westerink, J.J., Kennedy, A.B., Nong, S., 2018. Finite-Element Barotropic Model for the Indian and Western Pacific Oceans: Tidal Model-Data Comparisons and Sensitivities. Ocean Model. 129, 13–38. doi:10.1016/j.ocemod.2018.07.003</ref><ref name=Pringle2>Pringle, W.J., Gonzalez-lopez, J., Joyce, B., Westerink, J.J., van der Westhuysen, A.J., 2019. Baroclinic Coupling Improves Depth-Integrated Modeling of Coastal Sea Level Variations around Puerto Rico and the U.S. Virgin Islands. J. Geophys. Res. Ocean. 124, 2196–2217. doi:10.1029/2018JC014682</ref>. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_ideal_channel-2d-parallel GitHub test suite].<br />
<br />
== Mesh == <br />
The mesh is comprised of 64,415 vertices and 127,784 triangular elements, with resolution in the 10-60 m range. The mesh is symmetrical in the east-west direction so that the east and west lateral boundary vertices match for the application of the periodic lateral boundary conditions. An elevation specified boundary condition and absorption-generation sponge layer is prescribed at the southern end of the domain.<br />
<br />
== Options/Features Tested ==<br />
*<code>[[IM]]</code> = 111112: Uses the explicit scheme (computational time step is 2 seconds). <br />
*<code>[[A00, B00, C00]]</code> = 0.0, 1.0, 0.0: Must be used with explicit scheme.<br />
*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]]. <br />
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]]. <br />
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity). <br />
*[[Fort.13_file#Absorption-generation_Sponge_Layer|sponge_generator_layer]]: Applies a sponge layer to absorb outgoing waves while generating incoming waves. In this case incoming diurnal tidal waves are generated using the [[fort.53001]] and [[fort.54001]] input files.<br />
<br />
== References ==<br />
<references /></div>Wpringle