Fix for setting cumulative slip across earthquake ruptures (first set…

… to zero).
geodynamics · Jun 11, 2019 · 474bfdb · 474bfdb
1 parent 7bc6dff
commit 474bfdb
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/libsrc/pylith/faults/FaultCohesiveKin.cc b/libsrc/pylith/faults/FaultCohesiveKin.cc
@@ -258,6 +258,9 @@ pylith::faults::FaultCohesiveKin::updateAuxiliaryField(pylith::topology::Field*
     PYLITH_METHOD_BEGIN;
     PYLITH_COMPONENT_DEBUG("updateAuxiliaryField(auxiliaryField="<<auxiliaryField<<", t="<<t<<")");
 
+
+    auxiliaryField->zeroLocal(); // :KLUDGE: :TEMPORARY: Slip field is only auxiliary field, so we can zero entire field.
+
     // Compute slip field at current time step
     const srcs_type::const_iterator rupturesEnd = _ruptures.end();
     for (srcs_type::iterator r_iter = _ruptures.begin(); r_iter != rupturesEnd; ++r_iter) {