From c7950ff3417f9f89f89f469db038d0a1bd7f91af Mon Sep 17 00:00:00 2001
From: Ove Saevareid <osae@iris.no>
Date: Fri, 11 Feb 2011 11:32:52 +0100
Subject: [PATCH] Correcting (hopefully:)) compressibility computations.

---
 dune/porsol/blackoil/fluid/FluidSystemBlackoil.hpp | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/dune/porsol/blackoil/fluid/FluidSystemBlackoil.hpp b/dune/porsol/blackoil/fluid/FluidSystemBlackoil.hpp
index 05a4be5e..87f24a8b 100644
--- a/dune/porsol/blackoil/fluid/FluidSystemBlackoil.hpp
+++ b/dune/porsol/blackoil/fluid/FluidSystemBlackoil.hpp
@@ -217,8 +217,11 @@ public:
         double data_for_C[numComponents*numPhases];
         Dune::SharedFortranMatrix C(numComponents, numPhases, data_for_C);
         Dune::prod(Ai, dA, C);
-        CompVec ones(1.0);
-        cp = Dune::prod(C, ones);
+        //CompVec ones(1.0);
+        //cp = Dune::prod(C, ones); // Probably C' and not C; we want phasewise compressibilities:
+        cp[Aqua] = C(Water, Aqua);
+        cp[Liquid] = C(Oil, Liquid) + C(Gas, Liquid);
+        cp[Vapour] = C(Gas, Vapour) + C(Oil, Vapour);  
         fluid_state.total_compressibility_ = cp*s;
     }