[cig-commits] r19990 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210

[brad at geodynamics.org]

Author: brad
Date: 2012-04-24 15:11:37 -0700 (Tue, 24 Apr 2012)
New Revision: 19990

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_ruptime.py
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_sliprate.py
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_velth.py
Log:
Added multipart figure labels.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_ruptime.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_ruptime.py	2012-04-24 21:40:10 UTC (rev 19989)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_ruptime.py	2012-04-24 22:11:37 UTC (rev 19990)
@@ -11,6 +11,7 @@
 
 sim = "tpv13"
 griddata = False
+showAB = True
 
 # ----------------------------------------------------------------------
 import tables
@@ -29,6 +30,8 @@
              ("black", 'solid'),
              ]
 
+labelsAB = 'ab'
+
 # ----------------------------------------------------------------------
 def getval(v):
     try:
@@ -120,6 +123,11 @@
             ax.set_yticklabels([])
             ax.set_ylabel("")
 
+        if showAB:
+            ilabel = icol
+            ax.text(-15, -1, "(%s)" % labelsAB[ilabel],
+                     fontweight='bold')
+
 pyplot.show()
 pyplot.savefig("%s_ruptime" % (sim))
 

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_sliprate.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_sliprate.py	2012-04-24 21:40:10 UTC (rev 19989)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_sliprate.py	2012-04-24 22:11:37 UTC (rev 19990)
@@ -10,6 +10,7 @@
 # PREREQUISITES: matplotlib, numpy
 
 sim = "tpv13"
+showAB = True
 
 # ----------------------------------------------------------------------
 import tables
@@ -29,6 +30,8 @@
              ("black", (None, None)),
              ]
 
+labelsAB = 'abcdef'
+
 # ----------------------------------------------------------------------
 def getval(v):
     try:
@@ -66,7 +69,7 @@
     
         isim = 0
         for (label, simdir) in simdirs:
-            filename = "%s/faultst+%03ddp%03d.dat" % \
+            filename = "%s/faultst%+04ddp%03d.dat" % \
                 (simdir, int(locs[iloc][0]*10), int(locs[iloc][1]*10))
             data = numpy.loadtxt(filename, comments="#", usecols=(0,5),
                                  converters={0: getval,
@@ -106,6 +109,11 @@
         if icol > 0:
             ax.set_yticklabels([])
             ax.set_ylabel("")
+
+        if showAB:
+            ilabel = irow*ncols + icol
+            ax.text(0.15, 14.5, "(%s)" % labelsAB[ilabel],
+                    fontweight='bold')
     
         iloc += 1
 

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_velth.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_velth.py	2012-04-24 21:40:10 UTC (rev 19989)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_velth.py	2012-04-24 22:11:37 UTC (rev 19990)
@@ -10,6 +10,7 @@
 # PREREQUISITES: matplotlib, numpy
 
 sim = "tpv13"
+showAB = True
 
 # ----------------------------------------------------------------------
 import tables
@@ -29,6 +30,8 @@
              ("black", (None, None)),
              ]
 
+labelsAB = 'abcd'
+
 # ----------------------------------------------------------------------
 def getval(v):
     try:
@@ -111,6 +114,11 @@
                 ax.set_yticklabels([])
                 ax.set_ylabel("")
 
+            if showAB:
+                ilabel = irow*ncols + icol
+                ax.text(0.15, 2.4, "(%s)" % labelsAB[ilabel],
+                        fontweight='bold')
+
         isim += 1
     
 pyplot.show()