GeoClaw Sample Code

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

See $CLAW/doc/sphinx/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

Note: This fault model has be updated since this example was originally constructed, and the GeoClaw source code has been improved to use a combination of 180 subfaults specified at earthquake.usgs.gov/earthquakes/eqinthenews/2010/us2010tfan/finite_fault.php in place of the single fault plane used in the example included in Clawpack 4.5.1 and before.

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

To create all these directories of plots, use

  $ source make_all.sh

Plots created:

• _plots/_PlotIndex.html 3 level run
• _plots_noamr/_PlotIndex.html single level with same resolution
• _plots_4levels/_PlotIndex.html 4 levels with higher resolution

The script compare_gauges.py [.html] can be used plot a comparison of the gauge data for all three runs with the DART Buoy 32412 data, resulting in a plot dart_gauges.png.

The script plot_dart.py [.html] is a version that plots on top of computational results obtained by NOAA, 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.

subfaults.py [.html]

Used to create dtopo file from the file subfault.txt [.html], obtained from earthquake.usgs.gov/earthquakes/eqinthenews/2010/us2010tfan/finite_fault.php .

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

To run the 3 different cases for the figures in the paper,

  $ source make_all.sh

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