[cig-commits] r19980 - short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott

[willic3 at geodynamics.org]

Author: willic3
Date: 2012-04-23 20:05:11 -0700 (Mon, 23 Apr 2012)
New Revision: 19980

Modified:
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/plot_profiles.py
Log:
Fixed a couple of things in plotting script.



Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/plot_profiles.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/plot_profiles.py	2012-04-24 01:46:37 UTC (rev 19979)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/plot_profiles.py	2012-04-24 03:05:11 UTC (rev 19980)
@@ -17,7 +17,7 @@
 from pyre.units.time import year
 from pyre.units.length import km,m
 tcycle = 200.0*year
-elastThick = 20.0*km
+elastThick = 40.0*km
 eqslip = 4.0*m
 
 # ======================================================================
@@ -30,6 +30,9 @@
 sys.path.append("../../../figures")
 import matplotlibext
 
+# import pdb
+# pdb.set_trace()
+
 header = 0.45
 
 cells = ["Hex8",
@@ -100,7 +103,13 @@
     dist = vertices[indices,0].squeeze()
     disp = disp[:,indices,1].squeeze()
     
-    indices = numpy.argsort(dist)
+    # Remove value on negative side of fault and sort distances.
+    numSteps = disp.shape[0]
+    dispTest = disp[numSteps -1,:]
+    negInd = numpy.nonzero(dispTest < 0.0)[0][0]
+    distPos = numpy.delete(dist, negInd)
+    indices = numpy.argsort(distPos)
+    
     self.dist = dist[indices] / elastThick.value # Normalize by elastic thickness
     self.disp = disp[:,indices] / eqslip.value # Normalize by eqslip
     self.time = time