Fix cohesion bug by adding epsilon1 and epsilon2 (#707)

underworldcode · Oct 31, 2024 · 3f17aba · 3f17aba
1 parent 8ed37d7
commit 3f17aba
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/underworld/UWGeodynamics/_rheology.py b/underworld/UWGeodynamics/_rheology.py
@@ -304,7 +304,9 @@ def _cohesionFn(self):
             cohesion = self.cohesionWeakeningFn(
                 self.plasticStrain,
                 Cohesion=nd(self.cohesion),
-                CohesionSw=nd(self.cohesionAfterSoftening))
+                CohesionSw=nd(self.cohesionAfterSoftening),
+                epsilon1=self.epsilon1,
+                epsilon2=self.epsilon2)
         else:
             cohesion = fn.misc.constant(self.cohesion)
         return cohesion