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

Roff: replace -999 values for inactive cells.

Browse Source
This commit is contained in:
Kristian Bendiksen 2023-01-11 10:50:09 +01:00
parent 1d4140e854
commit 903dd34929
1 changed files with 14 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
ApplicationLibCode/FileInterface/RifRoffFileTools.cpp
View File

@ -598,7 +598,20 @@ bool RifRoffFileTools::appendNewInputPropertyResult( RigEclipseCaseData* caseDat
int ny = static_cast<int>( caseData->mainGrid()->cellCountJ() );
int nz = static_cast<int>( caseData->mainGrid()->cellCountK() );
std::vector<double> values = readAndConvertToDouble( nx, ny, nz, keyword, kind, reader );
if ( values.size() != caseData->mainGrid()->cellCount() ) return false;
auto mainGrid = caseData->mainGrid();
if ( values.size() != mainGrid->cellCount() ) return false;
// Set better invalid value for inactive cells: roff file has -999
auto activeCellInfo = caseData->activeCellInfo( RiaDefines::PorosityModelType::MATRIX_MODEL );
size_t cellCount = mainGrid->cellCount();
for ( size_t i = 0; i < cellCount; i++ )
{
if ( !activeCellInfo->isActive( mainGrid->reservoirCellIndex( i ) ) )
{
values[i] = HUGE_VAL;
}
}
RigEclipseResultAddress resAddr( RiaDefines::ResultCatType::INPUT_PROPERTY, resultName );
caseData->results( RiaDefines::PorosityModelType::MATRIX_MODEL )->createResultEntry( resAddr, false );