[cig-commits] r14192 - short/2.5D/benchmarks/savageprescott/parameters

willic3 at geodynamics.org willic3 at geodynamics.org
Mon Mar 2 16:40:52 PST 2009


Author: willic3
Date: 2009-03-02 16:40:52 -0800 (Mon, 02 Mar 2009)
New Revision: 14192

Added:
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded1_20km.cfg
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded1_20km.sh
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded2_20km.cfg
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded2_20km.sh
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif1_20km.cfg
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif1_20km.sh
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif2_20km.cfg
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif2_20km.sh
Removed:
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.cfg
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.sh
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.cfg
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.sh
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.cfg
   short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.sh
Modified:
   short/2.5D/benchmarks/savageprescott/parameters/finalslip_rupture.spatialdb
   short/2.5D/benchmarks/savageprescott/parameters/mat_maxwell.spatialdb
   short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg
   short/2.5D/benchmarks/savageprescott/parameters/sliprate_creep.spatialdb
   short/2.5D/benchmarks/savageprescott/parameters/timedep.cfg
Log:
Deleted obsolete .cfg and script files.
Added new .cfg and script files.
Fixed spatialdb files for new fault BC (slip only in elastic layer),
and fixed mistake in Maxwell spatialdb.



Modified: short/2.5D/benchmarks/savageprescott/parameters/finalslip_rupture.spatialdb
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/finalslip_rupture.spatialdb	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/finalslip_rupture.spatialdb	2009-03-03 00:40:52 UTC (rev 14192)
@@ -1,9 +1,9 @@
 // -*- 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 right-lateral slip, corresponding to
