OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

WellsManager(): Support grids that do not define "global_cell"

Browse Source 
Assume that the "global_cell" in that case corresponds to the identity
mapping.
This commit is contained in:
Bård Skaflestad 2012-05-15 11:54:59 +02:00
parent c0d43a69bd
commit 43bb6eac78
1 changed files with 10 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
opm/core/WellsManager.cpp
View File

@ -250,9 +250,17 @@ namespace Opm
const int* global_cell = grid.global_cell; const int* global_cell = grid.global_cell;
const int* cpgdim = grid.cartdims; const int* cpgdim = grid.cartdims;
std::map<int,int> cartesian_to_compressed; std::map<int,int> cartesian_to_compressed;
if (global_cell) {
for (int i = 0; i < grid.number_of_cells; ++i) { for (int i = 0; i < grid.number_of_cells; ++i) {
cartesian_to_compressed.insert(std::make_pair(global_cell[i], i)); cartesian_to_compressed.insert(std::make_pair(global_cell[i], i));
} }
}
else {
for (int i = 0; i < grid.number_of_cells; ++i) {
cartesian_to_compressed.insert(std::make_pair(i, i));
}
}
// Get COMPDAT data // Get COMPDAT data
const COMPDAT& compdat = deck.getCOMPDAT(); const COMPDAT& compdat = deck.getCOMPDAT();