From af76b390d2f03505452f555cccc8297fc6b517d6 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= Date: Mon, 21 May 2012 10:19:12 +0200 Subject: [PATCH] Added initialization of face pressures. --- opm/core/utility/initState_impl.hpp | 9 +++++++++ 1 file changed, 9 insertions(+) diff --git a/opm/core/utility/initState_impl.hpp b/opm/core/utility/initState_impl.hpp index 6a7d1202..389d5d40 100644 --- a/opm/core/utility/initState_impl.hpp +++ b/opm/core/utility/initState_impl.hpp @@ -410,6 +410,9 @@ namespace Opm const double ref_z = grid.cell_centroids[0 + grid.dimensions - 1]; initHydrostaticPressure(grid, dens, ref_z, gravity, ref_z, ref_p, state); } + + // Finally, init face pressures. + initFacePressure(grid, state); } @@ -488,6 +491,9 @@ namespace Opm const double woc = -1e100; initHydrostaticPressure(grid, props, woc, gravity, ref_z, ref_p, state); } + + // Finally, init face pressures. + initFacePressure(grid, state); } @@ -538,6 +544,9 @@ namespace Opm } else { THROW("initStateFromDeck(): we must either have EQUIL, or both SWAT and PRESSURE."); } + + // Finally, init face pressures. + initFacePressure(grid, state); }