+// 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 = -16 km.
+// 200 years. This slip is applied from the surface down to z = -20 km.
 //
 #SPATIAL.ascii 1
 SimpleDB {
@@ -25,7 +25,7 @@
 // (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  -13.50   -4.0  0.0  0.0
-0.0  0.0  -13.51    0.0  0.0  0.0
-0.0  0.0  -486.00   0.0  0.0  0.0
+0.0  0.0   0.00     4.0  0.0  0.0
+0.0  0.0  -20.00    4.0  0.0  0.0
+0.0  0.0  -20.01    0.0  0.0  0.0
+0.0  0.0  -1000.0   0.0  0.0  0.0

Modified: short/2.5D/benchmarks/savageprescott/parameters/mat_maxwell.spatialdb
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/mat_maxwell.spatialdb	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/mat_maxwell.spatialdb	2009-03-03 00:40:52 UTC (rev 14192)
@@ -5,7 +5,7 @@
 //
 #SPATIAL.ascii 1
 SimpleDB {
-  num-values = 3 // number of material property values
+  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 // units (conversions not supported)
   num-locs = 1 // number of locations

Modified: short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg	2009-03-03 00:40:52 UTC (rev 14192)
@@ -21,9 +21,9 @@
 #fiatlagrange = 1
 
 [pylithapp.journal.debug]
-pylithapp = 1
-problem = 1
-implicit = 1
+# pylithapp = 1
+# problem = 1
+# implicit = 1
 
 # ----------------------------------------------------------------------
 # mesh_generator
@@ -47,8 +47,9 @@
 # ----------------------------------------------------------------------
 [pylithapp.timedependent]
 # Set materials to an array of 2 materials:
-#   'elastic_upper' and 'elastic_lower'
+#   'elastic' and 'maxwell'
 materials = [elastic,maxwell]
+materials.maxwell = pylith.materials.MaxwellIsotropic3D
 
 [pylithapp.timedependent.materials.elastic]
 label = Upper elastic layer
@@ -65,6 +66,8 @@
 quadrature = pylith.feassemble.quadrature.Quadrature3D
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
+output.cell_data_fields = [total_strain,viscous_strain,stress]
+output.cell_info_fields = [mu,lambda,density,maxwell_time]
 
 # ----------------------------------------------------------------------
 # PETSc
@@ -77,7 +80,7 @@
 ksp_monitor = true
 ksp_view = true
 log_summary = true
-ksp_max_it = 100
+ksp_max_it = 800
 ksp_gmres_restart = 70
 ksp_monitor_singular_value = true
 # start_in_debugger = true

Modified: short/2.5D/benchmarks/savageprescott/parameters/sliprate_creep.spatialdb
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/sliprate_creep.spatialdb	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/sliprate_creep.spatialdb	2009-03-03 00:40:52 UTC (rev 14192)
@@ -11,7 +11,7 @@
   num-values = 3
   value-names =  left-lateral-slip  reverse-slip  fault-opening
   value-units =  m/s  m/s  m/s
-  num-locs = 4
+  num-locs = 6
   data-dim = 1 // Values are specified along a line.
   space-dim = 3
   cs-data = cartesian {
@@ -27,6 +27,8 @@
 // (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   -13.50     0.00  0.0  0.0
-0.0   0.0   -13.51    -6.337617562806e-10  0.0  0.0
-0.0   0.0  -486.00    -6.337617562806e-10  0.0  0.0
+0.0   0.0   -20.00     0.00  0.0  0.0
+0.0   0.0   -20.01     6.337617562806e-10  0.0  0.0
+0.0   0.0   -40.00     6.337617562806e-10  0.0  0.0
+0.0   0.0   -40.01     0.00  0.0  0.0
+0.0   0.0   -1000.0    0.00  0.0  0.0

Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded1_20km.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded1_20km.cfg	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded1_20km.cfg	2009-03-03 00:40:52 UTC (rev 14192)
@@ -0,0 +1,127 @@
+# -*- 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_graded1_20km.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+# Change the mesh importer to the CUBIT importer.
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+# Set filename of mesh to import.
+filename = ../meshes/spbm_hex8_graded1_20km.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
+bc = [x_pos,x_neg,z_neg]
+
+# Set interfaces to an array of 1 fault: 'fault'.
+interfaces = [fault]
+
+[pylithapp.timedependent.implicit]
+# Set the output to an array of 2 output managers.
+# We will output the solution over the domain and the ground surface.
+output = [domain,subdomain]
+
+# Set subdomain component to OutputSolnSubset (subset of domain).
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+
+# We change spatial database for the rate of change from FixedDOFDB
+# (which has zero velocity) to UniformDB (which allows nonzero
+# velocities).
+
+# Set the parameters for the desired boundary conditions.
+#
+# The label corresponds to the nodeset ID in CUBIT.
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0, 1]
+label = 11
+db.label = Dirichlet BC on +x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on +x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of 1 cm/year.
+rate_db.data = [0.0,3.1688087814e-10]
+
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1]
+label = 12
+db.label = Dirichlet BC on -x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on -x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of -1 cm/year.
+rate_db.data = [0.0,-3.1688087814e-10]
+
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [2]
+label = 16
+db.label = Dirichlet BC on -z
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+
+# Set the parameters for the fault interface conditions.
+
+[pylithapp.timedependent.interfaces.fault]
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+mat_db.iohandler.filename = mat_elastic.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain.writer]
+filename = ../results/spbm_hex8_graded1_20km/spbm_hex8_graded1_20km.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK domain output of solution over ground surface.
+[pylithapp.problem.formulation.output.subdomain]
+label = 17 ; nodeset for subdomain
+writer.filename = ../results/spbm_hex8_graded1_20km/spbm_hex8_graded1_20km-groundsurf.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for vtk fault rupture output.
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = ../results/spbm_hex8_graded1_20km/spbm_hex8_graded1_20km-fault.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_graded1_20km/spbm_hex8_graded1_20km-statevars-elastic.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.maxwell.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_graded1_20km/spbm_hex8_graded1_20km-statevars-maxwell.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+

Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded1_20km.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded1_20km.sh	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded1_20km.sh	2009-03-03 00:40:52 UTC (rev 14192)
@@ -0,0 +1 @@
+pylith spbm_hex8_graded1_20km.cfg timedep.cfg


