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

#4874 Dual Porosity : Missing active cells in some cases

Browse Source
This commit is contained in:
Magne Sjaastad 2019-10-16 11:11:01 +02:00
parent a53f1e1e91
commit 61db98ef04
1 changed files with 10 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
ApplicationCode/FileInterface/RifActiveCellsReader.cpp
View File

@ -89,19 +89,22 @@ std::vector<std::vector<int>> RifActiveCellsReader::activeCellsFromPorvKeyword(
if ( porvValues[poreValueIndex] > 0.0 )
{
if ( !dualPorosity || poreValueIndex < porvValues.size() / 2 )
if ( dualPorosity )
{
activeCellsOneGrid[indexToCell] = CELL_ACTIVE_MATRIX;
if ( poreValueIndex < activeCellCount )
{
activeCellsOneGrid[indexToCell] += CELL_ACTIVE_MATRIX;
}
else
{
activeCellsOneGrid[indexToCell] += CELL_ACTIVE_FRACTURE;
}
}
else
{
activeCellsOneGrid[indexToCell] += CELL_ACTIVE_FRACTURE;
activeCellsOneGrid[indexToCell] = CELL_ACTIVE_MATRIX;
}
}
else
{
activeCellsOneGrid[indexToCell] = 0;
}
}
activeCellsAllGrids.push_back( activeCellsOneGrid );