mirror of
https://github.com/OPM/ResInsight.git
synced 2025-02-25 18:55:39 -06:00
Thermal Fracture: Fix off-by-one in transmissibility calculation.
This commit is contained in:
parent
3dad3c65cb
commit
dbd8dd4ff4
@ -644,7 +644,7 @@ double RimThermalFractureTemplate::resultValueAtIJ( const RigFractureGrid* fract
|
||||
size_t adjustedI = i + 1;
|
||||
size_t adjustedJ = j + 1;
|
||||
|
||||
if ( adjustedI >= fractureGrid->iCellCount() || adjustedJ >= fractureGrid->jCellCount() )
|
||||
if ( adjustedI >= fractureGrid->iCellCount() + 1 || adjustedJ >= fractureGrid->jCellCount() + 1 )
|
||||
{
|
||||
return HUGE_VAL;
|
||||
}
|
||||
|
@ -149,7 +149,7 @@ void RigEclipseToStimPlanCellTransmissibilityCalculator::calculateStimPlanCellsM
|
||||
const std::set<size_t>& reservoirCellIndicesOpenForFlow )
|
||||
{
|
||||
// Not calculating flow into fracture if stimPlan cell cond value is 0 (assumed to be outside the fracture):
|
||||
if ( m_stimPlanCell.getConductivityValue() < 1e-7 ) return;
|
||||
if ( m_stimPlanCell.getConductivityValue() < 1e-7 || std::isinf( m_stimPlanCell.getConductivityValue() ) ) return;
|
||||
|
||||
const RigEclipseCaseData* eclipseCaseData = m_case->eclipseCaseData();
|
||||
|
||||
|
@ -72,7 +72,7 @@ size_t RigFractureCell::getJ() const
|
||||
//--------------------------------------------------------------------------------------------------
|
||||
bool RigFractureCell::hasNonZeroConductivity() const
|
||||
{
|
||||
return m_conductivityValue > 1e-7;
|
||||
return m_conductivityValue > 1e-7 && !std::isinf( m_conductivityValue );
|
||||
}
|
||||
|
||||
//--------------------------------------------------------------------------------------------------
|
||||
|
Loading…
Reference in New Issue
Block a user