Property changes on: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded1_20km.sh
___________________________________________________________________
Name: svn:executable
   + *

Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded2_20km.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded2_20km.cfg	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded2_20km.cfg	2009-03-03 00:40:52 UTC (rev 14192)
@@ -0,0 +1,127 @@
+# -*- 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_graded2_20km.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+# Change the mesh importer to the CUBIT importer.
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+# Set filename of mesh to import.
+filename = ../meshes/spbm_hex8_graded2_20km.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
+bc = [x_pos,x_neg,z_neg]
+
+# Set interfaces to an array of 1 fault: 'fault'.
+interfaces = [fault]
+
+[pylithapp.timedependent.implicit]
+# Set the output to an array of 2 output managers.
+# We will output the solution over the domain and the ground surface.
+output = [domain,subdomain]
+
+# Set subdomain component to OutputSolnSubset (subset of domain).
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+
+# We change spatial database for the rate of change from FixedDOFDB
+# (which has zero velocity) to UniformDB (which allows nonzero
+# velocities).
+
+# Set the parameters for the desired boundary conditions.
+#
+# The label corresponds to the nodeset ID in CUBIT.
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0, 1]
+label = 11
+db.label = Dirichlet BC on +x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on +x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of 1 cm/year.
+rate_db.data = [0.0,3.1688087814e-10]
+
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1]
+label = 12
+db.label = Dirichlet BC on -x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on -x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of -1 cm/year.
+rate_db.data = [0.0,-3.1688087814e-10]
+
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [2]
+label = 16
+db.label = Dirichlet BC on -z
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+
+# Set the parameters for the fault interface conditions.
+
+[pylithapp.timedependent.interfaces.fault]
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+mat_db.iohandler.filename = mat_elastic.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain.writer]
+filename = ../results/spbm_hex8_graded2_20km/spbm_hex8_graded2_20km.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK domain output of solution over ground surface.
+[pylithapp.problem.formulation.output.subdomain]
+label = 17 ; nodeset for subdomain
+writer.filename = ../results/spbm_hex8_graded2_20km/spbm_hex8_graded2_20km-groundsurf.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for vtk fault rupture output.
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = ../results/spbm_hex8_graded2_20km/spbm_hex8_graded2_20km-fault.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_graded2_20km/spbm_hex8_graded2_20km-statevars-elastic.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.maxwell.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_graded2_20km/spbm_hex8_graded2_20km-statevars-maxwell.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+

Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded2_20km.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded2_20km.sh	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded2_20km.sh	2009-03-03 00:40:52 UTC (rev 14192)
@@ -0,0 +1 @@
+pylith spbm_hex8_graded2_20km.cfg timedep.cfg


Property changes on: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded2_20km.sh
___________________________________________________________________
Name: svn:executable
   + *

