diff --git a/ApplicationLibCode/ReservoirDataModel/Completions/RigEclipseToStimPlanCellTransmissibilityCalculator.cpp b/ApplicationLibCode/ReservoirDataModel/Completions/RigEclipseToStimPlanCellTransmissibilityCalculator.cpp index 70ae0e861a..25856faef0 100644 --- a/ApplicationLibCode/ReservoirDataModel/Completions/RigEclipseToStimPlanCellTransmissibilityCalculator.cpp +++ b/ApplicationLibCode/ReservoirDataModel/Completions/RigEclipseToStimPlanCellTransmissibilityCalculator.cpp @@ -217,7 +217,6 @@ void RigEclipseToStimPlanCellTransmissibilityCalculator::calculateStimPlanCellsM } std::vector> polygonsForStimPlanCellInEclipseCell; - cvf::Vec3d areaVector; std::vector stimPlanPolygon = m_stimPlanCell.getPolygon(); for ( const std::vector& planeCellPolygon : planeCellPolygons ) @@ -233,7 +232,6 @@ void RigEclipseToStimPlanCellTransmissibilityCalculator::calculateStimPlanCellsM if ( polygonsForStimPlanCellInEclipseCell.empty() ) continue; std::vector areaOfFractureParts; - double length; std::vector lengthXareaOfFractureParts; double Ax = 0.0; double Ay = 0.0; @@ -241,11 +239,11 @@ void RigEclipseToStimPlanCellTransmissibilityCalculator::calculateStimPlanCellsM for ( const std::vector& fracturePartPolygon : polygonsForStimPlanCellInEclipseCell ) { - areaVector = cvf::GeometryTools::polygonAreaNormal3D( fracturePartPolygon ); - double area = areaVector.length(); + cvf::Vec3d areaVector = cvf::GeometryTools::polygonAreaNormal3D( fracturePartPolygon ); + double area = areaVector.length(); areaOfFractureParts.push_back( area ); - length = RigCellGeometryTools::polygonLengthInLocalXdirWeightedByArea( fracturePartPolygon ); + double length = RigCellGeometryTools::polygonLengthInLocalXdirWeightedByArea( fracturePartPolygon ); lengthXareaOfFractureParts.push_back( length * area ); cvf::Plane fracturePlane; diff --git a/ApplicationLibCode/ReservoirDataModel/Completions/RigFractureTransmissibilityEquations.cpp b/ApplicationLibCode/ReservoirDataModel/Completions/RigFractureTransmissibilityEquations.cpp index f20dc4997a..fc8dd3c510 100644 --- a/ApplicationLibCode/ReservoirDataModel/Completions/RigFractureTransmissibilityEquations.cpp +++ b/ApplicationLibCode/ReservoirDataModel/Completions/RigFractureTransmissibilityEquations.cpp @@ -132,17 +132,15 @@ double RigFractureTransmissibilityEquations::matrixToFractureTrans( double perm, double fractureAreaWeightedlength, double cDarcy ) { - double transmissibility; - double slDivPi = 0.0; if ( cvf::Math::abs( skinfactor ) > EPSILON ) { slDivPi = ( skinfactor * fractureAreaWeightedlength ) / cvf::PI_D; } - transmissibility = 8 * cDarcy * ( perm * NTG ) * A / ( cellSizeLength + slDivPi ); + double transmissibility = 8 * cDarcy * ( perm * NTG ) * A / ( cellSizeLength + slDivPi ); - CVF_ASSERT( transmissibility == transmissibility ); + CVF_ASSERT( !std::isnan( transmissibility ) ); return transmissibility; }