From 867b769fcbbc8ddebc5908a47ad68471857aa3b1 Mon Sep 17 00:00:00 2001
From: Kai Bao <kai.bao@sintef.no>
Date: Fri, 5 Feb 2021 11:22:49 +0100
Subject: [PATCH] only handling the interior elements for numerical aquifer

---
 opm/simulators/aquifers/AquiferNumerical.hpp | 7 +++++++
 1 file changed, 7 insertions(+)

diff --git a/opm/simulators/aquifers/AquiferNumerical.hpp b/opm/simulators/aquifers/AquiferNumerical.hpp
index 9782ea86c..a7b14e069 100644
--- a/opm/simulators/aquifers/AquiferNumerical.hpp
+++ b/opm/simulators/aquifers/AquiferNumerical.hpp
@@ -128,6 +128,9 @@ private:
         const auto& elemEndIt = gridView.template end</*codim=*/0>();
         for (; elemIt != elemEndIt; ++elemIt) {
             const auto& elem = *elemIt;
+            if (elem.partitionType() != Dune::InteriorEntity) {
+                continue;
+            }
             elem_ctx.updatePrimaryStencil(elem);
 
             const size_t cell_index = elem_ctx.globalSpaceIndex(/*spaceIdx=*/0, /*timeIdx=*/0);
@@ -135,6 +138,7 @@ private:
             if (idx < 0) {
                 continue;
             }
+
             elem_ctx.updatePrimaryIntensiveQuantities(/*timeIdx=*/0);
             const auto& iq0 = elem_ctx.intensiveQuantities(/*spaceIdx=*/0, /*timeIdx=*/0);
             const auto& fs = iq0.fluidState();
@@ -168,6 +172,9 @@ private:
         const auto& elemEndIt = gridView.template end</*codim=*/0>();
         for (; elemIt != elemEndIt; ++elemIt) {
             const auto &elem = *elemIt;
+            if (elem.partitionType() != Dune::InteriorEntity) {
+                continue;
+            }
             // elem_ctx.updatePrimaryStencil(elem);
             elem_ctx.updateStencil(elem);