Default equil region should be one not zero ...

Otherwise problems when kw EQLNUM is used.
This commit is contained in:
osae 2014-03-28 17:18:50 +01:00 committed by Andreas Lauser
parent a2cd03197c
commit abb8a63e6d

View File

@ -241,8 +241,8 @@ namespace Opm
} }
else { else {
// No explicit equilibration region. // No explicit equilibration region.
// All cells in region zero. // All cells in region one.
eqlnum.assign(G.number_of_cells, 0); eqlnum.assign(G.number_of_cells, 1);
} }
return eqlnum; return eqlnum;
@ -384,7 +384,7 @@ namespace Opm
r = 0, nr = reg.numRegions(); r = 0, nr = reg.numRegions();
r < nr; ++r) 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 int repcell = *cells.begin();
const RhoCalc calc(props, repcell); const RhoCalc calc(props, repcell);