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(& x.temperature()[0], nc, 1); + state.temperature = ADB::constant(T); + // Saturation. const std::vector& bpat = vars[0].blockPattern(); {