[cig-commits] r18510 - short/3D/PyLith/trunk/tests/2d/frictionslide

[brad at geodynamics.org]

Author: brad
Date: 2011-05-31 13:55:24 -0700 (Tue, 31 May 2011)
New Revision: 18510

Modified:
   short/3D/PyLith/trunk/tests/2d/frictionslide/plot_friction.py
   short/3D/PyLith/trunk/tests/2d/frictionslide/ratestate_weak.cfg
   short/3D/PyLith/trunk/tests/2d/frictionslide/velocitysteps.timedb
Log:
Updated weakening test case.

Modified: short/3D/PyLith/trunk/tests/2d/frictionslide/plot_friction.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/frictionslide/plot_friction.py	2011-05-31 18:22:47 UTC (rev 18509)
+++ short/3D/PyLith/trunk/tests/2d/frictionslide/plot_friction.py	2011-05-31 20:55:24 UTC (rev 18510)
@@ -1,6 +1,6 @@
 #!/usr/bin/env python
 
-sim = "ratestate_stable"
+sim = "ratestate_weak"
 
 # ======================================================================
 import tables
@@ -12,8 +12,14 @@
 dt = 0.01
 t = numpy.arange(0.0, 12.001, dt)
 mu0 = 0.6
-a = 0.016
-b = 0.012
+if sim == "ratestate_stable":
+    a = 0.016
+    b = 0.012
+elif sim == "ratestate_weak":
+    a = 0.008
+    b = 0.012
+else:
+    raise ValueError("Unknown sim '%s'." % sim)
 L = 2.0e-6
 V0 = 1.0e-6
 
@@ -40,7 +46,7 @@
 integrateStateVar(theta, V3, 4.0)
 
 mask4 = numpy.bitwise_and(t >= 8.0, t < 12.0)
-V4 = 2.0e-7
+V4 = 1.0e-7
 integrateStateVar(theta, V4, 8.0)
 
 slipRateE = mask1*V1 + mask2*V2 + mask3*V3 + mask4*V4

Modified: short/3D/PyLith/trunk/tests/2d/frictionslide/ratestate_weak.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/frictionslide/ratestate_weak.cfg	2011-05-31 18:22:47 UTC (rev 18509)
+++ short/3D/PyLith/trunk/tests/2d/frictionslide/ratestate_weak.cfg	2011-05-31 20:55:24 UTC (rev 18510)
@@ -23,23 +23,23 @@
 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.01*m,0.008,0.012,0.0*Pa]
+friction.db_properties.data = [0.6,1.0e-6*m/s,2.0e-6*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 = [1.0e+4*s]
+friction.db_initial_state.data = [2.0e+0*s]
 
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
 # Set filenames for output.
 [pylithapp.problem.formulation.output.output]
-writer.filename = output/ratestate_weak.vtk
+writer.filename = output/ratestate_weak.h5
 
 [pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/ratestate_weak-fault.vtk
+writer.filename = output/ratestate_weak-fault.h5
 
 [pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/ratestate_weak-statevars.vtk
+writer.filename = output/ratestate_weak-statevars.h5

Modified: short/3D/PyLith/trunk/tests/2d/frictionslide/velocitysteps.timedb
===================================================================
--- short/3D/PyLith/trunk/tests/2d/frictionslide/velocitysteps.timedb	2011-05-31 18:22:47 UTC (rev 18509)
+++ short/3D/PyLith/trunk/tests/2d/frictionslide/velocitysteps.timedb	2011-05-31 20:55:24 UTC (rev 18510)
@@ -1,11 +1,12 @@
 #TIME HISTORY ascii
 TimeHistory {
-  num-points = 6
+  num-points = 7
   time-units = second
 }
  0.0  0.0
- 2.0  2.0e-6
- 4.0  4.2e-5
- 8.0  6.2e-05
-12.0  6.28e-05
-99.9  6.28e-05
+ 2.0  2.0e-6 // slip rate = 1.0e-6
+ 4.0  4.2e-5 // slip rate = 2.0e-5
+ 8.0  6.2e-05 // slip rate = 5.0e-6
+12.0  6.24e-05 // slip rate = 1.0e-7
+16.0  6.244e-05 // slip rate = 1.0e-8
+99.9  6.244e-05 // slip rate = 0