Deleted: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.cfg	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.cfg	2009-03-03 00:40:52 UTC (rev 14192)
@@ -1,127 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies parameters specific to the boundary and
-# interface conditions. The general parameters are specificed in the
-# pylithapp.cfg file which PyLith reads by default. Information controlling
-# time stepping is contained in the file timedep.cfg.
-#
-# To run the simulation:
-# pylith spbm_hex8_unif_8000m.cfg timedep.cfg
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-# Change the mesh importer to the CUBIT importer.
-importer = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.importer]
-# Set filename of mesh to import.
-filename = ../meshes/spbm_hex8_graded_1000m.exo
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
-bc = [x_pos,x_neg,z_neg]
-
-# Set interfaces to an array of 1 fault: 'fault'.
-interfaces = [fault]
-
-[pylithapp.timedependent.implicit]
-# Set the output to an array of 2 output managers.
-# We will output the solution over the domain and the ground surface.
-output = [domain,subdomain]
-
-# Set subdomain component to OutputSolnSubset (subset of domain).
-output.subdomain = pylith.meshio.OutputSolnSubset
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-
-# We change spatial database for the rate of change from FixedDOFDB
-# (which has zero velocity) to UniformDB (which allows nonzero
-# velocities).
-
-# Set the parameters for the desired boundary conditions.
-#
-# The label corresponds to the nodeset ID in CUBIT.
-
-[pylithapp.timedependent.bc.x_pos]
-fixed_dof = [0, 1]
-label = 11
-db.label = Dirichlet BC on +x
-rate_db = spatialdata.spatialdb.UniformDB
-rate_db.label = Dirichlet rate BC on +x
-rate_db.values = [dof-0,dof-1]
-# This rate (in m/s) gives a value of -1 cm/year.
-rate_db.data = [0.0,-3.1688087814e-10]
-
-[pylithapp.timedependent.bc.x_neg]
-fixed_dof = [0, 1]
-label = 12
-db.label = Dirichlet BC on -x
-rate_db = spatialdata.spatialdb.UniformDB
-rate_db.label = Dirichlet rate BC on -x
-rate_db.values = [dof-0,dof-1]
-# This rate (in m/s) gives a value of 1 cm/year.
-rate_db.data = [0.0,3.1688087814e-10]
-
-[pylithapp.timedependent.bc.z_neg]
-fixed_dof = [2]
-label = 16
-db.label = Dirichlet BC on -z
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces]
-
-# Set the parameters for the fault interface conditions.
-
-[pylithapp.timedependent.interfaces.fault]
-label = 10
-quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-mat_db.iohandler.filename = mat_elastic.spatialdb
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.domain.writer]
-filename = ../results/spbm_hex8_graded_1000m.vtk
-time_format = %04.0f
-time_constant = 1.0*year
-
-# Give basename for VTK domain output of solution over ground surface.
-[pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
-writer.filename = ../results/spbm_hex8_graded_1000m-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Give basename for vtk fault rupture output.
-[pylithapp.timedependent.interfaces.fault.output.writer]
-filename = ../results/spbm_hex8_graded_1000m-fault.vtk
-time_format = %04.0f
-time_constant = 1.0*year
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = ../results/spbm_hex8_graded_1000m-statevars-elastic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-[pylithapp.timedependent.materials.maxwell.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = ../results/spbm_hex8_graded_1000m-statevars-maxwell.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-

Deleted: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.sh	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.sh	2009-03-03 00:40:52 UTC (rev 14192)
@@ -1 +0,0 @@
-pylith spbm_hex8_graded_1000m.cfg timedep.cfg

Deleted: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.cfg	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.cfg	2009-03-03 00:40:52 UTC (rev 14192)
@@ -1,127 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies parameters specific to the boundary and
-# interface conditions. The general parameters are specificed in the
-# pylithapp.cfg file which PyLith reads by default. Information controlling
-# time stepping is contained in the file timedep.cfg.
-#
-# To run the simulation:
-# pylith spbm_hex8_unif_8000m.cfg timedep.cfg
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-# Change the mesh importer to the CUBIT importer.
-importer = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.importer]
-# Set filename of mesh to import.
-filename = ../meshes/spbm_hex8_graded_3000m.exo
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
-bc = [x_pos,x_neg,z_neg]
-
-# Set interfaces to an array of 1 fault: 'fault'.
-interfaces = [fault]
-
-[pylithapp.timedependent.implicit]
-# Set the output to an array of 2 output managers.
-# We will output the solution over the domain and the ground surface.
-output = [domain,subdomain]
-
-# Set subdomain component to OutputSolnSubset (subset of domain).
-output.subdomain = pylith.meshio.OutputSolnSubset
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-
-# We change spatial database for the rate of change from FixedDOFDB
-# (which has zero velocity) to UniformDB (which allows nonzero
-# velocities).
-
-# Set the parameters for the desired boundary conditions.
-#
-# The label corresponds to the nodeset ID in CUBIT.
-
-[pylithapp.timedependent.bc.x_pos]
-fixed_dof = [0, 1]
-label = 11
-db.label = Dirichlet BC on +x
-rate_db = spatialdata.spatialdb.UniformDB
-rate_db.label = Dirichlet rate BC on +x
-rate_db.values = [dof-0,dof-1]
-# This rate (in m/s) gives a value of -1 cm/year.
-rate_db.data = [0.0,-3.1688087814e-10]
-
-[pylithapp.timedependent.bc.x_neg]
-fixed_dof = [0, 1]
-label = 12
-db.label = Dirichlet BC on -x
-rate_db = spatialdata.spatialdb.UniformDB
-rate_db.label = Dirichlet rate BC on -x
-rate_db.values = [dof-0,dof-1]
-# This rate (in m/s) gives a value of 1 cm/year.
-rate_db.data = [0.0,3.1688087814e-10]
-
-[pylithapp.timedependent.bc.z_neg]
-fixed_dof = [2]
-label = 16
-db.label = Dirichlet BC on -z
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces]
-
-# Set the parameters for the fault interface conditions.
-
-[pylithapp.timedependent.interfaces.fault]
-label = 10
-quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-mat_db.iohandler.filename = mat_elastic.spatialdb
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.domain.writer]
-filename = ../results/spbm_hex8_graded_3000m.vtk
-time_format = %04.0f
-time_constant = 1.0*year
-
-# Give basename for VTK domain output of solution over ground surface.
-[pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
-writer.filename = ../results/spbm_hex8_graded_3000m-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Give basename for vtk fault rupture output.
-[pylithapp.timedependent.interfaces.fault.output.writer]
-filename = ../results/spbm_hex8_graded_3000m-fault.vtk
-time_format = %04.0f
-time_constant = 1.0*year
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = ../results/spbm_hex8_graded_3000m-statevars-elastic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-[pylithapp.timedependent.materials.maxwell.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = ../results/spbm_hex8_graded_3000m-statevars-maxwell.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-

Deleted: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.sh	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.sh	2009-03-03 00:40:52 UTC (rev 14192)
@@ -1 +0,0 @@
-pylith spbm_hex8_graded_3000m.cfg timedep.cfg

Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif1_20km.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif1_20km.cfg	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif1_20km.cfg	2009-03-03 00:40:52 UTC (rev 14192)
@@ -0,0 +1,127 @@
+# -*- 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_unif1_20km.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+# Change the mesh importer to the CUBIT importer.
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+# Set filename of mesh to import.
+filename = ../meshes/spbm_hex8_unif1_20km.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
+bc = [x_pos,x_neg,z_neg]
+
+# Set interfaces to an array of 1 fault: 'fault'.
+interfaces = [fault]
+
+[pylithapp.timedependent.implicit]
+# Set the output to an array of 2 output managers.
+# We will output the solution over the domain and the ground surface.
+output = [domain,subdomain]
+
+# Set subdomain component to OutputSolnSubset (subset of domain).
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+
+# We change spatial database for the rate of change from FixedDOFDB
+# (which has zero velocity) to UniformDB (which allows nonzero
+# velocities).
+
+# Set the parameters for the desired boundary conditions.
+#
+# The label corresponds to the nodeset ID in CUBIT.
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0, 1]
+label = 11
+db.label = Dirichlet BC on +x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on +x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of 1 cm/year.
+rate_db.data = [0.0,3.1688087814e-10]
+
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1]
+label = 12
+db.label = Dirichlet BC on -x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on -x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of -1 cm/year.
+rate_db.data = [0.0,-3.1688087814e-10]
+
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [2]
+label = 16
+db.label = Dirichlet BC on -z
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+
+# Set the parameters for the fault interface conditions.
+
+[pylithapp.timedependent.interfaces.fault]
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+mat_db.iohandler.filename = mat_elastic.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain.writer]
+filename = ../results/spbm_hex8_unif1_20km/spbm_hex8_unif1_20km.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK domain output of solution over ground surface.
+[pylithapp.problem.formulation.output.subdomain]
+label = 17 ; nodeset for subdomain
+writer.filename = ../results/spbm_hex8_unif1_20km/spbm_hex8_unif1_20km-groundsurf.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for vtk fault rupture output.
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = ../results/spbm_hex8_unif1_20km/spbm_hex8_unif1_20km-fault.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_unif1_20km/spbm_hex8_unif1_20km-statevars-elastic.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.maxwell.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_unif1_20km/spbm_hex8_unif1_20km-statevars-maxwell.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+

Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif1_20km.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif1_20km.sh	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif1_20km.sh	2009-03-03 00:40:52 UTC (rev 14192)
@@ -0,0 +1 @@
+pylith spbm_hex8_unif1_20km.cfg timedep.cfg


