From 0ac0c97612e620fd4039bbae8c02ff539f06172f Mon Sep 17 00:00:00 2001 From: Anders Matheson Date: Mon, 19 Sep 2016 16:43:41 +0200 Subject: [PATCH] Output rates for present phases --- opm/autodiff/WellStateFullyImplicitBlackoil.hpp | 9 +++++++-- 1 file changed, 7 insertions(+), 2 deletions(-) diff --git a/opm/autodiff/WellStateFullyImplicitBlackoil.hpp b/opm/autodiff/WellStateFullyImplicitBlackoil.hpp index ac7fbf428..c1310c10b 100644 --- a/opm/autodiff/WellStateFullyImplicitBlackoil.hpp +++ b/opm/autodiff/WellStateFullyImplicitBlackoil.hpp @@ -230,10 +230,15 @@ namespace Opm const auto& wv = this->wellRates(); data::Rates wellrates; - if( np == 3 ) { - /* only write if 3-phase solution */ + if( np > 0 ) { wellrates.set( rt::wat, wv[ wellrate_index + 0 ] ); + } + + if( np > 1 ) { wellrates.set( rt::oil, wv[ wellrate_index + 1 ] ); + } + + if( np > 2 ) { wellrates.set( rt::gas, wv[ wellrate_index + 2 ] ); }