[cig-commits] r16189 - short/3D/PyLith/trunk/libsrc/faults

[brad at geodynamics.org]

Author: brad
Date: 2010-01-27 17:19:29 -0800 (Wed, 27 Jan 2010)
New Revision: 16189

Modified:
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDynL.cc
Log:
Started setting up calculating slip rate.

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDynL.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDynL.cc	2010-01-28 00:14:37 UTC (rev 16188)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDynL.cc	2010-01-28 01:19:29 UTC (rev 16189)
@@ -656,6 +656,15 @@
   const ALE::Obj<RealSection>& slipSection =  _fields->get("slip").section();
   assert(!slipSection.isNull());
 
+  // :TODO: Get slip rate from velocity(t)
+  //_updateSlipRate(*fields);
+  double_array slipRateVertex(spaceDim);
+#if 0
+  const ALE::Obj<RealSection>& slipRateSection =
+    _fields->get("slip_rate").section();
+  assert(!slipRateSection.isNull());
+#endif
+
   const ALE::Obj<RealSection>& areaSection =  _fields->get("area").section();
   assert(!areaSection.isNull());
 
@@ -713,6 +722,12 @@
       slipSection->restrictPoint(vertexFault, 
 				 &slipVertex[0], slipVertex.size());
 
+#if 0
+      // Get slip rate
+      slipRateSection->restrictPoint(vertexFault, 
+				     &slipRateVertex[0], slipRateVertex.size());
+#endif
+
       // Get total fault area asssociated with vertex (assembled over all cells)
       const double* areaVertex = areaSection->restrictPoint(vertexFault);
       assert(0 != areaVertex);