A numerical algorithm for computing tsunami wave amplitude

A numerical algorithm for computing tsunami wave front amplitude is proposed. The first step consists in solving an appropriate eikonal equation. The eikonal equation is solved by the Godunov approach and the bicharacteristic method. The qualitative comparison of the two above methods is described. Then a change in variables associated with the eikonal solution is introduced. At the last step, using the expansion of the fundamental solution of shallow water equations in the sum of singular and regular parts, we obtain the Cauchy problem for the wave amplitude. This approach allows one to reduce computer costs. The numerical results are presented.