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

Ignore empty branches and invalid cells

Browse Source
This commit is contained in:
Jon Jenssen
2024-10-29 18:02:05 +01:00
committed by jonjenssen
parent 2154c0ae3b
commit 4468e89ded
1 changed files with 9 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
ApplicationLibCode/ReservoirDataModel/RigMswCenterLineCalculator.cpp
View File

@@ -148,8 +148,9 @@ std::vector<SimulationWellCellBranch>
const auto& [gridIndex, cellIndex] = gridAndCellIndices.front();
if ( gridIndex < eclipseCaseData->gridCount() && cellIndex < eclipseCaseData->grid( gridIndex )->cellCount() )
{
const RigCell& cell = eclipseCaseData->grid( gridIndex )->cell( cellIndex );
cvf::Vec3d whStartPos = cell.faceCenter( cvf::StructGridInterface::NEG_K );
const RigCell& cell = eclipseCaseData->grid( gridIndex )->cell( cellIndex );
if ( cell.isInvalid() ) continue;
cvf::Vec3d whStartPos = cell.faceCenter( cvf::StructGridInterface::NEG_K );
// Add extra coordinate between cell face and cell center
// to make sure the well pipe terminated in a segment parallel to z-axis
@@ -209,9 +210,13 @@ std::vector<SimulationWellCellBranch>
}
}
const auto simWellBranch = addCoordsAtCellFaceIntersectionsAndCreateBranch( cellCenterCoords, cellCenterResultPoints, eclipseCaseData );
if ( !cellCenterCoords.empty() )
{
const auto simWellBranch =
addCoordsAtCellFaceIntersectionsAndCreateBranch( cellCenterCoords, cellCenterResultPoints, eclipseCaseData );
simWellBranches.emplace_back( simWellBranch );
simWellBranches.emplace_back( simWellBranch );
}
}
return simWellBranches;