From ebcaea3faa5f7a6c918e31bcf270c4e16491fb96 Mon Sep 17 00:00:00 2001
From: Tor Harald Sandve <tosa@norceresearch.no>
Date: Wed, 14 Apr 2021 15:58:58 +0200
Subject: [PATCH] get bhp from thp without solving

---
 opm/simulators/wells/StandardWell_impl.hpp | 15 +++++----------
 1 file changed, 5 insertions(+), 10 deletions(-)

diff --git a/opm/simulators/wells/StandardWell_impl.hpp b/opm/simulators/wells/StandardWell_impl.hpp
index b96956b09..6099fecb5 100644
--- a/opm/simulators/wells/StandardWell_impl.hpp
+++ b/opm/simulators/wells/StandardWell_impl.hpp
@@ -1547,17 +1547,12 @@ namespace Opm
             }
             case Well::ProducerCMode::THP:
             {
-                well_state.thp()[well_index] = this->getTHPConstraint(summaryState);
-                gliftDebug(
-                    "computing BHP from THP to update well state",
-                    deferred_logger);
-                auto bhp = computeBhpAtThpLimitProd(well_state, ebos_simulator, summaryState, deferred_logger);
-                if (bhp) {
-                    well_state.bhp()[well_index] = *bhp;
-                } else {
-                    deferred_logger.warning("FAILURE_GETTING_CONVERGED_BHP",
-                                            "Failed to find BHP when switching to THP control for well " + name());
+                std::vector<double> rates(3, 0.0);
+                for (int p = 0; p<np; ++p) {
+                    rates[p] = well_state.wellRates()[well_index*np + p];
                 }
+                double bhp = calculateBhpFromThp(well_state, rates, well, summaryState, deferred_logger);
+                well_state.bhp()[well_index] = bhp;
                 break;
             }
             case Well::ProducerCMode::GRUP: