From 56bcd0c52c8e34dab916cd215928b403630115a0 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Fri, 19 Nov 2010 15:33:52 +0100
Subject: [PATCH] Fixed 'UNSET' bcs for interiour faces. Made total
 compressibilty nonzero for experiment.

---
 dune/porsol/mimetic/TpfaCompressible.hpp | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/dune/porsol/mimetic/TpfaCompressible.hpp b/dune/porsol/mimetic/TpfaCompressible.hpp
index 5d1590e..f697941 100644
--- a/dune/porsol/mimetic/TpfaCompressible.hpp
+++ b/dune/porsol/mimetic/TpfaCompressible.hpp
@@ -183,6 +183,10 @@ namespace Dune
             std::vector<double> bcvalues(num_faces, 0.0);
             for (int face = 0; face < num_faces; ++face) {
                 int bid = pgrid_->grid().boundaryId(face);
+                if (bid == 0) {
+                    bctypes[face] = PressureSolver::FBC_UNSET;
+                    continue;
+                }
                 FlowBC face_bc = bc.flowCond(bid);
                 if (face_bc.isDirichlet()) {
                     bctypes[face] = PressureSolver::FBC_PRESSURE;
@@ -206,7 +210,7 @@ namespace Dune
             const int np = 3;
             const int nc = 3;
             BOOST_STATIC_ASSERT(np == nc);
-            std::vector<double> totcompr(num_cells, 0.0);
+            std::vector<double> totcompr(num_cells, 1.0e-6);
             std::vector<double> voldiscr(num_cells, 0.0);
             std::vector<double> cellA(num_cells*nc*np, 0.0);
             std::vector<double> faceA(num_faces*nc*np, 0.0);