From 37677fe03291190ca4ab4f3ba2245885e5e98df9 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Thu, 16 May 2013 14:49:12 +0200
Subject: [PATCH] Print face fluxes after solve.

Also define the HACK_INCOMPRESSIBLE_GRAVITY symbol, since test case
is meaningless without it.
---
 examples/test_impestpfa_ad.cpp | 7 +++++++
 1 file changed, 7 insertions(+)

diff --git a/examples/test_impestpfa_ad.cpp b/examples/test_impestpfa_ad.cpp
index 1b15a0e1a..08343d99b 100644
--- a/examples/test_impestpfa_ad.cpp
+++ b/examples/test_impestpfa_ad.cpp
@@ -20,6 +20,8 @@
 
 #include <config.h>
 
+#define HACK_INCOMPRESSIBLE_GRAVITY 1
+
 #include <opm/autodiff/ImpesTPFAAD.hpp>
 
 #include <opm/core/grid.h>
@@ -140,8 +142,13 @@ main(int argc, char* argv[])
 
     ps.solve(1.0, state, well_state);
 
+    std::cout << "Cell pressure:" << std::endl;
     std::copy(state.pressure().begin(), state.pressure().end(), std::ostream_iterator<double>(std::cout, " "));
     std::cout << std::endl;
+    std::cout << "Face flux:" << std::endl;
+    std::copy(state.faceflux().begin(), state.faceflux().end(), std::ostream_iterator<double>(std::cout, " "));
+    std::cout << std::endl;
+    std::cout << "Well bhp pressure:" << std::endl;
     std::copy(well_state.bhp().begin(), well_state.bhp().end(), std::ostream_iterator<double>(std::cout, " "));
     std::cout << std::endl;