diff --git a/geolib_plus/robertson_cpt_interpretation/robertson_cpt_interpretation.py b/geolib_plus/robertson_cpt_interpretation/robertson_cpt_interpretation.py
index c9d56de..103d146 100644
--- a/geolib_plus/robertson_cpt_interpretation/robertson_cpt_interpretation.py
+++ b/geolib_plus/robertson_cpt_interpretation/robertson_cpt_interpretation.py
@@ -386,7 +386,7 @@ def stress_calc(self):
         z_aux = np.min(
             [self.data.pwp, abs(self.data.depth_to_reference[0]) + self.data.depth[0]]
         )
-        pwp = (z_aux - abs(self.data.depth_to_reference)) * self.data.g
+        pwp = (z_aux - self.data.depth_to_reference) * self.data.g
         # no suction is allowed
         pwp[pwp <= 0] = 0
         # compute effective stress