diff --git a/src/kinematics_slipplane_opening.f90 b/src/kinematics_slipplane_opening.f90
index 018a5b4b5..a736338f9 100644
--- a/src/kinematics_slipplane_opening.f90
+++ b/src/kinematics_slipplane_opening.f90
@@ -146,30 +146,24 @@ subroutine kinematics_slipplane_opening_LiAndItsTangent(Ld, dLd_dTstar, S, ipc,
 
     udotd = sign(1.0_pReal,traction_d)* prm%sdot0* (  abs(traction_d)/traction_crit &
                                                     - abs(traction_d)/prm%critLoad(i))**prm%n
-    if (abs(udotd) > tol_math_check) then
-      Ld = Ld + udotd*projection_d
-      dudotd_dt = udotd*prm%n/traction_d
-      forall (k=1:3,l=1:3,m=1:3,n=1:3) &
-        dLd_dTstar(k,l,m,n) = dLd_dTstar(k,l,m,n) + dudotd_dt*projection_d(k,l)*projection_d(m,n)
-    endif                
+    Ld = Ld + udotd*projection_d
+    dudotd_dt = udotd*prm%n/traction_d
+    forall (k=1:3,l=1:3,m=1:3,n=1:3) &
+      dLd_dTstar(k,l,m,n) = dLd_dTstar(k,l,m,n) + dudotd_dt*projection_d(k,l)*projection_d(m,n)
   
     udott = sign(1.0_pReal,traction_t)* prm%sdot0* (  abs(traction_t)/traction_crit &
                                                     - abs(traction_t)/prm%critLoad(i))**prm%n
-    if (abs(udott) > tol_math_check) then
-      Ld = Ld + udott*projection_t
-      dudott_dt = udott*prm%n/traction_t
-      forall (k=1:3,l=1:3,m=1:3,n=1:3) &
-        dLd_dTstar(k,l,m,n) = dLd_dTstar(k,l,m,n) + dudott_dt*projection_t(k,l)*projection_t(m,n)
-    endif
+    Ld = Ld + udott*projection_t
+    dudott_dt = udott*prm%n/traction_t
+    forall (k=1:3,l=1:3,m=1:3,n=1:3) &
+      dLd_dTstar(k,l,m,n) = dLd_dTstar(k,l,m,n) + dudott_dt*projection_t(k,l)*projection_t(m,n)
 
     udotn = prm%sdot0* (  max(0.0_pReal,traction_n)/traction_crit &
                         - max(0.0_pReal,traction_n)/prm%critLoad(i))**prm%n
-    if (abs(udotn) > tol_math_check) then
-      Ld = Ld + udotn*projection_n
-      dudotn_dt = udotn*prm%n/traction_n
-      forall (k=1:3,l=1:3,m=1:3,n=1:3) &
-        dLd_dTstar(k,l,m,n) = dLd_dTstar(k,l,m,n) + dudotn_dt*projection_n(k,l)*projection_n(m,n)
-    endif 
+    Ld = Ld + udotn*projection_n
+    dudotn_dt = udotn*prm%n/traction_n
+    forall (k=1:3,l=1:3,m=1:3,n=1:3) &
+      dLd_dTstar(k,l,m,n) = dLd_dTstar(k,l,m,n) + dudotn_dt*projection_n(k,l)*projection_n(m,n)
   enddo
   
   end associate