[cig-commits] r21426 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24

[brad at geodynamics.org]

Author: brad
Date: 2013-03-01 08:42:15 -0800 (Fri, 01 Mar 2013)
New Revision: 21426

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_stresses.py
Log:
Small normal stress at surface.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_stresses.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_stresses.py	2013-03-01 13:23:54 UTC (rev 21425)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_stresses.py	2013-03-01 16:42:15 UTC (rev 21426)
@@ -11,7 +11,7 @@
 # Create spatial variation of initial stresses and fault tractions for
 # TPV24-25.
 
-sim = "tpv25"
+sim = "tpv24"
 
 # ----------------------------------------------------------------------
 # Parameters from benchmark description
@@ -36,7 +36,7 @@
 from spatialdata.spatialdb.SimpleIOAscii import SimpleIOAscii
 from spatialdata.geocoords.CSCart import CSCart
 
-z = numpy.array([0.0, -15.0e+3, -50e+3], dtype=numpy.float64)
+z = numpy.array([0.0, -50.0, -15.0e+3], dtype=numpy.float64)
 points = numpy.zeros( (z.shape[0], 3), dtype=numpy.float64)
 points[:,2] = z
 
@@ -45,10 +45,11 @@
 cs.initialize()
 
 
-mask1 = z >= -15.0e+3
+mask1 = z >= -15.6e+3
 
 Pf = -densityW*gacc*z
 Szz = density*gacc*z
+Szz[z > -1.0] = (density-densityW)*gacc*-10.0
 
 Syy = mask1*(b22*(Szz+Pf)-0*Pf) + ~mask1*Szz
 Sxx = mask1*(b33*(Szz+Pf)-0*Pf) + ~mask1*Szz