[cig-commits] r19083 - short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide

[brad at geodynamics.org]

Author: brad
Date: 2011-10-17 10:54:06 -0700 (Mon, 17 Oct 2011)
New Revision: 19083

Modified:
   short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/plot_friction.py
   short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/pylithapp.cfg
   short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/velocitysteps.timedb
Log:
Small tweaks to friction slide rate and state test.

Modified: short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/plot_friction.py
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/plot_friction.py	2011-10-15 22:56:47 UTC (rev 19082)
+++ short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/plot_friction.py	2011-10-17 17:54:06 UTC (rev 19083)
@@ -1,6 +1,6 @@
 #!/usr/bin/env python
 
-sim = "ratestate_stable"
+sim = "ratestate_weak"
 
 # ======================================================================
 import tables
@@ -10,7 +10,7 @@
 
 # ----------------------------------------------------------------------
 dt = 0.01
-t = numpy.arange(0.0, 12.001, dt)
+t = numpy.arange(0.0, 14.001, dt)
 mu0 = 0.6
 if sim == "ratestate_stable":
     a = 0.016
@@ -31,25 +31,33 @@
     return
 
 
-theta = L/V0*numpy.ones(t.shape)
 
 mask1 = t < 2.0
 V1 = 1.0e-6
+theta = L/V1*numpy.ones(t.shape)
 integrateStateVar(theta, V1, 0.0)
 
 mask2 = numpy.bitwise_and(t >= 2.0, t < 4.0)
-V2 = 2.0e-5
+V2 = 1.0e-5
 integrateStateVar(theta, V2, 2.0)
 
-mask3 = numpy.bitwise_and(t >= 4.0, t < 8.0)
-V3 = 5.0e-6
+mask3 = numpy.bitwise_and(t >= 4.0, t < 6.0)
+V3 = 4.0e-6
 integrateStateVar(theta, V3, 4.0)
 
-mask4 = numpy.bitwise_and(t >= 8.0, t < 12.0)
-V4 = 2.0e-7
-integrateStateVar(theta, V4, 8.0)
+mask4 = numpy.bitwise_and(t >= 6.0, t < 8.0)
+V4 = 2.0e-5
+integrateStateVar(theta, V4, 6.0)
 
-slipRateE = mask1*V1 + mask2*V2 + mask3*V3 + mask4*V4
+mask5 = numpy.bitwise_and(t >= 8.0, t < 12.0)
+V5 = 5.0e-6
+integrateStateVar(theta, V5, 8.0)
+
+mask6 = t >= 12.0
+V6 = 1.0e-6
+integrateStateVar(theta, V6, 12.0)
+
+slipRateE = mask1*V1 + mask2*V2 + mask3*V3 + mask4*V4 + mask5*V5 + mask6*V6
 stateVarE = theta
 
 muE = mu0 + a*numpy.log(slipRateE/V0) + b*numpy.log(V0*stateVarE/L)
@@ -69,7 +77,7 @@
 p = 2
 
 ax = pylab.subplot(1, 4, 1)
-ax.plot(time, slip[:,p,0])
+ax.plot(time, slip[:,p,0], 'r--')
 
 ax = pylab.subplot(1, 4, 2)
 ax.plot(t, numpy.log10(numpy.abs(slipRateE)), 'b-',

Modified: short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/pylithapp.cfg	2011-10-15 22:56:47 UTC (rev 19082)
+++ short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/pylithapp.cfg	2011-10-17 17:54:06 UTC (rev 19083)
@@ -46,8 +46,8 @@
 solver = pylith.problems.SolverNonlinear
 
 [pylithapp.timedependent.implicit.time_step]
-total_time = 12.0*s
-dt = 0.10*s
+total_time = 14.0*s
+dt = 0.05*s
 
 # ----------------------------------------------------------------------
 # materials

Modified: short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/velocitysteps.timedb
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/velocitysteps.timedb	2011-10-15 22:56:47 UTC (rev 19082)
+++ short/3D/PyLith/branches/v1.6-revisedfault/tests/2d/frictionslide/velocitysteps.timedb	2011-10-17 17:54:06 UTC (rev 19083)
@@ -5,8 +5,8 @@
 }
  0.0  0.0
  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.28e-05 // slip rate = 2.0e-7
-16.0  6.244e-05 // slip rate = 1.0e-8
-99.9  6.244e-05 // slip rate = 0
+ 4.0  2.2e-5 // slip rate = 1.0e-5
+ 6.0  3.0e-05 // slip rate = 4.0e-6
+ 8.0  7.0e-05 // slip rate = 2.0e-5
+12.0  9.0e-05 // slip rate = 5.0e-6
+22.0  1.0e-04 // slip rate = 1.0e-6