From d8830c721b48209a50c7243b7a55d8002c89f73e Mon Sep 17 00:00:00 2001
From: Joakim Hove <joakim.hove@gmail.com>
Date: Sun, 28 Feb 2021 15:48:02 +0100
Subject: [PATCH] Use only active elements when updating reference porosity

---
 ebos/eclproblem.hh | 5 ++---
 1 file changed, 2 insertions(+), 3 deletions(-)

diff --git a/ebos/eclproblem.hh b/ebos/eclproblem.hh
index eba3a2d3e..43ab9a41a 100644
--- a/ebos/eclproblem.hh
+++ b/ebos/eclproblem.hh
@@ -2669,11 +2669,10 @@ private:
         referencePorosity_[/*timeIdx=*/0].resize(numDof);
 
         const auto& fp = eclState.fieldProps();
-        const std::vector<double> porvData = fp.porv(true);
+        const std::vector<double> porvData = fp.porv(false);
         const std::vector<int> actnumData = fp.actnum();
         for (size_t dofIdx = 0; dofIdx < numDof; ++ dofIdx) {
-            unsigned cartElemIdx = vanguard.cartesianIndex(dofIdx);
-            Scalar poreVolume = porvData[cartElemIdx];
+            Scalar poreVolume = porvData[dofIdx];
 
             // we define the porosity as the accumulated pore volume divided by the
             // geometric volume of the element. Note that -- in pathetic cases -- it can