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

#557, #560, #561) WLP: Allign WL-extraction code Geom-Eclipse

Browse Source
This commit is contained in:
Jacob Støren 2015-10-14 09:39:03 +02:00
parent 4c73b31919
commit 8d2fee9bf4
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
ApplicationCode/ReservoirDataModel/RigEclipseWellLogExtractor.cpp
View File

@ -84,6 +84,14 @@ void RigEclipseWellLogExtractor::calculateIntersection()
int intersectionCount = RigHexIntersector::lineHexCellIntersection(p1, p2, hexCorners, closeCells[cIdx], &intersections); int intersectionCount = RigHexIntersector::lineHexCellIntersection(p1, p2, hexCorners, closeCells[cIdx], &intersections);
} }
if (!isCellFaceNormalsOut)
{
for (size_t intIdx = 0; intIdx < intersections.size(); ++intIdx)
{
intersections[intIdx].m_isIntersectionEntering = !intersections[intIdx].m_isIntersectionEntering ;
}
}
// Now, with all the intersections of this piece of line, we need to // Now, with all the intersections of this piece of line, we need to
// sort them in order, and set the measured depth and corresponding cell index // sort them in order, and set the measured depth and corresponding cell index
@ -95,8 +103,6 @@ void RigEclipseWellLogExtractor::calculateIntersection()
for (size_t intIdx = 0; intIdx < intersections.size(); ++intIdx) for (size_t intIdx = 0; intIdx < intersections.size(); ++intIdx)
{ {
if (!isCellFaceNormalsOut) intersections[intIdx].m_isIntersectionEntering = !intersections[intIdx].m_isIntersectionEntering ;
double lenghtAlongLineSegment1 = (intersections[intIdx].m_intersectionPoint - p1).length(); double lenghtAlongLineSegment1 = (intersections[intIdx].m_intersectionPoint - p1).length();
double lenghtAlongLineSegment2 = (p2 - intersections[intIdx].m_intersectionPoint).length(); double lenghtAlongLineSegment2 = (p2 - intersections[intIdx].m_intersectionPoint).length();
double measuredDepthDiff = md2 - md1; double measuredDepthDiff = md2 - md1;