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

Zero-porosity cells and MINPV.

Browse Source 
Model 2 partly uses PORO=0.0 to signal inactive cells. This PR
prevents such cells from causing unwanted zcorn modifications in
the current min-pv algoritm. (The problem is however more general as
zero PORO or NTG should kill a cell whether MINPV is active or not,
refer kw ACTNUM.)
This commit is contained in:
osae 2015-08-27 17:49:10 +02:00
parent 6579a150f3
commit f81a4d62a4
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -109,7 +109,7 @@ namespace Opm
} }
setCellZcorn(ii, jj, kk, cz, zcorn); setCellZcorn(ii, jj, kk, cz, zcorn);
// Check if there is a cell below. // Check if there is a cell below.
if (kk < dims_[2] - 1) { if (pv[c] > 0.0 && kk < dims_[2] - 1) {
// Set lower k coordinates of cell below to upper cells's coordinates. // Set lower k coordinates of cell below to upper cells's coordinates.
std::array<double, 8> cz_below = getCellZcorn(ii, jj, kk + 1, zcorn); std::array<double, 8> cz_below = getCellZcorn(ii, jj, kk + 1, zcorn);
for (int count = 0; count < 4; ++count) { for (int count = 0; count < 4; ++count) {