Property changes on: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif1_20km.sh
___________________________________________________________________
Name: svn:executable
   + *

Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif2_20km.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif2_20km.cfg	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif2_20km.cfg	2009-03-03 00:40:52 UTC (rev 14192)
@@ -0,0 +1,127 @@
+# -*- 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_unif2_20km.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+# Change the mesh importer to the CUBIT importer.
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+# Set filename of mesh to import.
+filename = ../meshes/spbm_hex8_unif2_20km.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
+bc = [x_pos,x_neg,z_neg]
+
+# Set interfaces to an array of 1 fault: 'fault'.
+interfaces = [fault]
+
+[pylithapp.timedependent.implicit]
+# Set the output to an array of 2 output managers.
+# We will output the solution over the domain and the ground surface.
+output = [domain,subdomain]
+
+# Set subdomain component to OutputSolnSubset (subset of domain).
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+
+# We change spatial database for the rate of change from FixedDOFDB
+# (which has zero velocity) to UniformDB (which allows nonzero
+# velocities).
+
+# Set the parameters for the desired boundary conditions.
+#
+# The label corresponds to the nodeset ID in CUBIT.
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0, 1]
+label = 11
+db.label = Dirichlet BC on +x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on +x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of 1 cm/year.
+rate_db.data = [0.0,3.1688087814e-10]
+
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1]
+label = 12
+db.label = Dirichlet BC on -x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on -x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of -1 cm/year.
+rate_db.data = [0.0,-3.1688087814e-10]
+
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [2]
+label = 16
+db.label = Dirichlet BC on -z
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+
+# Set the parameters for the fault interface conditions.
+
+[pylithapp.timedependent.interfaces.fault]
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+mat_db.iohandler.filename = mat_elastic.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain.writer]
+filename = ../results/spbm_hex8_unif2_20km/spbm_hex8_unif2_20km.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK domain output of solution over ground surface.
+[pylithapp.problem.formulation.output.subdomain]
+label = 17 ; nodeset for subdomain
+writer.filename = ../results/spbm_hex8_unif2_20km/spbm_hex8_unif2_20km-groundsurf.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for vtk fault rupture output.
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = ../results/spbm_hex8_unif2_20km/spbm_hex8_unif2_20km-fault.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_unif2_20km/spbm_hex8_unif2_20km-statevars-elastic.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.maxwell.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_unif2_20km/spbm_hex8_unif2_20km-statevars-maxwell.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+

Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif2_20km.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif2_20km.sh	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif2_20km.sh	2009-03-03 00:40:52 UTC (rev 14192)
@@ -0,0 +1 @@
+pylith spbm_hex8_unif2_20km.cfg timedep.cfg


