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
parent 782a5934c9
commit 34f1209b45

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);