GeoClaw Sample Code

Sample earthquake-generated tsunami based on the 2010 Chile event. Primarily to illustrate how set up topo and dtopo files.

See www.clawpack.org/users/topo.html for more information.

See maketopo.py [.html]. This downloads topography data from the GeoClaw topography database that originally came from the NOAA National Geophysical Data Center (NGDC) using Design-a-grid .

USGS site for earthquake data: earthquake.usgs.gov/earthquakes/eqinthenews/2010/us2010tfan

One gauge is included at the location of DART Buoy 32412: www.ndbc.noaa.gov/station_page.php?station=32412

To download other DART buoy data, see www.ndbc.noaa.gov/dart.shtml

See also the NOAA/PMEL Center for Tsunami Research site for this tsunami, containing DART data, tide gauges, MOST simulations: nctr.pmel.noaa.gov/chile20100227/

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 .

The script plot_dart.py [.html] can be used plot a comparison of the gauge data with the DART Buoy 32412 data, resulting in a plot dart.png.

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.

setprob.f [.html]

Standard for for GeoClaw

Python files

maketopo.py [.html]

Used to create topo file and dtopo 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.

setplot_imshow.py [.html]

Version of setplot illustrating how the '2d_imshow' plot type can be used in place of '2d_pcolor'. This plots much faster.

plot_dart.py [.html]

This file contains a script that can be used plot a comparison of the gauge data with the DART Buoy 32412 data.

Data files

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

Instructions

To make topo and dtopo data files:

  $ make topo

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

  $ make .data

To run code:

  $ make .output

To plot results, either generate html pages via:

  $ make .plots

All of this can be done with:

  $ source make_all.py

For more documentation, see www.clawpack.org/users .