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

Ensemble Fracture Statitics: Handle missing P10/P90 area.

Browse Source
This commit is contained in:
Kristian Bendiksen 2021-11-10 12:14:02 +01:00
parent acd36761cf
commit a3e54d5573
1 changed files with 6 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
ApplicationLibCode/ProjectDataModel/Completions/RimEnsembleFractureStatistics.cpp
View File

@ -1213,6 +1213,11 @@ std::shared_ptr<RigSlice2D> RimEnsembleFractureStatistics::setCellsToFillTargetA
const RigSlice2D& distanceGrid, const RigSlice2D& distanceGrid,
double targetArea ) double targetArea )
{ {
std::shared_ptr<RigSlice2D> outputGrid = std::make_shared<RigSlice2D>( grid.nx(), grid.ny() );
// Invalid target area: can happen for P10/P90 grids.
if ( std::isinf( targetArea ) ) return outputGrid;
// Internal cell data class for ordering cells. // Internal cell data class for ordering cells.
class CellData class CellData
{ {
@ -1265,8 +1270,7 @@ std::shared_ptr<RigSlice2D> RimEnsembleFractureStatistics::setCellsToFillTargetA
// Fill cells in the output grid until the target area is reached. // Fill cells in the output grid until the target area is reached.
// This ensures that the statistics fracture grids have representantive sizes. // This ensures that the statistics fracture grids have representantive sizes.
std::shared_ptr<RigSlice2D> outputGrid = std::make_shared<RigSlice2D>( grid.nx(), grid.ny() ); double area = 0.0;
double area = 0.0;
for ( const CellData& cellData : cells ) for ( const CellData& cellData : cells )
{ {
if ( area < targetArea ) if ( area < targetArea )