This uses a very simple source model with a single fault plane and constant slip, from an early inversion by the USGS. The input parameters are specified in usgs100227.cfg and converted into the dtopo file usgs100227.txt by the script maketopo.py. This can be run via:
make topo
which also downloads a topo file for the ocean bathymetry. This bathymetry originally came from the NOAA National Geophysical Data Center (NGDC) using Design-a-grid .
A single gauge captures the sea surface elevation at the location of DART buoy 32412.
setplot_speeds.py is a version of setplot that also plots velocities.
Creating kml files to view on Google Earth
setplot_kml.py is a version of the setplot file that produces kml files suitable for viewing on Google Earth or other platforms supporting kml.
To use:
make .output make plots SETPLOT_FILE=setplot_kml.py
and then open _output/Chile_2010.kmz from Google Earth. Note that this is a zipped directory containing a number of kml files and opening this way should open them all. You can then select which ones to view (e.g. the simultion, the grid patches, regions specified in setrun.py, and/or gauges).
You can unzip the kmz file if you want to extract individual kml files.
See http://www.clawpack.org/googleearth_plotting.html for documentation.