From 2e91d2f353a20d36c1393b880621b3cf28f432a9 Mon Sep 17 00:00:00 2001 From: Tor Harald Sandve Date: Wed, 30 Mar 2022 08:23:41 +0200 Subject: [PATCH] Do not throw when d = 0 instead continue with rs/rv = 0 --- opm/simulators/wells/MultisegmentWellEval.cpp | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/opm/simulators/wells/MultisegmentWellEval.cpp b/opm/simulators/wells/MultisegmentWellEval.cpp index 3f1267f2be..d568783591 100644 --- a/opm/simulators/wells/MultisegmentWellEval.cpp +++ b/opm/simulators/wells/MultisegmentWellEval.cpp @@ -1194,9 +1194,10 @@ getSegmentSurfaceVolume(const EvalWell& temperature, sstr << "Problematic d value " << d << " obtained for well " << baseif_.name() << " during conversion to surface volume with rs " << rs << ", rv " << rv << " and pressure " << seg_pressure - << " obtaining d " << d; + << " obtaining d " << d + << " Continue without dissolution (rs = 0) and vaporization (rv = 0) " + << " for this connection."; OpmLog::debug(sstr.str()); - OPM_THROW_NOLOG(NumericalIssue, sstr.str()); } if (rs > 0.0) { // rs > 0.0?