[cig-commits] r18954 - short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters

willic3 at geodynamics.org willic3 at geodynamics.org
Tue Sep 20 21:13:06 PDT 2011


Author: willic3
Date: 2011-09-20 21:13:05 -0700 (Tue, 20 Sep 2011)
New Revision: 18954

Added:
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock40.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/finalslip_rupture_lock40.spatialdb
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8_lock20_tau0-1.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-1.spatialdb
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliprate_creep_lock40.spatialdb
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/tet4.cfg
Removed:
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_graded_20km.cfg
Modified:
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock20.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-2.spatialdb
Log:
More changes to S&P benchmark.



Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock20.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock20.cfg	2011-09-21 01:00:48 UTC (rev 18953)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock20.cfg	2011-09-21 04:13:05 UTC (rev 18954)
@@ -32,59 +32,59 @@
 
 # one
 [pylithapp.timedependent.interfaces.fault.eq_srcs.one.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # two
 [pylithapp.timedependent.interfaces.fault.eq_srcs.two.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # three
 [pylithapp.timedependent.interfaces.fault.eq_srcs.three.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # four
 [pylithapp.timedependent.interfaces.fault.eq_srcs.four.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # five
 [pylithapp.timedependent.interfaces.fault.eq_srcs.five.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # six
 [pylithapp.timedependent.interfaces.fault.eq_srcs.six.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # seven
 [pylithapp.timedependent.interfaces.fault.eq_srcs.seven.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # eight
 [pylithapp.timedependent.interfaces.fault.eq_srcs.eight.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # nine
 [pylithapp.timedependent.interfaces.fault.eq_srcs.nine.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # ten
 [pylithapp.timedependent.interfaces.fault.eq_srcs.ten.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # creep
 [pylithapp.timedependent.interfaces.fault.eq_srcs.creep]
 slip_function = pylith.faults.ConstRateSlipFn
 origin_time = 0.0*s
-slip_function.slip_rate.iohandler.filename = sliprate_creep.spatialdb
+slip_function.slip_rate.iohandler.filename = sliprate_creep_lock20.spatialdb
 slip_function.slip_time.iohandler.filename = sliptime.spatialdb
 
 # ----------------------------------------------------------------------

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock40.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock40.cfg	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock40.cfg	2011-09-21 04:13:05 UTC (rev 18954)
@@ -0,0 +1,116 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to the time-dependent
+# interface conditions. The general parameters are specificed in the
+# pylithapp.cfg file which PyLith reads by default, and the
+# problem-specific parameters are specified in the spbm_xxxx.cfg files.
+#
+# To run the simulation:
+# pylith spbm_xxxx.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+output.vertex_info_fields = [normal_dir,strike_dir,dip_dir,final_slip_creep,final_slip_one,final_slip_two,final_slip_three,final_slip_four,final_slip_five,final_slip_six,final_slip_seven,final_slip_eight,final_slip_nine,final_slip_ten,slip_time_creep,slip_time_one,slip_time_two,slip_time_three,slip_time_four,slip_time_five,slip_time_six,slip_time_seven,slip_time_eight,slip_time_nine,slip_time_ten]
+
+# Set earthquake ruptures
+eq_srcs = [creep,one,two,three,four,five,six,seven,eight,nine,ten]
+eq_srcs.creep.origin_time = 00.0*year
+eq_srcs.one.origin_time = 200.0*year
+eq_srcs.two.origin_time = 400.0*year
+eq_srcs.three.origin_time = 600.0*year
+eq_srcs.four.origin_time = 800.0*year
+eq_srcs.five.origin_time = 1000.0*year
+eq_srcs.six.origin_time = 1200.0*year
+eq_srcs.seven.origin_time = 1400.0*year
+eq_srcs.eight.origin_time = 1600.0*year
+eq_srcs.nine.origin_time = 1800.0*year
+eq_srcs.ten.origin_time = 2000.0*year
+
+# one
+[pylithapp.timedependent.interfaces.fault.eq_srcs.one.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# two
+[pylithapp.timedependent.interfaces.fault.eq_srcs.two.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# three
+[pylithapp.timedependent.interfaces.fault.eq_srcs.three.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# four
+[pylithapp.timedependent.interfaces.fault.eq_srcs.four.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# five
+[pylithapp.timedependent.interfaces.fault.eq_srcs.five.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# six
+[pylithapp.timedependent.interfaces.fault.eq_srcs.six.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# seven
+[pylithapp.timedependent.interfaces.fault.eq_srcs.seven.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# eight
+[pylithapp.timedependent.interfaces.fault.eq_srcs.eight.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# nine
+[pylithapp.timedependent.interfaces.fault.eq_srcs.nine.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# ten
+[pylithapp.timedependent.interfaces.fault.eq_srcs.ten.slip_function]
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# creep
+[pylithapp.timedependent.interfaces.fault.eq_srcs.creep]
+slip_function = pylith.faults.ConstRateSlipFn
+origin_time = 0.0*s
+slip_function.slip_rate.iohandler.filename = sliprate_creep_lock40.spatialdb
+slip_function.slip_time.iohandler.filename = sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Specify that output for the entire domain occurs every 50 years.
+[pylithapp.problem.formulation.output.domain]
+output_freq = time_step
+time_step = 50*year
+
+# Specify that output for the ground surface occurs every 10 years.
+[pylithapp.problem.formulation.output.subdomain]
+output_freq = time_step
+time_step = 10*year
+
+# Specify that output for the fault occurs every 10 years.
+[pylithapp.timedependent.interfaces.fault.output]
+output_freq = time_step
+time_step = 10*year
+
+# Specify that output for both materials occurs every 50 years.
+[pylithapp.timedependent.materials.elastic.output]
+output_freq = time_step
+time_step = 50*year
+
+[pylithapp.timedependent.materials.maxwell.output]
+output_freq = time_step
+time_step = 50*year
+

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/finalslip_rupture_lock40.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/finalslip_rupture_lock40.spatialdb	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/finalslip_rupture_lock40.spatialdb	2011-09-21 04:13:05 UTC (rev 18954)
@@ -0,0 +1,32 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of slip on the
+// fault surface. We apply 4.0 m of left-lateral slip, corresponding to
+// the accumulated slip deficit assuming a slip rate of 2 cm/yr over
+// 200 years. This slip is applied from the surface down to z > -40 km.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+  num-values = 3
+  value-names =  left-lateral-slip  reverse-slip  fault-opening
+  value-units =  m  m  m
+  num-locs = 5
+  data-dim = 1 // Data is specified along a line.
+  space-dim = 3
+  cs-data = cartesian {
+    to-meters = 1.0e+3 // Specify coordinates in km for convenience.
+    space-dim = 3
+  } // cs-data
+} // SimpleDB
+// Columns are
+// (1) x coordinate (km)
+// (2) y coordinate (km)
+// (3) z coordinate (km)
+// (4) left-lateral-slip (m) (right-lateral is negative)
+// (5) reverse-slip (m)
+// (6) fault-opening (m)
+0.0  0.0   0.00     4.0  0.0  0.0
+0.0  0.0  -39.00    4.0  0.0  0.0
+0.0  0.0  -40.00    0.0  0.0  0.0
+0.0  0.0  -41.00    0.0  0.0  0.0
+0.0  0.0  -1000.0   0.0  0.0  0.0

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8.cfg	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8.cfg	2011-09-21 04:13:05 UTC (rev 18954)
@@ -0,0 +1,34 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to the mesh used.
+# The general parameters are specificed in the pylithapp.cfg file which
+# PyLith reads by default. Information controlling time-dependent fault
+# behavior is contained in the file timedep.cfg.
+#
+# To run the simulation:
+# pylith spbm_hex8_graded_20km.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+# Set filename of mesh to import.
+filename = ../meshes/hex8_refined_20km.exo
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+elastic.quadrature.cell = pylith.feassemble.FIATLagrange
+elastic.quadrature.cell.dimension = 3
+maxwell.quadrature.cell = pylith.feassemble.FIATLagrange
+maxwell.quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2

Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8_lock20_tau0-1.cfg (from rev 18953, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_graded_20km.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8_lock20_tau0-1.cfg	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8_lock20_tau0-1.cfg	2011-09-21 04:13:05 UTC (rev 18954)
@@ -0,0 +1,54 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to the mesh used.
+# The general parameters are specificed in the pylithapp.cfg file which
+# PyLith reads by default. Information controlling time-dependent fault
+# behavior is contained in the file timedep.cfg.
+#
+# To run the simulation:
+# pylith spbm_hex8_graded_20km.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+# Set filename of mesh to import.
+filename = ../meshes/spbm_hex8_graded_20km.exo
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+elastic.quadrature.cell = pylith.feassemble.FIATLagrange
+elastic.quadrature.cell.dimension = 3
+maxwell.quadrature.cell = pylith.feassemble.FIATLagrange
+maxwell.quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK files.
+[pylithapp.problem.formulation.output.domain.writer]
+filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km.vtk
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-groundsurf.vtk
+
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-fault.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-elastic.vtk
+
+[pylithapp.timedependent.materials.maxwell.output]
+writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-maxwell.vtk
+

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-1.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-1.spatialdb	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-1.spatialdb	2011-09-21 04:13:05 UTC (rev 18954)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties. In this case, the material properties are uniform.
+// The properties here correspond to lambda = mu = 30 GPa, and a viscosity
+// that provides a relaxation time (2 * vis/mu) of 200 years. This
+// corresponds to a value of tau0 (recurrence_time/relaxation_time) of 1.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+  num-values = 4 // number of material property values
+  value-names =  density vs vp viscosity // names of the material property values
+  value-units =  kg/m**3  m/s  m/s Pa*s
+  num-locs = 1 // number of locations
+  data-dim = 0
+  space-dim = 3
+  cs-data = cartesian {
+    to-meters = 1.0
+    space-dim = 3
+  }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m**3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0  0.0  0.0   2700.0  3333.333333333333  5773.502691896258 9.46728e19

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-2.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-2.spatialdb	2011-09-21 01:00:48 UTC (rev 18953)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-2.spatialdb	2011-09-21 04:13:05 UTC (rev 18954)
@@ -2,6 +2,9 @@
 //
 // This spatial database specifies the distribution of material
 // properties. In this case, the material properties are uniform.
+// The properties here correspond to lambda = mu = 30 GPa, and a viscosity
+// that provides a relaxation time (2 * vis/mu) of 100 years. This
+// corresponds to a value of tau0 (recurrence_time/relaxation_time) of 2.
 //
 #SPATIAL.ascii 1
 SimpleDB {

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliprate_creep_lock40.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliprate_creep_lock40.spatialdb	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliprate_creep_lock40.spatialdb	2011-09-21 04:13:05 UTC (rev 18954)
@@ -0,0 +1,33 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of slip rate on
+// the fault surface for a constant slip rate time function. In this
+// case we specify a uniform slip rate of -2.0 cm/year for z = -40.0 km and
+// zero slip for z > -40.0 km. There is no reverse-slip or fault
+// opening.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+  num-values = 3
+  value-names =  left-lateral-slip  reverse-slip  fault-opening
+  value-units =  cm/year  cm/year  cm/year
+  num-locs = 5
+  data-dim = 1 // Values are specified along a line.
+  space-dim = 3
+  cs-data = cartesian {
+    to-meters = 1.0e+3 // Specify coordinates in km for convenience.
+    space-dim = 3
+  } // cs-data
+} // SimpleDB
+// Columns are
+// (1) x coordinate (km)
+// (2) y coordinate (km)
+// (3) z coordinate (km)
+// (4) left-lateral-slip rate (m/s) (right-lateral is negative)
+// (5) reverse-slip rate (m/s)
+// (6) fault-opening rate (m/s)
+0.0   0.0     0.00     0.00  0.0  0.0
+0.0   0.0   -39.00     0.00  0.0  0.0
+0.0   0.0   -40.00     2.00  0.0  0.0
+0.0   0.0   -41.00     0.00  0.0  0.0
+0.0   0.0   -1000.0    0.00  0.0  0.0

Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_graded_20km.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_graded_20km.cfg	2011-09-21 01:00:48 UTC (rev 18953)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_graded_20km.cfg	2011-09-21 04:13:05 UTC (rev 18954)
@@ -1,54 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies parameters specific to the mesh used.
-# The general parameters are specificed in the pylithapp.cfg file which
-# PyLith reads by default. Information controlling time-dependent fault
-# behavior is contained in the file timedep.cfg.
-#
-# To run the simulation:
-# pylith spbm_hex8_graded_20km.cfg timedep.cfg
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.reader]
-# Set filename of mesh to import.
-filename = ../meshes/spbm_hex8_graded_20km.exo
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials]
-elastic.quadrature.cell = pylith.feassemble.FIATLagrange
-elastic.quadrature.cell.dimension = 3
-maxwell.quadrature.cell = pylith.feassemble.FIATLagrange
-maxwell.quadrature.cell.dimension = 3
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK files.
-[pylithapp.problem.formulation.output.domain.writer]
-filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km.vtk
-
-[pylithapp.problem.formulation.output.subdomain]
-writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-groundsurf.vtk
-
-[pylithapp.timedependent.interfaces.fault.output.writer]
-filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-fault.vtk
-
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-elastic.vtk
-
-[pylithapp.timedependent.materials.maxwell.output]
-writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-maxwell.vtk
-

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/tet4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/tet4.cfg	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/tet4.cfg	2011-09-21 04:13:05 UTC (rev 18954)
@@ -0,0 +1,34 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to the mesh used.
+# The general parameters are specificed in the pylithapp.cfg file which
+# PyLith reads by default. Information controlling time-dependent fault
+# behavior is contained in the file timedep.cfg.
+#
+# To run the simulation:
+# pylith spbm_hex8_graded_20km.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+# Set filename of mesh to import.
+filename = ../meshes/tet4_graded_4km.exo
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+elastic.quadrature.cell = pylith.feassemble.FIATSimplex
+elastic.quadrature.cell.shape = tetrahedron
+maxwell.quadrature.cell = pylith.feassemble.FIATSimplex
+maxwell.quadrature.cell.shape = tetrahedron
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle



More information about the CIG-COMMITS mailing list