This one dimensional test problem consists of a flat ocean floor, linear continental slope, flat continental shelf, and a solid wall reflecting boundary.
It is designed to illustrate how a tsunami wave is modified as it moves from the deep ocean onto the continental shelf, and the manner in which some of the energy can be trapped on the shelf and bounce back and forth.
Note: For more about shoaling of tsunami waves on continental shelves, and the manner in which the width of the continental slope affects the transmission and reflection, see the recent paper:
See also: http://faculty.washington.edu/rjl/pubs/Shoaling2019/index.html
Additional examples from this paper, and Jupyter notebooks, can be found in the GitHub repository https://github.com/rjleveque/shoaling_paper_figures.
Try the following bathymetry by changing lines in the file setrun.py and the rerunning the code and plotting the results via:
make .plots
Original: a step discontinuity between deep ocean and shallow shelf:
Bocean = -4000. Bshelf = -200. width = 1. start = -30.e3
With a shallower shelf, note that the wave moves slower but is more amplified:
Bocean = -4000. Bshelf = -50. width = 1. start = -30.e3
With a wide continental shelf rather than a discontinuity, note that there is less energy trapped on the shelf:
Bocean = -4000. Bshelf = -200. width = 100.e3 start = -130.e3
The IPython notebook Shelf_1d.ipynb illustrates these cases.