[cig-commits] r19949 - in short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott: . parameters
willic3 at geodynamics.org
willic3 at geodynamics.org
Mon Apr 16 21:02:38 PDT 2012
Author: willic3
Date: 2012-04-16 21:02:38 -0700 (Mon, 16 Apr 2012)
New Revision: 19949
Added:
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/fault_lock20.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/fault_lock40.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/hex8.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/mat_elastic.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/mat_maxwell_tau0-10.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/sliprate_creep_lock20.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/sliptime.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/tet4.cfg
Removed:
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock20.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock40.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_elastic.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-2.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliprate_creep_lock20.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliptime.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/tet4.cfg
Log:
Partially finished cleaning up files for S&P benchmark.
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/fault_lock20.cfg (from rev 19947, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock20.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/fault_lock20.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/fault_lock20.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -0,0 +1,108 @@
+# -*- 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. Mesh information is
+# contained in hex8.cfg and tet4.cfg and solution parameters are contained
+# in asm.cfg and fieldsplit.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.label = Slip for rupture source 1
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 1
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# two
+[pylithapp.timedependent.interfaces.fault.eq_srcs.two.slip_function]
+slip.label = Slip for rupture source 2
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 2
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# three
+[pylithapp.timedependent.interfaces.fault.eq_srcs.three.slip_function]
+slip.label = Slip for rupture source 3
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 3
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# four
+[pylithapp.timedependent.interfaces.fault.eq_srcs.four.slip_function]
+slip.label = Slip for rupture source 4
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 4
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# five
+[pylithapp.timedependent.interfaces.fault.eq_srcs.five.slip_function]
+slip.label = Slip for rupture source 5
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 5
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# six
+[pylithapp.timedependent.interfaces.fault.eq_srcs.six.slip_function]
+slip.label = Slip for rupture source 6
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 6
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# seven
+[pylithapp.timedependent.interfaces.fault.eq_srcs.seven.slip_function]
+slip.label = Slip for rupture source 7
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 7
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# eight
+[pylithapp.timedependent.interfaces.fault.eq_srcs.eight.slip_function]
+slip.label = Slip for rupture source 8
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 8
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# nine
+[pylithapp.timedependent.interfaces.fault.eq_srcs.nine.slip_function]
+slip.label = Slip for rupture source 9
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 9
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# ten
+[pylithapp.timedependent.interfaces.fault.eq_srcs.ten.slip_function]
+slip.label = Slip for rupture source 10
+slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
+slip_time.label = Slip time for rupture source 10
+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.label = Slip rate for creep slip function
+slip_function.slip_rate.iohandler.filename = sliprate_creep_lock20.spatialdb
+slip_function.slip_time.label = Slip time for creep slip function
+slip_function.slip_time.iohandler.filename = sliptime.spatialdb
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/fault_lock40.cfg (from rev 19947, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock40.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/fault_lock40.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/fault_lock40.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -0,0 +1,108 @@
+# -*- 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. Mesh information is
+# contained in hex8.cfg and tet4.cfg and solution parameters are contained
+# in asm.cfg and fieldsplit.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.label = Slip for rupture source 1
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 1
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# two
+[pylithapp.timedependent.interfaces.fault.eq_srcs.two.slip_function]
+slip.label = Slip for rupture source 2
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 2
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# three
+[pylithapp.timedependent.interfaces.fault.eq_srcs.three.slip_function]
+slip.label = Slip for rupture source 3
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 3
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# four
+[pylithapp.timedependent.interfaces.fault.eq_srcs.four.slip_function]
+slip.label = Slip for rupture source 4
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 4
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# five
+[pylithapp.timedependent.interfaces.fault.eq_srcs.five.slip_function]
+slip.label = Slip for rupture source 5
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 5
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# six
+[pylithapp.timedependent.interfaces.fault.eq_srcs.six.slip_function]
+slip.label = Slip for rupture source 6
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 6
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# seven
+[pylithapp.timedependent.interfaces.fault.eq_srcs.seven.slip_function]
+slip.label = Slip for rupture source 7
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 7
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# eight
+[pylithapp.timedependent.interfaces.fault.eq_srcs.eight.slip_function]
+slip.label = Slip for rupture source 8
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 8
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# nine
+[pylithapp.timedependent.interfaces.fault.eq_srcs.nine.slip_function]
+slip.label = Slip for rupture source 9
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 9
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# ten
+[pylithapp.timedependent.interfaces.fault.eq_srcs.ten.slip_function]
+slip.label = Slip for rupture source 10
+slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
+slip_time.label = Slip time for rupture source 10
+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.label = Slip rate for creep slip function
+slip_function.slip_rate.iohandler.filename = sliprate_creep_lock40.spatialdb
+slip_function.slip_time.label = Slip time for creep slip function
+slip_function.slip_time.iohandler.filename = sliptime.spatialdb
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/hex8.cfg (from rev 19947, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/hex8.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/hex8.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -0,0 +1,31 @@
+# -*- 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 files fault_lock20.cfg and fault_lock40.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/mat_elastic.spatialdb (from rev 19947, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_elastic.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/mat_elastic.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/mat_elastic.spatialdb 2012-04-17 04:02:38 UTC (rev 19949)
@@ -0,0 +1,26 @@
+// -*- 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.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3 // number of material property values
+ value-names = density vs vp // names of the material property values
+ value-units = kg/m**3 m/s m/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)
+0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/mat_maxwell_tau0-10.spatialdb (from rev 19947, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-2.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/mat_maxwell_tau0-10.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/mat_maxwell_tau0-10.spatialdb 2012-04-17 04:02:38 UTC (rev 19949)
@@ -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 100 years. This
+// corresponds to a value of tau0 (recurrence_time/relaxation_time) of 2.
+//
+#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 4.73364e19
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock20.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock20.cfg 2012-04-16 21:29:35 UTC (rev 19948)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock20.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -1,108 +0,0 @@
-# -*- 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. Mesh information is
-# contained in hex8.cfg and tet4.cfg and solution parameters are contained
-# in asm.cfg and fieldsplit.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.label = Slip for rupture source 1
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 1
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# two
-[pylithapp.timedependent.interfaces.fault.eq_srcs.two.slip_function]
-slip.label = Slip for rupture source 2
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 2
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# three
-[pylithapp.timedependent.interfaces.fault.eq_srcs.three.slip_function]
-slip.label = Slip for rupture source 3
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 3
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# four
-[pylithapp.timedependent.interfaces.fault.eq_srcs.four.slip_function]
-slip.label = Slip for rupture source 4
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 4
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# five
-[pylithapp.timedependent.interfaces.fault.eq_srcs.five.slip_function]
-slip.label = Slip for rupture source 5
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 5
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# six
-[pylithapp.timedependent.interfaces.fault.eq_srcs.six.slip_function]
-slip.label = Slip for rupture source 6
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 6
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# seven
-[pylithapp.timedependent.interfaces.fault.eq_srcs.seven.slip_function]
-slip.label = Slip for rupture source 7
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 7
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# eight
-[pylithapp.timedependent.interfaces.fault.eq_srcs.eight.slip_function]
-slip.label = Slip for rupture source 8
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 8
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# nine
-[pylithapp.timedependent.interfaces.fault.eq_srcs.nine.slip_function]
-slip.label = Slip for rupture source 9
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 9
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# ten
-[pylithapp.timedependent.interfaces.fault.eq_srcs.ten.slip_function]
-slip.label = Slip for rupture source 10
-slip.iohandler.filename = finalslip_rupture_lock20.spatialdb
-slip_time.label = Slip time for rupture source 10
-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.label = Slip rate for creep slip function
-slip_function.slip_rate.iohandler.filename = sliprate_creep_lock20.spatialdb
-slip_function.slip_time.label = Slip time for creep slip function
-slip_function.slip_time.iohandler.filename = sliptime.spatialdb
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock40.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock40.cfg 2012-04-16 21:29:35 UTC (rev 19948)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fault_lock40.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -1,108 +0,0 @@
-# -*- 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. Mesh information is
-# contained in hex8.cfg and tet4.cfg and solution parameters are contained
-# in asm.cfg and fieldsplit.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.label = Slip for rupture source 1
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 1
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# two
-[pylithapp.timedependent.interfaces.fault.eq_srcs.two.slip_function]
-slip.label = Slip for rupture source 2
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 2
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# three
-[pylithapp.timedependent.interfaces.fault.eq_srcs.three.slip_function]
-slip.label = Slip for rupture source 3
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 3
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# four
-[pylithapp.timedependent.interfaces.fault.eq_srcs.four.slip_function]
-slip.label = Slip for rupture source 4
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 4
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# five
-[pylithapp.timedependent.interfaces.fault.eq_srcs.five.slip_function]
-slip.label = Slip for rupture source 5
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 5
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# six
-[pylithapp.timedependent.interfaces.fault.eq_srcs.six.slip_function]
-slip.label = Slip for rupture source 6
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 6
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# seven
-[pylithapp.timedependent.interfaces.fault.eq_srcs.seven.slip_function]
-slip.label = Slip for rupture source 7
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 7
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# eight
-[pylithapp.timedependent.interfaces.fault.eq_srcs.eight.slip_function]
-slip.label = Slip for rupture source 8
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 8
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# nine
-[pylithapp.timedependent.interfaces.fault.eq_srcs.nine.slip_function]
-slip.label = Slip for rupture source 9
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 9
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# ten
-[pylithapp.timedependent.interfaces.fault.eq_srcs.ten.slip_function]
-slip.label = Slip for rupture source 10
-slip.iohandler.filename = finalslip_rupture_lock40.spatialdb
-slip_time.label = Slip time for rupture source 10
-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.label = Slip rate for creep slip function
-slip_function.slip_rate.iohandler.filename = sliprate_creep_lock40.spatialdb
-slip_function.slip_time.label = Slip time for creep slip function
-slip_function.slip_time.iohandler.filename = sliptime.spatialdb
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8.cfg 2012-04-16 21:29:35 UTC (rev 19948)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -1,31 +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 files fault_lock20.cfg and fault_lock40.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
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_elastic.spatialdb 2012-04-16 21:29:35 UTC (rev 19948)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_elastic.spatialdb 2012-04-17 04:02:38 UTC (rev 19949)
@@ -1,26 +0,0 @@
-// -*- 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.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 3 // number of material property values
- value-names = density vs vp // names of the material property values
- value-units = kg/m**3 m/s m/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)
-0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258
Deleted: 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 2012-04-16 21:29:35 UTC (rev 19948)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/mat_maxwell_tau0-2.spatialdb 2012-04-17 04:02:38 UTC (rev 19949)
@@ -1,30 +0,0 @@
-// -*- 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 100 years. This
-// corresponds to a value of tau0 (recurrence_time/relaxation_time) of 2.
-//
-#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 4.73364e19
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg 2012-04-16 21:29:35 UTC (rev 19948)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -1,200 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This file contains the problem parameters common to all the examples,
-# and is read automatically. To run a specific example, it is also
-# necessary to specify the file describing the mesh, the file describing
-# the time-dependent fault behavior, the file describing the solution
-# parameters, and the problem-specific file.
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-# Turn on some journals to show progress.
-[pylithapp.journal.info]
-timedependent = 1
-implicit = 1
-petsc = 1
-solverlinear = 1
-meshimporter = 1
-meshiocubit = 1
-meshiolagrit = 1
-gmvfile = 1
-psetfile = 1
-implicitelasticity = 1
-faultcohesivekin = 1
-#quadrature3d = 1
-#fiatlagrange = 1
-
-[pylithapp.journal.debug]
-# pylithapp = 1
-# problem = 1
-# implicit = 1
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-# Change the mesh reader to the CUBIT reader.
-reader = pylith.meshio.MeshIOCubit
-# reorder_mesh = true
-
-[pylithapp.mesh_generator.reader]
-coordsys.space_dim = 3
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-bc = [x_pos,x_neg,z_neg,y_pos,y_neg]
-
-interfaces = [fault]
-
-[pylithapp.timedependent.implicit]
-# 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
-
-# ----------------------------------------------------------------------
-# time stepping
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.formulation.time_step]
-# Define the total time for the simulation and the default time step size.
-total_time = 2000.0*year ; total time of simulation
-dt = 5.0*year ; time step
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-# Set materials to an array of 2 materials:
-# 'elastic' and 'maxwell'
-materials = [elastic,maxwell]
-materials.maxwell = pylith.materials.MaxwellIsotropic3D
-
-[pylithapp.timedependent.materials.elastic]
-label = Upper elastic layer
-id = 1
-db_properties.iohandler.filename = mat_elastic.spatialdb
-db_properties.label = Elastic material properties
-
-[pylithapp.timedependent.materials.maxwell]
-label = Lower Maxwell material
-id = 2
-db_properties.iohandler.filename = mat_maxwell_tau0-2.spatialdb
-db_properties.label = Viscoelastic material properties
-output.cell_data_fields = [total_strain,viscous_strain,stress]
-output.cell_info_fields = [mu,lambda,density,maxwell_time]
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-bc_dof = [0, 1]
-label = face_xpos
-db_initial.label = Dirichlet BC on +x
-# We change spatial database for the rate of change from FixedDOFDB
-# (which has zero velocity) to UniformDB (which allows nonzero
-# velocities).
-db_rate = spatialdata.spatialdb.UniformDB
-db_rate.label = Dirichlet rate BC on +x
-db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
-# Apply a rate of +1 cm/year.
-db_rate.data = [0.0*cm/year,1.0*cm/year,0.0*year]
-
-[pylithapp.timedependent.bc.x_neg]
-bc_dof = [0, 1]
-label = face_xneg
-db_initial.label = Dirichlet BC on -x
-# We change spatial database for the rate of change from FixedDOFDB
-# (which has zero velocity) to UniformDB (which allows nonzero
-# velocities).
-db_rate = spatialdata.spatialdb.UniformDB
-db_rate.label = Dirichlet rate BC on -x
-db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
-# Apply a rate of -1 cm/year.
-db_rate.data = [0.0*cm/year,-1.0*cm/year,0.0*year]
-
-# Fix z-displacements on bottom boundary.
-[pylithapp.timedependent.bc.z_neg]
-bc_dof = [2]
-label = face_zneg
-db_initial.label = Dirichlet BC on -z
-
-# Fix x-displacements on y-boundaries to prevent movement perpendicular to
-# the fault.
-[pylithapp.timedependent.bc.y_pos]
-bc_dof = [0]
-label = face_ypos_nofault
-db_initial.label = Dirichlet BC on +y
-
-[pylithapp.timedependent.bc.y_neg]
-bc_dof = [0]
-label = face_yneg_nofault
-db_initial.label = Dirichlet BC on -y
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-# Set the parameters for the fault interface conditions.
-[pylithapp.timedependent.interfaces.fault]
-label = fault
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.domain]
-output_freq = time_step
-time_step = 50*year
-writer = pylith.meshio.DataWriterHDF5Mesh
-
-# Give basename for VTK domain output of solution over ground surface.
-[pylithapp.problem.formulation.output.subdomain]
-label = face_zpos ; nodeset for subdomain
-vertex_data_fields = [displacement,velocity]
-output_freq = time_step
-time_step = 10*year
-writer = pylith.meshio.DataWriterHDF5SubMesh
-
-# Give basename for vtk fault rupture output.
-[pylithapp.timedependent.interfaces.fault.output]
-output_freq = time_step
-time_step = 10*year
-writer = pylith.meshio.DataWriterHDF5SubSubMesh
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-output_freq = time_step
-time_step = 50*year
-cell_filter = pylith.meshio.CellFilterAvgMesh
-writer = pylith.meshio.DataWriterHDF5Mesh
-
-[pylithapp.timedependent.materials.maxwell.output]
-output_freq = time_step
-time_step = 50*year
-cell_filter = pylith.meshio.CellFilterAvgMesh
-writer = pylith.meshio.DataWriterHDF5Mesh
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-# Set the solver options.
-
-[pylithapp.petsc]
-
-ksp_rtol = 1.0e-8
-ksp_max_it = 1000
-ksp_gmres_restart = 100
-ksp_monitor_singular_value = true
-ksp_converged_reason = true
-
-ksp_monitor = true
-ksp_view = true
-
-log_summary = true
-# start_in_debugger = true
-
-# End of file
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliprate_creep_lock20.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliprate_creep_lock20.spatialdb 2012-04-16 21:29:35 UTC (rev 19948)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliprate_creep_lock20.spatialdb 2012-04-17 04:02:38 UTC (rev 19949)
@@ -1,36 +0,0 @@
-// -*- 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 < -16.0 km and
-// zero slip for z >= -16.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 = 8
- 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 -19.00 0.00 0.0 0.0
-0.0 0.0 -20.00 1.00 0.0 0.0
-0.0 0.0 -21.00 2.00 0.0 0.0
-0.0 0.0 -39.00 2.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/sliptime.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliptime.spatialdb 2012-04-16 21:29:35 UTC (rev 19948)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliptime.spatialdb 2012-04-17 04:02:38 UTC (rev 19949)
@@ -1,26 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of the slip
-// initiation time over the fault surface for the kinematic fault
-// rupture. We specify a slip time of 0.0 s so that the slip occurs at
-// the origin time of the earthquake rupture.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 1
- value-names = slip-time
- value-units = s
- num-locs = 1
- data-dim = 0 // data is uniform (1 location)
- 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) Slip initiation time (s)
-0.0 0.0 0.0 0.0
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/tet4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/tet4.cfg 2012-04-16 21:29:35 UTC (rev 19948)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/tet4.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -1,31 +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 files fault_lock20.cfg and fault_lock40.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
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/pylithapp.cfg (from rev 19947, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/pylithapp.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -0,0 +1,201 @@
+# -*- Python -*-
+[pylithapp]
+
+# This file contains the problem parameters common to all the examples,
+# and is read automatically. To run a specific example, it is also
+# necessary to specify the file describing the mesh, the file describing
+# the time-dependent fault behavior, the file describing the solution
+# parameters, and the problem-specific file.
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on some journals to show progress.
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshimporter = 1
+meshiocubit = 1
+meshiolagrit = 1
+gmvfile = 1
+psetfile = 1
+implicitelasticity = 1
+faultcohesivekin = 1
+#quadrature3d = 1
+#fiatlagrange = 1
+
+[pylithapp.journal.debug]
+# pylithapp = 1
+# problem = 1
+# implicit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+# Change the mesh reader to the CUBIT reader.
+reader = pylith.meshio.MeshIOCubit
+reorder_mesh = true
+
+[pylithapp.mesh_generator.reader]
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+bc = [x_pos,x_neg,z_neg,y_pos,y_neg]
+
+interfaces = [fault]
+
+[pylithapp.timedependent.implicit]
+# 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
+
+# ----------------------------------------------------------------------
+# time stepping
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+# Define the total time for the simulation and the default time step size.
+total_time = 2000.0*year ; total time of simulation
+dt = 5.0*year ; time step
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set materials to an array of 2 materials:
+# 'elastic' and 'maxwell'
+materials = [elastic,maxwell]
+materials.maxwell = pylith.materials.MaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.elastic]
+label = Upper elastic layer
+id = 1
+db_properties.iohandler.filename = mat_elastic.spatialdb
+db_properties.label = Elastic material properties
+
+[pylithapp.timedependent.materials.maxwell]
+label = Lower Maxwell material
+id = 2
+db_properties.iohandler.filename = mat_maxwell_tau0-10.spatialdb
+db_properties.label = Viscoelastic material properties
+output.cell_data_fields = [total_strain,viscous_strain,stress]
+output.cell_info_fields = [mu,lambda,density,maxwell_time]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0, 1]
+label = face_xpos
+db_initial.label = Dirichlet BC on +x
+# We change spatial database for the rate of change from FixedDOFDB
+# (which has zero velocity) to UniformDB (which allows nonzero
+# velocities).
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
+# Apply a rate of +1 cm/year.
+db_rate.data = [0.0*cm/year,1.0*cm/year,0.0*year]
+
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0, 1]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+# We change spatial database for the rate of change from FixedDOFDB
+# (which has zero velocity) to UniformDB (which allows nonzero
+# velocities).
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on -x
+db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
+# Apply a rate of -1 cm/year.
+db_rate.data = [0.0*cm/year,-1.0*cm/year,0.0*year]
+
+# Fix z-displacements on bottom boundary.
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# Fix x-displacements on y-boundaries to prevent movement perpendicular to
+# the fault.
+[pylithapp.timedependent.bc.y_pos]
+bc_dof = [0]
+label = face_ypos_nofault
+db_initial.label = Dirichlet BC on +y
+
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0]
+label = face_yneg_nofault
+db_initial.label = Dirichlet BC on -y
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+# Set the parameters for the fault interface conditions.
+[pylithapp.timedependent.interfaces.fault]
+label = fault
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+output_freq = time_step
+time_step = 50*year
+writer = pylith.meshio.DataWriterHDF5Mesh
+
+# Give basename for VTK domain output of solution over ground surface.
+[pylithapp.problem.formulation.output.subdomain]
+label = face_zpos ; nodeset for subdomain
+vertex_data_fields = [displacement,velocity]
+output_freq = time_step
+time_step = 10*year
+writer = pylith.meshio.DataWriterHDF5SubMesh
+
+# Give basename for vtk fault rupture output.
+[pylithapp.timedependent.interfaces.fault.output]
+output_freq = time_step
+time_step = 10*year
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+output_freq = time_step
+time_step = 50*year
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer = pylith.meshio.DataWriterHDF5Mesh
+
+[pylithapp.timedependent.materials.maxwell.output]
+output_freq = time_step
+time_step = 50*year
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer = pylith.meshio.DataWriterHDF5Mesh
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+# Set the solver options.
+
+[pylithapp.petsc]
+
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-20
+ksp_max_it = 500
+
+ksp_monitor = true
+ksp_view = true
+ksp_monitor_singular_value = true
+ksp_converged_reason = true
+ksp_monitor_true_residual = true
+
+log_summary = true
+# start_in_debugger = true
+
+# End of file
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/sliprate_creep_lock20.spatialdb (from rev 19947, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliprate_creep_lock20.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/sliprate_creep_lock20.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/sliprate_creep_lock20.spatialdb 2012-04-17 04:02:38 UTC (rev 19949)
@@ -0,0 +1,36 @@
+// -*- 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 < -16.0 km and
+// zero slip for z >= -16.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 = 8
+ 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 -19.00 0.00 0.0 0.0
+0.0 0.0 -20.00 1.00 0.0 0.0
+0.0 0.0 -21.00 2.00 0.0 0.0
+0.0 0.0 -39.00 2.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
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/sliptime.spatialdb (from rev 19947, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/sliptime.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/sliptime.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/sliptime.spatialdb 2012-04-17 04:02:38 UTC (rev 19949)
@@ -0,0 +1,26 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the slip
+// initiation time over the fault surface for the kinematic fault
+// rupture. We specify a slip time of 0.0 s so that the slip occurs at
+// the origin time of the earthquake rupture.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 1
+ value-names = slip-time
+ value-units = s
+ num-locs = 1
+ data-dim = 0 // data is uniform (1 location)
+ 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) Slip initiation time (s)
+0.0 0.0 0.0 0.0
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/tet4.cfg (from rev 19947, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/tet4.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/tet4.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/tet4.cfg 2012-04-17 04:02:38 UTC (rev 19949)
@@ -0,0 +1,31 @@
+# -*- 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 files fault_lock20.cfg and fault_lock40.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