OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Bugfix. Fix face direction in the pinchprocessor

Browse Source
This commit is contained in:
Tor Harald Sandve 2015-12-14 10:27:49 +01:00
parent 43abdf5946
commit 8827ca0794
2 changed files with 3 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/core/grid/PinchProcessor.hpp
View File

@ -396,7 +396,7 @@ namespace Opm
}
}
}
pinFaces.push_back(interface_(grid, topCell, Opm::FaceDir::ZMinus));
pinFaces.push_back(interface_(grid, topCell, Opm::FaceDir::ZPlus));
pinCells.push_back(topCell);
tmp.insert(tmp.begin(), topCell);
@ -410,7 +410,7 @@ namespace Opm
}
}
}
pinFaces.push_back(interface_(grid, botCell, Opm::FaceDir::ZPlus));
pinFaces.push_back(interface_(grid, botCell, Opm::FaceDir::ZMinus));
pinCells.push_back(botCell);
}

6
tests/test_pinchprocessor.cpp
View File

@ -112,10 +112,6 @@ BOOST_AUTO_TEST_CASE(Processing)
/ Opm::unit::barsa;
double trans = unit::convert::to(nncdata[0].trans,factor);
std::cout << "WARNING. The opmfil option is hardcoded i.e. the calculated transmissibility ";
std::cout << "is half the correct value due to merging of cells \n";
BOOST_CHECK(std::fabs(trans*2-4.26350022) < 1e-3);
std::cout << trans << std::endl;
std::cout << "WARNING. The opmfil option is hardcoded i.e. the calculated transmissibility is wrong";
//BOOST_CHECK(std::fabs(trans-4.26350022) < 1e-3);
}