Property changes on: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif2_20km.sh
___________________________________________________________________
Name: svn:executable
   + *

Deleted: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.cfg	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.cfg	2009-03-03 00:40:52 UTC (rev 14192)
@@ -1,127 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies parameters specific to the boundary and
-# interface conditions. The general parameters are specificed in the
-# pylithapp.cfg file which PyLith reads by default. Information controlling
-# time stepping is contained in the file timedep.cfg.
-#
-# To run the simulation:
-# pylith spbm_hex8_unif_6750m.cfg timedep.cfg
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-# Change the mesh importer to the CUBIT importer.
-importer = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.importer]
-# Set filename of mesh to import.
-filename = ../meshes/spbm_hex8_unif_6750m.exo
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
-bc = [x_pos,x_neg,z_neg]
-
-# Set interfaces to an array of 1 fault: 'fault'.
-interfaces = [fault]
-
-[pylithapp.timedependent.implicit]
-# Set the output to an array of 2 output managers.
-# We will output the solution over the domain and the ground surface.
-output = [domain,subdomain]
-
-# Set subdomain component to OutputSolnSubset (subset of domain).
-output.subdomain = pylith.meshio.OutputSolnSubset
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-
-# We change spatial database for the rate of change from FixedDOFDB
-# (which has zero velocity) to UniformDB (which allows nonzero
-# velocities).
-
-# Set the parameters for the desired boundary conditions.
-#
-# The label corresponds to the nodeset ID in CUBIT.
-
-[pylithapp.timedependent.bc.x_pos]
-fixed_dof = [0, 1]
-label = 11
-db.label = Dirichlet BC on +x
-rate_db = spatialdata.spatialdb.UniformDB
-rate_db.label = Dirichlet rate BC on +x
-rate_db.values = [dof-0,dof-1]
-# This rate (in m/s) gives a value of -1 cm/year.
-rate_db.data = [0.0,-3.1688087814e-10]
-
-[pylithapp.timedependent.bc.x_neg]
-fixed_dof = [0, 1]
-label = 12
-db.label = Dirichlet BC on -x
-rate_db = spatialdata.spatialdb.UniformDB
-rate_db.label = Dirichlet rate BC on -x
-rate_db.values = [dof-0,dof-1]
-# This rate (in m/s) gives a value of 1 cm/year.
-rate_db.data = [0.0,3.1688087814e-10]
-
-[pylithapp.timedependent.bc.z_neg]
-fixed_dof = [2]
-label = 16
-db.label = Dirichlet BC on -z
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces]
-
-# Set the parameters for the fault interface conditions.
-
-[pylithapp.timedependent.interfaces.fault]
-label = 10
-quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-mat_db.iohandler.filename = mat_elastic.spatialdb
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.domain.writer]
-filename = ../results/spbm_hex8_unif_6750m.vtk
-time_format = %04.0f
-time_constant = 1.0*year
-
-# Give basename for VTK domain output of solution over ground surface.
-[pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
-writer.filename = ../results/spbm_hex8_unif_6750m-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Give basename for vtk fault rupture output.
-[pylithapp.timedependent.interfaces.fault.output.writer]
-filename = ../results/spbm_hex8_unif_6750m-fault.vtk
-time_format = %04.0f
-time_constant = 1.0*year
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = ../results/spbm_hex8_unif_6750m-statevars-elastic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-[pylithapp.timedependent.materials.maxwell.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = ../results/spbm_hex8_unif_6750m-statevars-maxwell.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-

Deleted: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.sh	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.sh	2009-03-03 00:40:52 UTC (rev 14192)
@@ -1 +0,0 @@
-pylith spbm_hex8_unif_6750m.cfg timedep.cfg

Modified: short/2.5D/benchmarks/savageprescott/parameters/timedep.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/timedep.cfg	2009-03-03 00:30:34 UTC (rev 14191)
+++ short/2.5D/benchmarks/savageprescott/parameters/timedep.cfg	2009-03-03 00:40:52 UTC (rev 14192)
@@ -2,12 +2,13 @@
 [pylithapp]
 
 # This is not a self-contained simulation configuration file. This
-# file only specifies parameters specific to the boundary and
+# 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.
+# 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_hex8_unif_8000m.cfg timedep.cfg
+# pylith spbm_xxxx.cfg timedep.cfg
 
 # ----------------------------------------------------------------------
 # faults
@@ -89,27 +90,27 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-# Specify that output for the entire domain occurs every 100 years.
+# Specify that output for the entire domain occurs every 50 years.
 [pylithapp.problem.formulation.output.domain]
 output_freq = time_step
-time_step = 100*year
+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 100 years.
+# Specify that output for the fault occurs every 50 years.
 [pylithapp.timedependent.interfaces.fault.output]
 output_freq = time_step
-time_step = 100*year
+time_step = 50*year
 
-# Specify that output for both materials occurs every 100 years.
+# Specify that output for both materials occurs every 50 years.
 [pylithapp.timedependent.materials.elastic.output]
 output_freq = time_step
-time_step = 100*year
+time_step = 50*year
 
 [pylithapp.timedependent.materials.maxwell.output]
 output_freq = time_step
-time_step = 100*year
+time_step = 50*year
 



More information about the CIG-COMMITS mailing list