From bbf1862023998f0d968e171cd5b94733a0545a8b Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Fri, 13 Nov 2015 13:30:01 +0100
Subject: [PATCH] Update for new surfaceDensity() interface.

---
 .../FullyImplicitCompressiblePolymerSolver.cpp              | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/opm/polymer/fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp b/opm/polymer/fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp
index 02e553875..54fab8935 100644
--- a/opm/polymer/fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp
+++ b/opm/polymer/fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp
@@ -951,9 +951,9 @@ namespace {
                                                          const std::vector<PhasePresence>& cond,
                                               		     const std::vector<int>&           cells) const
     {
-        const double* rhos = fluid_.surfaceDensity();
-        ADB b = fluidReciprocFVF(phase, p, T, cond, cells);
-        ADB rho = V::Constant(p.size(), 1, rhos[phase]) * b;
+        const V rhos = fluid_.surfaceDensity(phase, cells);
+        const ADB b = fluidReciprocFVF(phase, p, T, cond, cells);
+        const ADB rho = rhos * b;
         return rho;
     }