From 937b85e88f22e666b7913dd304abc66e3ed380df Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Thu, 19 Jul 2012 14:45:15 +0200
Subject: [PATCH] Bugfix in CompressibleTpfa for gravity case. Gravity
 acceleration should be a factor in the face gravity potential.

---
 opm/core/pressure/CompressibleTpfa.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/opm/core/pressure/CompressibleTpfa.cpp b/opm/core/pressure/CompressibleTpfa.cpp
index 7d936d6e9..154f06336 100644
--- a/opm/core/pressure/CompressibleTpfa.cpp
+++ b/opm/core/pressure/CompressibleTpfa.cpp
@@ -341,7 +341,7 @@ namespace Opm
                         const double depth_diff = face_depth - grid_.cell_centroids[c[j]*dim + dim - 1];
                         props_.density(1, &cell_A_[np*np*c[j]], &gravcontrib[j][0]);
                         for (int p = 0; p < np; ++p) {
-                            gravcontrib[j][p] *= depth_diff;
+                            gravcontrib[j][p] *= depth_diff*grav;
                         }
                     } else {
                         std::fill(gravcontrib[j].begin(), gravcontrib[j].end(), 0.0);