OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-01-01 03:37:15 -06:00
Code Issues Packages Projects Releases Wiki Activity

#10411 Surface Extraction: Include inactive cells in active cell area

Browse Source
This commit is contained in:
Magne Sjaastad 2023-06-22 14:11:13 +02:00
parent 0685078ab3
commit f8fd72c823
1 changed files with 9 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
ApplicationLibCode/ProjectDataModel/Surfaces/RimGridCaseSurface.cpp
View File

@ -302,6 +302,12 @@ void RimGridCaseSurface::extractGridDataUsingFourVerticesPerCell()
std::vector<unsigned> triangleIndices;
std::vector<cvf::Vec3d> vertices;
cvf::Vec3st min, max;
if ( activeCells )
{
activeCells->IJKBoundingBox( min, max );
}
for ( size_t i = minI; i <= maxI; i++ )
{
for ( size_t j = minJ; j <= maxJ; j++ )
@ -314,8 +320,9 @@ void RimGridCaseSurface::extractGridDataUsingFourVerticesPerCell()
if ( !m_includeInactiveCells() && activeCells )
{
auto reservoirCellIndex = grid->reservoirCellIndex( currentCellIndex );
if ( !activeCells->isActive( reservoirCellIndex ) ) continue;
if ( i < min.x() || i > max.x() ) continue;
if ( j < min.y() || j > max.y() ) continue;
if ( zeroBasedLayerIndex < min.z() || zeroBasedLayerIndex > max.z() ) continue;
}
cvf::Vec3d currentCornerVerts[8];