[cig-commits] r22291 - short/3D/PyLith/trunk/examples/bar_shearwave/quad4

[brad at geodynamics.org]

Author: brad
Date: 2013-06-14 15:56:53 -0700 (Fri, 14 Jun 2013)
New Revision: 22291

Modified:
   short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
   short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup.cfg
   short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup_ratestateageing.cfg
   short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup_slipweakening.cfg
Log:
Manual patch from v1.7-trunk.

Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg	2013-06-14 22:53:42 UTC (rev 22290)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg	2013-06-14 22:56:53 UTC (rev 22291)
@@ -43,11 +43,10 @@
 [pylithapp.timedependent]
 # Set the spatial dimension of the problem.
 dimension = 2
-elastic_prestep = False
 
-# Change to an explicit time stepping formulation.  We use explicit
-# time-stepping with a consistent (full) Jacobian matrix.
+# Change to an explicit time stepping formulation with no elastic prestep.
 formulation = pylith.problems.Explicit
+elastic_prestep = False
 
 # Nondimensionalize problem using wave propagation parameters.
 normalizer = spatialdata.units.NondimElasticDynamic

Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup.cfg	2013-06-14 22:53:42 UTC (rev 22290)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup.cfg	2013-06-14 22:56:53 UTC (rev 22291)
@@ -19,17 +19,31 @@
 fault = pylith.faults.FaultCohesiveDyn
 
 [pylithapp.timedependent.interfaces.fault]
-# Specify the initial tractions on the fault using a uniform DB.
-#  shear: 6.1 MPa (right-lateral)
-#  normal 10 MPa (compressive)
+# Specify the initial tractions and a temporal perturbation on the
+# fault using a uniform DB.
+#
+# Initial
+#   shear: 75.0 MPa (right-lateral)
+#   normal 120 MPa (compressive)
+# Perturbation
+#   shear: 25.0 MPa (right-lateral)
 traction_perturbation = pylith.faults.TractPerturbation
 
 [pylithapp.timedependent.interfaces.fault.traction_perturbation]
 db_initial = spatialdata.spatialdb.UniformDB
 db_initial.label = Initial fault tractions
 db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [-6.1*MPa, -10.0*MPa]
+db_initial.data = [-75.0*MPa, -120.0*MPa]
 
+db_change = spatialdata.spatialdb.UniformDB
+db_change.label = Traction perturbation
+db_change.values = [traction-shear,traction-normal,change-start-time]
+db_change.data = [-25.0*MPa, 0.0*MPa, 0.0*s]
+
+th_change = spatialdata.spatialdb.TimeHistory
+th_change.label = Traction perturbation time history
+th_change.filename = traction_change.timedb
+
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------

Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup_ratestateageing.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup_ratestateageing.cfg	2013-06-14 22:53:42 UTC (rev 22290)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup_ratestateageing.cfg	2013-06-14 22:56:53 UTC (rev 22291)
@@ -44,13 +44,13 @@
 friction.db_properties = spatialdata.spatialdb.UniformDB
 friction.db_properties.label = Rate State Ageing
 friction.db_properties.values = [reference-friction-coefficient,reference-slip-rate,characteristic-slip-distance,constitutive-parameter-a,constitutive-parameter-b,cohesion]
-friction.db_properties.data = [0.6,1.0e-6*m/s,0.0370*m,0.0125,0.0172,0.0*Pa]
+friction.db_properties.data = [0.6,1.0e-6*m/s,0.02*m,0.008,0.012,0.0*Pa]
 
 # Set spatial database for the initial value of the state variable.
 friction.db_initial_state = spatialdata.spatialdb.UniformDB
 friction.db_initial_state.label = Rate State Ageing State
 friction.db_initial_state.values = [state-variable]
-friction.db_initial_state.data = [92.7*s]
+friction.db_initial_state.data = [1.60623899921e+9*s]
 
 # ----------------------------------------------------------------------
 # output
@@ -60,7 +60,8 @@
 writer.filename = output/ratestateageing.vtk
 
 [pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/ratestateageing-fault.vtk
+vertex_data_fields = [slip, slip_rate, traction, state_variable]
+writer.filename = output/ratestateageing-faut.vtk
 
 [pylithapp.timedependent.materials.elastic.output]
 writer.filename = output/ratestateageing-statevars.vtk

Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup_slipweakening.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup_slipweakening.cfg	2013-06-14 22:53:42 UTC (rev 22290)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/spontaneousrup_slipweakening.cfg	2013-06-14 22:56:53 UTC (rev 22291)
@@ -56,9 +56,7 @@
 # default values.
 vertex_info_fields = [strike_dir,normal_dir,traction_initial_value,static_coefficient,dynamic_coefficient,slip_weakening_parameter,cohesion]
 
-# Add output of state variables cumulative slip and previous slip to default
-# values.
-vertex_data_fields = [slip,traction,cumulative_slip,previous_slip]
+vertex_data_fields = [slip,slip_rate,traction]
 
 
 [pylithapp.timedependent.materials.elastic.output]