[cig-commits] r21394 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22

[brad at geodynamics.org]

Author: brad
Date: 2013-02-25 16:21:10 -0800 (Mon, 25 Feb 2013)
New Revision: 21394

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/create_friction.py
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/pylithapp.cfg
Log:
Added missing t0 parameter for new slip-time weakening friction.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/create_friction.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/create_friction.py	2013-02-26 00:20:32 UTC (rev 21393)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/create_friction.py	2013-02-26 00:21:10 UTC (rev 21394)
@@ -11,7 +11,7 @@
 # Create spatial variation of friction parameters for TPV22-23 on the
 # faults.
 
-hypoy = -5.0e+3
+hypoy = -10.0e+3
 hypoz = -10.0e+3
 vs = 3464.0
 
@@ -44,6 +44,7 @@
 coefStatic = 0.548*numpy.ones( (nptsy*nptsz), dtype=numpy.float64)
 coefDyn = 0.373*numpy.ones( (nptsy*nptsz), dtype=numpy.float64)
 d0 = 0.30*numpy.ones( (nptsy*nptsz), dtype=numpy.float64)
+t0 = 0.50*numpy.ones( (nptsy*nptsz), dtype=numpy.float64)
 
 mask = points[:,2] > -5000.0
 cohesion = mask*0.0014*(+5.0e+3+points[:,2]) + ~mask*0.0
@@ -76,9 +77,12 @@
                       {'name': 'cohesion', 
                        'units': 'Pa',
                        'data': cohesion},
-                      {'name': 'weakening-time', 
+                      {'name': 'time-weakening-time', 
                        'units': 's',
                        'data': weakTime},
+                      {'name': 'time-weakening-parameter', 
+                       'units': 's',
+                       'data': t0},
                       ],
            }
 writer.write(dataOut)
@@ -109,6 +113,7 @@
 coefDyn = 0.373*numpy.ones( (nptsy*nptsz), dtype=numpy.float64)
 d0 = 0.30*numpy.ones( (nptsy*nptsz), dtype=numpy.float64)
 weakTime = 1.0e+9*numpy.ones( (nptsy*nptsz), dtype=numpy.float64)
+t0 = 0.50*numpy.ones( (nptsy*nptsz), dtype=numpy.float64)
 
 mask = points[:,2] > -5000.0
 cohesion = mask*0.0014*(+5.0e+3+points[:,2]) + ~mask*0.0
@@ -133,9 +138,12 @@
                       {'name': 'cohesion', 
                        'units': 'Pa',
                        'data': cohesion},
-                      {'name': 'weakening-time', 
+                      {'name': 'time-weakening-time', 
                        'units': 's',
                        'data': weakTime},
+                      {'name': 'time-weakening-parameter', 
+                       'units': 's',
+                       'data': t0},
                       ],
            }
 writer.write(dataOut)

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/pylithapp.cfg	2013-02-26 00:20:32 UTC (rev 21393)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/pylithapp.cfg	2013-02-26 00:21:10 UTC (rev 21394)
@@ -115,7 +115,7 @@
 id = 100
 label = fault_main
 
-friction = pylith.friction.SlipWeakeningTime
+friction = pylith.friction.SlipWeakeningTimeStable
 friction.label = Slip Weakening on main fault
 friction.db_properties = spatialdata.spatialdb.SimpleDB
 friction.db_properties.iohandler.filename = friction_main.spatialdb
@@ -134,7 +134,7 @@
 id = 101
 label = fault_ext
 
-friction = pylith.friction.SlipWeakeningTime
+friction = pylith.friction.SlipWeakeningTimeStable
 friction.label = Slip Weakening on main fault
 friction.db_properties = spatialdata.spatialdb.SimpleDB
 friction.db_properties.iohandler.filename = friction_stepover.spatialdb
@@ -164,14 +164,14 @@
 writer = pylith.meshio.DataWriterHDF5ExtSubMesh
 
 [pylithapp.timedependent.interfaces.fault_main.output]
-vertex_info_fields = [traction_initial_value,static_coefficient,dynamic_coefficient,slip_weakening_parameter,weakening_time]
+vertex_info_fields = [traction_initial_value,static_coefficient,dynamic_coefficient,slip_weakening_parameter,time_weakening_time,time_weakening_parameter]
 vertex_data_fields = [slip,slip_rate,traction]
 output_freq = time_step
 time_step = 0.04999*s
 writer = pylith.meshio.DataWriterHDF5ExtSubSubMesh
 
 [pylithapp.timedependent.interfaces.fault_stepover.output]
-vertex_info_fields = [traction_initial_value,static_coefficient,dynamic_coefficient,slip_weakening_parameter]
+vertex_info_fields = [traction_initial_value,static_coefficient,dynamic_coefficient,slip_weakening_parameter,time_weakening_time,time_weakening_parameter]
 vertex_data_fields = [slip,slip_rate,traction]
 output_freq = time_step
 time_step = 0.04999*s