The goal of the simulation is to study the generation and open-ocean propagation of the tsunami waves triggered by the Tōhoku earthquake, with particular focus on their evolution over realistic bathymetry.

The underlying mathematical model is the two-dimensional nonlinear shallow-water equations, formulated in terms of the water depth \( h \), the depth-averaged momenta \( hu \) and \( hv \), and the sea-floor topography \( z \). The computational domain is a rectangular ocean region of size \( 4 \times 10^6 m \) by \( 4 \times 10^6 m \), discretised on an adaptively refined Cartesian grid. The bathymetry before the earthquake and the coseismic seafloor displacement are read from NetCDF data sets.