From 4e2066e83c3b9a476c02120f9795dbf359982b2e Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Einar=20=C3=96rn=20=C3=93lason?= Date: Tue, 22 Jul 2025 09:58:51 +0200 Subject: [PATCH] Correctly skip land points in the stress update The previous version skipped _ocean_ points and not land points. So I added a not (!) to the commit in PR #884. --- dynamics/src/include/BBMStressUpdateStep.hpp | 2 +- dynamics/src/include/MEVPStressUpdateStep.hpp | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) 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 { #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 { #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.