OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fix actnum with cartesian grid issue.

Browse Source
This commit is contained in:
Liu Ming 2015-10-09 08:43:15 +08:00
parent 28a793c448
commit 8788563088
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/core/grid/MinpvProcessor.hpp
View File

@ -92,7 +92,7 @@ namespace Opm
if (pv.size() != log_size) { if (pv.size() != log_size) {
OPM_THROW(std::runtime_error, "Wrong size of PORV input, must have one element per logical cartesian cell."); OPM_THROW(std::runtime_error, "Wrong size of PORV input, must have one element per logical cartesian cell.");
} }
if (actnum.size() != log_size) { if (!actnum.empty() && actnum.size() != log_size) {
OPM_THROW(std::runtime_error, "Wrong size of ACTNUM input, must have one element per logical cartesian cell."); OPM_THROW(std::runtime_error, "Wrong size of ACTNUM input, must have one element per logical cartesian cell.");
} }
@ -101,7 +101,7 @@ namespace Opm
for (int jj = 0; jj < dims_[1]; ++jj) { for (int jj = 0; jj < dims_[1]; ++jj) {
for (int ii = 0; ii < dims_[0]; ++ii) { for (int ii = 0; ii < dims_[0]; ++ii) {
const int c = ii + dims_[0] * (jj + dims_[1] * kk); const int c = ii + dims_[0] * (jj + dims_[1] * kk);
if (pv[c] < minpv && actnum[c]) { if (pv[c] < minpv && (actnum.empty() || actnum[c])) {
// Move deeper (higher k) coordinates to lower k coordinates. // Move deeper (higher k) coordinates to lower k coordinates.
std::array<double, 8> cz = getCellZcorn(ii, jj, kk, zcorn); std::array<double, 8> cz = getCellZcorn(ii, jj, kk, zcorn);
for (int count = 0; count < 4; ++count) { for (int count = 0; count < 4; ++count) {