From 3f17aba92fe35bd6b37b8420849d1c6a29c9638d Mon Sep 17 00:00:00 2001
From: tiannh7 <156781848+tiannh7@users.noreply.github.com>
Date: Thu, 31 Oct 2024 10:52:27 +0800
Subject: [PATCH] Fix cohesion bug by adding epsilon1 and epsilon2 (#707)

---
 underworld/UWGeodynamics/_rheology.py | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/underworld/UWGeodynamics/_rheology.py b/underworld/UWGeodynamics/_rheology.py
index 55f99283e..3bbe87d03 100644
--- 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