diff --git a/dynamics/src/include/BBMStressUpdateStep.hpp b/dynamics/src/include/BBMStressUpdateStep.hpp
index 188488de6..93cb47f7c 100644
--- a/dynamics/src/include/BBMStressUpdateStep.hpp
+++ b/dynamics/src/include/BBMStressUpdateStep.hpp
@@ -53,7 +53,7 @@ class BBMStressUpdateStep : public StressUpdateStep<DGadvection, DGstress> {
 #pragma omp parallel for
         for (size_t i = 0; i < smesh.nelements; ++i) {
 
-            if (smesh.landmask[i])
+            if (!smesh.landmask[i])
                 continue;
 
             //! Evaluate values in Gauss points (3 point Gauss rule in 2d => 9 points)
diff --git a/dynamics/src/include/MEVPStressUpdateStep.hpp b/dynamics/src/include/MEVPStressUpdateStep.hpp
index 58e3868e1..d88eed95e 100644
--- a/dynamics/src/include/MEVPStressUpdateStep.hpp
+++ b/dynamics/src/include/MEVPStressUpdateStep.hpp
@@ -49,7 +49,7 @@ class MEVPStressUpdateStep : public StressUpdateStep<DGadvection, DGstress> {
 #pragma omp parallel for
         for (size_t i = 0; i < smesh.nelements; ++i) {
 
-            if (smesh.landmask[i])
+            if (!smesh.landmask[i])
                 continue;
 
             // Here, one should check if it is enough to use a 2-point Gauss rule.