[cig-commits] r20068 - in short/3D/PyLith/branches/v1.7-trunk/pylith: faults feassemble

[brad at geodynamics.org]

Author: brad
Date: 2012-05-10 13:02:09 -0700 (Thu, 10 May 2012)
New Revision: 20068

Modified:
   short/3D/PyLith/branches/v1.7-trunk/pylith/faults/FaultCohesiveDyn.py
   short/3D/PyLith/branches/v1.7-trunk/pylith/feassemble/FIATLagrange.py
   short/3D/PyLith/branches/v1.7-trunk/pylith/feassemble/FIATSimplex.py
Log:
Fixed validation of quadrature dimension (allow 0). Fixed default setting for precsribed tractions with spontaneous rupture.

Modified: short/3D/PyLith/branches/v1.7-trunk/pylith/faults/FaultCohesiveDyn.py
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/pylith/faults/FaultCohesiveDyn.py	2012-05-10 20:01:15 UTC (rev 20067)
+++ short/3D/PyLith/branches/v1.7-trunk/pylith/faults/FaultCohesiveDyn.py	2012-05-10 20:02:09 UTC (rev 20068)
@@ -63,7 +63,7 @@
                                        validator=pyre.inventory.greaterEqual(0.0))
   zeroTolerance.meta['tip'] = "Tolerance for detecting zero values."
 
-  openFreeSurf = pyre.inventory.bool("open_free_surface", default=False)
+  openFreeSurf = pyre.inventory.bool("open_free_surface", default=True)
   openFreeSurf.meta['tip'] = "If True, enforce traction free surface when " \
     "the fault opens, otherwise use initial tractions even when the " \
     "fault opens."
@@ -115,7 +115,6 @@
       self._info.log("Pre-initializing fault '%s'." % self.label())
     FaultCohesive.preinitialize(self, mesh)
     Integrator.preinitialize(self, mesh)
-    self.tract.preinitialize(mesh)
 
     ModuleFaultCohesiveDyn.quadrature(self, self.faultQuadrature)
 
@@ -126,7 +125,7 @@
                                                  "dip_dir"]
 
     if not isinstance(self.tract, NullComponent):
-      print "HERE",self.tract.availableFields
+      self.tract.preinitialize(mesh)
       self.availableFields['vertex']['info'] += self.tract.availableFields['vertex']['info']
 
     self.availableFields['vertex']['info'] += \

Modified: short/3D/PyLith/branches/v1.7-trunk/pylith/feassemble/FIATLagrange.py
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/pylith/feassemble/FIATLagrange.py	2012-05-10 20:01:15 UTC (rev 20067)
+++ short/3D/PyLith/branches/v1.7-trunk/pylith/feassemble/FIATLagrange.py	2012-05-10 20:02:09 UTC (rev 20068)
@@ -31,8 +31,8 @@
 import numpy
 
 def validateDimension(dim):
-  if dim < 1 or dim > 3:
-    raise ValueError("Dimension of Lagrange cell must be 1, 2, or 3.")
+  if dim < 0 or dim > 3:
+    raise ValueError("Dimension of Lagrange cell must be 0, 1, 2, or 3.")
   return dim
 
 

Modified: short/3D/PyLith/branches/v1.7-trunk/pylith/feassemble/FIATSimplex.py
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/pylith/feassemble/FIATSimplex.py	2012-05-10 20:01:15 UTC (rev 20067)
+++ short/3D/PyLith/branches/v1.7-trunk/pylith/feassemble/FIATSimplex.py	2012-05-10 20:02:09 UTC (rev 20068)
@@ -28,8 +28,8 @@
 import numpy
 
 def validateDimension(dim):
-  if dim < 1 or dim > 3:
-    raise ValueError("Dimension of simplex cell must be 1, 2, or 3.")
+  if dim < 0 or dim > 3:
+    raise ValueError("Dimension of simplex cell must be 0, 1, 2, or 3.")
   return dim