[cig-commits] r14197 - in short/2.5D/benchmarks/savageprescott: . results

willic3 at geodynamics.org willic3 at geodynamics.org
Mon Mar 2 18:21:14 PST 2009

Author: willic3
Date: 2009-03-02 18:21:14 -0800 (Mon, 02 Mar 2009)
New Revision: 14197

Added top-level README (more will be added in subdirectories), and
directories to hold PyLith results.

Added: short/2.5D/benchmarks/savageprescott/README
--- short/2.5D/benchmarks/savageprescott/README	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/README	2009-03-03 02:21:14 UTC (rev 14197)
@@ -0,0 +1,36 @@
+This directory contains the files necessary to run a simulation of the
+Savage and Prescott [1978] solution for a strike-slip fault in an elastic
+layer overlying a viscoelastic half-space. In this model, coseismic slip
+occurs at regular recurrence intervals on the upper locked part of the
+fault, while constant-velocity creep occurs on the lower part of the fault.
+This directory contains the files necessary to run the problem using
+PyLith, as well as a utils directory that contains a Python script and
+.cfg file to compute the analytical solution.
+To run the PyLith simulation, go to the parameters directory and run one
+of the shell scripts, or just run PyLith directly using the appropriate
+.cfg files on the command line (see the README in the parameters directory).
+Note that if you run using one of the graded meshes, the problem will take
+quite a while to run. Once the problem has finished running, there will be a
+number of .vtk files in the results directory, which may be viewed with a
+package such as ParaView.
+The simplest way to compare the numerical and analytical solutions is to
+look at profile results using a plotting package. The savpres_ss.py
+script will produce a comma-separated-values file (savageprescott.csv) that
+is easy to use with a plotting package. For results from PyLith, it is
+probably easiest to generate profile results using ParaView (e.g., the
+"Plot Over Line" filter). By loading the entire time series for the
+ground surface results, it is possible to generate profile results for
+all time steps. For comparison with the analytical results, it is easiest
+to use a profile coincident with the x-axis that extends from x=0 to the
+right-hand side of the mesh. The profiles can also be saved as CSV files
+from ParaView.
+Note that it is not possible to directly compare the PyLith results with the
+analytical solution. This is because the analytical solution assumes steady
+asthenospheric flow at the beginning of computations, while the numerical
+solution requires several earthquake cycles to reach this state. For direct
+comparison, it is therefore more reasonable to reference all results to
+the displacement field immediately following an earthquake. This is what was
+done for the figures contained in the "figures" directory.

More information about the CIG-COMMITS mailing list