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

(#401) Finally got the basic comparison to work and thus Ecl->geomech mapping

Browse Source 
Now weed to clean up, make it fast, include other mappings
Geom/geom Ecl/ecl etc..
and handle inserted/removed layers
This commit is contained in:
Jacob Støren 2015-09-17 15:28:38 +02:00
parent bf84d7d261
commit 3d09ceb92b
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/RigCaseToCaseCellMapper.cpp
View File

@ -402,6 +402,12 @@ bool isEclFemCellsMatching(RigMainGrid* eclGrid, size_t reservoirCellIndex, Rig
eclShallowQuad[3] = gomConvertedEclCell[localElmNodeIndicesForBotZFace[3]];
}
if (!eclGrid->isFaceNormalsOutwards())
{
flipQuadWinding(femShallowQuad);
flipQuadWinding(femDeepestQuad);
}
// Now the top/bottom have opposite winding. To make the comparisons and index rotations simpler
// flip the winding of the bottom face:
@ -488,8 +494,8 @@ RigCaseToCaseCellMapper::RigCaseToCaseCellMapper(RigMainGrid* masterEclGrid, Rig
double cellSizeI, cellSizeJ, cellSizeK;
masterEclGrid->characteristicCellSizes(&cellSizeI, &cellSizeJ, &cellSizeK);
double xyTolerance = cellSizeI* 2;
double zTolerance = cellSizeK* 2;
double xyTolerance = cellSizeI* 0.1;
double zTolerance = cellSizeK* 0.1;
int elementCount = dependentFemPart->elementCount();
cvf::Vec3d elmCorners[8];