--- 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.

Difference between revisions of "Alaskan Winter Storm with Ice"

From ADCIRCWiki
Jump to: navigation, search
(Mesh)
(Mesh)
Line 2: Line 2:
  
 
== Mesh ==  
 
== Mesh ==  
The mesh is a coarse representation of Alaska (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 generated using the OceanMesh2D  Alaska [https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/Examples/Example_8_AK.m Example_8_AK.m].
+
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.
  
 
== Options/Features Tested ==
 
== Options/Features Tested ==

Revision as of 22:31, 8 June 2020

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 (affecting the surface wind drag)[1]. 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 GitHub test suite.

Mesh

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 = -22: Uses the Mercator projection with a coordinate rotation to remove the pole singularity (need to provide a fort.rotm).
  • IM = 513113: Uses the fully implicit scheme for the gravity wave term (computational time step is 12 minutes).
  • NTIP = 2: Equilibrium tide + self-attraction and loading tide (read from a fort.24 file) forcing for 10 tidal constituents.
  • NWS = -14: Reads from GRIB2 files that specify the global atmospheric forcing (6-hourly CFS reanalysis data) in addition to OWI ASCII files that specify the 3-hourly atmospheric forcing in the Hurricane Katrina landfall region.
  • WTIMINC = 21600, 10800: First value gives the temporal interval of the GRIB2 met data (6 hours), second value gives the temporal interval of the OWI met data (3 hours) in seconds.
  • A00, B00, C00 = 0.5, 0.5, 0: Ensures that the fully implicit scheme is stable with a 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.
  • quadratic_friction_coefficient: Spatially varying quadratic bottom friction fort.13 file attribute.

References

  1. 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