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://kingkong.amath.washington.edu/clawpack/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