From 7c3020822a6b577f45d4981bb2eb7fd376ffeca8 Mon Sep 17 00:00:00 2001
From: Andreas Lauser <and@poware.org>
Date: Mon, 22 Feb 2016 15:14:09 +0100
Subject: [PATCH] properly create state.temperature in the
 FullyImplicitCompressiblePolymerSolver class

---
 .../fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp  | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/opm/polymer/fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp b/opm/polymer/fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp
index 93d7ab625..1ca4df5dd 100644
--- a/opm/polymer/fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp
+++ b/opm/polymer/fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp
@@ -440,6 +440,10 @@ namespace {
         int nextvar = 0;
         state.pressure = vars[ nextvar++ ];
 
+        // Temperature. (this is always a constant so far)
+        const V T = Eigen::Map<const V>(& x.temperature()[0], nc, 1);
+        state.temperature = ADB::constant(T);
+
         // Saturation.
         const std::vector<int>& bpat = vars[0].blockPattern();
         {