README.html

README.html
Source file: README.txt
Directory: /home/rjl/git/rjleveque/clawpack-4.x/apps/tsunami/bowl-slosh
Converted: Tue Jul 26 2011 at 12:58:49 using clawcode2html
This documentation file will not reflect any later changes in the source file.

$\phantom{******** If you see this on the webpage then the browser could not locate *********}$
$\phantom{******** the jsMath file load.js *********}$

$\newcommand{\vector}[1]{\left[\begin{array}{c} #1 \end{array}\right]}$ $\newenvironment{matrix}{\left[\begin{array}{cccccccccc}} {\end{array}\right]}$ $\newcommand{\A}{{\cal A}}$ $\newcommand{\W}{{\cal W}}$

GeoClaw Sample Code

Waves in a parabolic bowl with a flat surface sloshing around. An exact analytic solution is known in which the surface stays flat.

In this code, $x$ and $y$ are in meters (iccordsys=1 in setrun.py [.html]).

Topography: $B(x,y) = h_0((x^2 + y^2)/a^2 -1)$,

Depth: $h(x,y,t) = \max\left(0,~~ (\sigma h_0/a^2)(2x\cos(\omega t) + 2y\sin(\omega t) - \sigma) - B(x,y)\right)$

Velocities: $u(x,y,t) = -\sigma \omega \sin(\omega t),\qquad v(x,y,t) = \sigma \omega \cos(\omega t).$

where $\omega = \sqrt{2gh_0} / a$.

The period of oscillation is $T = 2\pi / \omega$.

The following parameters are currently hardwired several places:

$a = 1, ~~\sigma = 0.5, ~~h = 0.1,~~g = 9.81$

This should be cleaned up: better to put them in a setprob.data file that is read in where needed.

References

W. C. Thacker, Some exact solutions to the nonlinear shallow water wave equations, J. Fluid Mech. 107 (1981), 499-508.
J.M. Gallardo, C. Pares, and M. Castro, On a well-balanced high-order finite volume scheme for shallow water equations with topography and dry areas, J. Comput. Phys. 227(2007) 574-601.
Y. Xing, X. Zhang and C.-W. Shu, Positivity preserving high order well balanced discontinuous Galerkin methods for the shallow water equations , submitted to Advances in Water Resources. preprint

This test problem has been used in several other papers too.

To run the code, see Instructions

Plots of results

After running this code and creating plots via "make .plots", you should be able to view the plots in _plots/_PlotIndex.html .

Fortran files

Makefile [.html]

Determines which version of fortran files are used when compiling the code with make and specifies where output and plots should be directed. Type "make .help" at the Unix prompt for options.

Python files

maketopo.py [.html]: Used to create topo file and qinit data file.
setrun.py [.html]: This file contains a function that specifies what run-time parameters will be used.
setplot.py [.html]: This file contains a function that specifies what plots will be done and sets various plotting parameters.

Data files

The .data files are automatically generated using the information in setrun.py [.html].

Instructions

To make topo and qinit data files:

  $ make topo

To make all data files, edit setrun.py and then

  $ make .data

To run code:

You may need to type

  $ make new

to make sure the modules are accessible.

Then run the code with

  $ make .output

To plot results, either generate html pages via:

  $ make .plots

or view interactively using ipython and Iplotclaw.

All of this can be done with:

  $ source make_all.py