[cig-commits] r19733 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d

[brad at geodynamics.org]

Author: brad
Date: 2012-03-06 14:17:56 -0800 (Tue, 06 Mar 2012)
New Revision: 19733

Added:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/plot_tpv13.py
Log:
Added diagnostic plotting routine.

Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/plot_tpv13.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/plot_tpv13.py	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/plot_tpv13.py	2012-03-06 22:17:56 UTC (rev 19733)
@@ -0,0 +1,43 @@
+#!/usr/bin/env python
+#
+# ======================================================================
+#
+#                           Brad T. Aagaard
+#                        U.S. Geological Survey
+#
+# ======================================================================
+#
+
+cell = "tri3"
+dx = 200
+
+sim = "tpv13"
+inputRoot = "output/%s_%s_%03dm-fault" % (sim,cell,dx)
+
+# ----------------------------------------------------------------------
+import tables
+import numpy
+import pylab
+
+# ----------------------------------------------------------------------
+h5 = tables.openFile("%s.h5" % inputRoot, 'r')
+vertices = h5.root.geometry.vertices[:]
+traction = h5.root.vertex_fields.traction[1,:,:].squeeze()
+h5.close()
+
+from math import pi, sin
+dipDist = vertices[:,1]/sin(60.0*pi/180.0)
+
+# Expected tractions
+if sim == "tpv13":
+    mask = dipDist > -13.8e+3
+    normalE = mask*7390.01*dipDist + ~mask*14427.98*dipDist
+    shearE = mask*-0.549847*normalE + ~mask*0.0
+
+pylab.plot(shearE, dipDist, 'k*',
+           traction[:,0], dipDist, 'rx',
+           normalE, dipDist, 'k*',
+           traction[:,1], dipDist, 'rx')
+pylab.show()
+
+# End of file


Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/plot_tpv13.py
___________________________________________________________________
Name: svn:executable
   + *