From 9ecd153099ec0e62964157e101b0e9469b7aee3d Mon Sep 17 00:00:00 2001
From: Paul Egberts <paul.egberts@tno.nl>
Date: Wed, 14 Oct 2020 17:52:04 +0200
Subject: [PATCH] fix temperature reporting for injectors

---
 opm/simulators/wells/WellState.hpp | 10 +++++++---
 1 file changed, 7 insertions(+), 3 deletions(-)

diff --git a/opm/simulators/wells/WellState.hpp b/opm/simulators/wells/WellState.hpp
index bea26c092..7380a7ad9 100644
--- a/opm/simulators/wells/WellState.hpp
+++ b/opm/simulators/wells/WellState.hpp
@@ -68,13 +68,13 @@ namespace Opm
                 open_for_output_.assign(nw, true);
                 bhp_.resize(nw, 0.0);
                 thp_.resize(nw, 0.0);
-                temperature_.resize(nw, 273.15 + 20); // standard temperature for now
+                temperature_.resize(nw, 273.15 + 15.56); // standard condition temperature
                 wellrates_.resize(nw * np, 0.0);
                 int connpos = 0;
                 for (int w = 0; w < nw; ++w) {
                     const Well& well = wells_ecl[w];
-
-                    // Initialize bhp(), thp(), wellRates().
+              
+                    // Initialize bhp(), thp(), wellRates(), temperature().
                     initSingleWell(cellPressures, w, well, pu, summary_state);
 
                     // Setup wellname -> well index mapping.
@@ -248,6 +248,10 @@ namespace Opm
                 wellrates_[np*w + p] = 0.0;
             }
 
+            if ( well.isInjector() ) { 
+                temperature_[w] = well.injectionControls(summary_state).temperature;
+            }
+
             const int num_perf_this_well = well_perf_data_[w].size();
             if ( num_perf_this_well == 0 ) {
                 // No perforations of the well. Initialize to zero.