Perfectly matched layer absorbing boundary condition for simulation of tsunami wave propagation

Corresponding Author: Takuto Maeda (maeda (at) eri.u-tokyo.ac.jp)

This software simulates tsunami wave propagation under the linear long wave (LLW) and the liner dispersive wave (LDW) approximation with the perfectly matched layer absorbing boundary condition. The theory and example are described in the accompanying paper:

Maeda, T., H. Tsushima, and T. Furumura, An effective absorbing boundary condition for linear long-wave and linear dispersive wave tsunami simulations, Earth, Planets, and Space, 68, 63, doi:10.1186/s40623-016-0436-y, 2016.

The authors request that the user to cite the above accompanying paper in any publications that result from the use of this software, although this is not an obligation.

MIT

Fortran2008

The author confirms it successfully work at the following compilers: gfortran 5.1.0 (Mac OSX) Intel fortran 2015 (Linux)

The LLW and LDW models are given in separated programs (llw.f90 and ldw.f90).

They include the following three block files when they are compiled:

• blk_param.f90: parameter for model size, grid width and the title for output
• blk_initheight.f90: define the initial sea height
• blk_bathymetry.f90: define the bathymetry
• blk_station.f90: define the station information for output

After modification of these files, compile llw.f90 and/or ldw.f90. Block files (blk_*.f90) are automatically included into the main program files. It is not necessary to compile them separately.