From ace2930603f9fdfb683e7b4f2faf487d4faf5921 Mon Sep 17 00:00:00 2001 From: Paul Egberts Date: Mon, 22 Feb 2021 14:23:07 +0100 Subject: [PATCH] use of normalizeKrValues() to account for TOLCRIT --- .../fluidmatrixinteractions/EclMaterialLawManager.hpp | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/opm/material/fluidmatrixinteractions/EclMaterialLawManager.hpp b/opm/material/fluidmatrixinteractions/EclMaterialLawManager.hpp index 4b1aa1a8c..66753571e 100644 --- a/opm/material/fluidmatrixinteractions/EclMaterialLawManager.hpp +++ b/opm/material/fluidmatrixinteractions/EclMaterialLawManager.hpp @@ -960,6 +960,9 @@ private: auto& effParams = *dest[satRegionIdx]; + const auto tolcrit = eclState.runspec().saturationFunctionControls() + .minimumRelpermMobilityThreshold(); + const auto& tableManager = eclState.getTableManager(); switch (eclState.runspec().saturationFunctionControls().family()) { @@ -977,10 +980,10 @@ private: std::vector SwColumn = swfnTable.getColumn("SW").vectorCopy(); - effParams.setKrwSamples(SwColumn, swfnTable.getColumn("KRW").vectorCopy()); - effParams.setKrnSamples(SwColumn, sgfnTable.getColumn("KRG").vectorCopy()); + effParams.setKrwSamples(SwColumn, normalizeKrValues_(tolcrit, swfnTable.getColumn("KRW"))); + effParams.setKrnSamples(SwColumn, normalizeKrValues_(tolcrit, sgfnTable.getColumn("KRG"))); //Capillary pressure is read from SWFN. - //For gas-water system the capillary pressure column valeas are set to 0 in SGFN + //For gas-water system the capillary pressure column values are set to 0 in SGFN effParams.setPcnwSamples(SwColumn, swfnTable.getColumn("PCOW").vectorCopy()); effParams.finalize();