--- 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: Difference between revisions

From ADCIRCWiki
Jump to navigation Jump to search
No edit summary
 
(14 intermediate revisions by the same user not shown)
Line 1: Line 1:
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>). 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].
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.


== Mesh ==  
== Mesh ==  
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 generating-absorbing sponge layer is prescribed at the southern end of the domain.
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.
 
[[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.]]
 
[[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]]


== Options/Features Tested ==
== Options/Features Tested ==
Line 10: Line 14:
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]].  
*<code>[[NOUTGV]]</code> = 5: Outputs the global velocities into a netCDF4 [[fort.64 file]].  
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity).  
*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity).  
*[[fort.13_file#Sponge|sponge_generator_layer]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.
*[[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]).
*[[Fort.14_file_format|IBTYPE=94]]: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.


== References ==
== References ==
<references />
<references />

Latest revision as of 23:01, 7 April 2021

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[1]). It tests lateral periodic boundary conditions and the absorption-generation sponge layer[2][3]. 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 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.

Mesh

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.

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.
Elevation time series for the idealized channel problem
North-south velocity time series for the idealized channel problem

Options/Features Tested

  • IM = 111112: Uses the explicit scheme (computational time step is 2 seconds).
  • A00, B00, C00 = 0.0, 1.0, 0.0: Must be used with explicit scheme.
  • 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).
  • 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 can be used to automatically generate the sponge_generator_layer attribute (Calc_Sponge) and the input files (Make_f5354).
  • IBTYPE=94: Node pairs are matched along opposite lateral boundaries where a periodic (repeating) boundary condition is applied.

References

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