--- Welcome to the official ADCIRCWiki site! The site is currently under construction, with limited information. ---

For general information, see the ADCIRC site, ADCIRC FAQ, or the ADCIRC Wikipedia page. For model documentation not yet available on the wiki, see the ADCIRC site. New content is being continuously added to the ADCIRCWiki, and material from the main ADCIRC site will be gradually transitioned over to the wiki.

Idealized Channel Problem

Revision as of 18:52, 8 June 2020 by Wpringle (talk | contribs) (Created page with "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<ref name=Keith>Roberts, K.J...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

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[1]. The test finishes in about 8 minutes in parallel ADCIRC (2 processors) for 6 hours of simulation. Find the test at the GitHub test suite.


The mesh was generated using the OceanMesh2D Alaska 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.

Options/Features Tested

  • ICS = 20: Equal-Area cylindrical projection.
  • IM = 513111: Uses the implicit scheme for the linear component of the gravity wave term (computational time step is 4 minutes).
  • NTIP = 2: Equilibrium tide + self-attraction and loading tide (read from a fort.24 file) forcing for 8 tidal constituents.
  • NWS = 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[2].
  • WTIMINC = 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.
  • A00, B00, C00 = 0.4, 0.4, 0.2: Ensures that the implicit scheme is stable with a fairly large time step.
  • ESLM = -0.2: Enables the Smagorinsky turbulence closure with a coefficient of 0.2.
  • NOUTGE = 5: Outputs the global elevations into a netCDF4 fort.63 file.
  • NOUTGV = 5: Outputs the global velocities into a netCDF4 fort.64 file.
  • NOUTGM = 5: Outputs the global meteorology into a netCDF4 fort.73 file (pressure) and a netCDF4 fort.74 file (velocity).
  • 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.
  • &WarnElevControl namelist: Set "WarnElev", the warning elevation level, to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).
  • &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.


  1. 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.
  2. Cite error: Invalid <ref> tag; no text was provided for refs named Brian