From ab1b57b994d5cd5c920666ebc9caa1cced5d1c64 Mon Sep 17 00:00:00 2001 From: Kai Bao Date: Mon, 31 Jul 2017 12:43:50 +0200 Subject: [PATCH] fixing the if condition for 2p + dummy phase from if (numComp < numEq) to if (numComp < numWellEq) --- opm/autodiff/StandardWellsDense_impl.hpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/opm/autodiff/StandardWellsDense_impl.hpp b/opm/autodiff/StandardWellsDense_impl.hpp index c9ba026a3..39100538a 100644 --- a/opm/autodiff/StandardWellsDense_impl.hpp +++ b/opm/autodiff/StandardWellsDense_impl.hpp @@ -238,8 +238,8 @@ namespace Opm { } } - // add trivial equation for 2p cases (Only support water + oil) - if (numComp < numEq ) { + // add a trivial equation for the dummy phase for 2p cases (Only support water + oil) + if ( numComp < numWellEq ) { assert(!active_[ Gas ]); invDuneD_[w][w][Gas][Gas] = 1.0; }