findSideFaces(): Support grids that don't define an explicit global_cell.

Assume that global_cell==0 denotes the identity mapping (typically a
Cartesian grid).
This commit is contained in:
Bård Skaflestad 2012-03-08 20:53:32 +01:00
parent 52cda7cbb9
commit 501a1d9a51

View File

@ -198,7 +198,8 @@ namespace Opm
const int correct_ijk = (side % 2) ? grid.cartdims[side/2] : 0;
for (int c = 0; c < grid.number_of_cells; ++c) {
int ijk[3] = { -1, -1, -1 };
cartCoord(grid.global_cell[c], grid.cartdims, ijk);
int gc = (grid.global_cell != 0) ? grid.global_cell[c] : c;
cartCoord(gc, grid.cartdims, ijk);
if (ijk[side/2] != correct_ijk) {
continue;
}