From abb8a63e6d6b34da69b9d4727838a794ba72bfc7 Mon Sep 17 00:00:00 2001 From: osae Date: Fri, 28 Mar 2014 17:18:50 +0100 Subject: [PATCH] Default equil region should be one not zero ... Otherwise problems when kw EQLNUM is used. --- opm/core/simulator/initStateEquil.hpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/opm/core/simulator/initStateEquil.hpp b/opm/core/simulator/initStateEquil.hpp index 566f2b00f..cfcb56bae 100644 --- a/opm/core/simulator/initStateEquil.hpp +++ b/opm/core/simulator/initStateEquil.hpp @@ -241,8 +241,8 @@ namespace Opm } else { // No explicit equilibration region. - // All cells in region zero. - eqlnum.assign(G.number_of_cells, 0); + // All cells in region one. + eqlnum.assign(G.number_of_cells, 1); } return eqlnum; @@ -384,7 +384,7 @@ namespace Opm r = 0, nr = reg.numRegions(); r < nr; ++r) { - const typename RMap::CellRange cells = reg.cells(r); + const typename RMap::CellRange cells = reg.cells(r+1); const int repcell = *cells.begin(); const RhoCalc calc(props, repcell);