begin_html [use: doc/doc.css] [use:jsMath]

GeoClaw Sample Code

Waves in a parabolic bowl generated by a Gaussian hump of water at the origin. The solution should remain circularly symmetric. In this code, $x$ and $y$ are in meters (iccordsys=1 in [code:setrun.py]). Refinement is allowed only in the upper right quadrant, using the regions specified in [code:setrun.py], see [http://www.clawpack.org/users/regions.html]. Gauges are placed near the shore along the x-axis and along the diagonal. Ideally these gauges would show identical results. To run the code, see [link: #instructions Instructions]

Plots of results

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

Fortran files

[code: Makefile]
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.
[code: setprob.f]
Standard for for GeoClaw

Python files

[code: maketopo.py]
Used to create topo file and qinit data file.
[code: setrun.py]
This file contains a function that specifies what run-time parameters will be used.
[code: setplot.py]
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 [code: setrun.py]. [name: instructions]

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 }}} end_html