OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-18 21:43:27 -06:00
Code Issues Packages Projects Releases Wiki Activity

Initialisation of EGLNUM (analogous to e.g SATNUM)

Browse Source
This commit is contained in:
osae 2014-05-06 12:47:24 +02:00
parent ca5b0eb16f
commit d648917fcf
1 changed files with 8 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
opm/core/simulator/initStateEquil.hpp
View File

@ -311,24 +311,15 @@ namespace Opm
const UnstructuredGrid& G ) const UnstructuredGrid& G )
{ {
std::vector<int> eqlnum; std::vector<int> eqlnum;
if (newParserDeck->hasKeyword("EQLNUM")) { if (newParserDeck->hasKeyword("EQLNUM")) {
const std::vector<int>& e = newParserDeck->getKeyword("EQLNUM")->getIntData(); eqlnum.resize(G.number_of_cells);
eqlnum.reserve(e.size()); const std::vector<int>& e =
if (newParserDeck->hasKeyword("ACTNUM")) { newParserDeck->getKeyword("EQLNUM")->getIntData();
std::vector<int>::const_iterator actnum_begin = const int* gc = G.global_cell;
newParserDeck->getKeyword("ACTNUM")->getIntData().begin(); for (int cell = 0; cell < G.number_of_cells; ++cell) {
std::transform(e.begin(), e.end(), actnum_begin, std::back_inserter(eqlnum), const int deck_pos = (gc == NULL) ? cell : gc[cell];
[](int eql, int act){ return (act == 0) ? -1 : eql-1;}); eqlnum[cell] = e[deck_pos] - 1;
eqlnum.erase(std::remove_if(eqlnum.begin(), eqlnum.end(),
std::bind2nd(std::less<int>(), 0)),
eqlnum.end());
} }
else {
std::transform(e.begin(), e.end(), std::back_inserter(eqlnum),
std::bind2nd(std::minus<int>(), 1));
}
if (eqlnum.size() != G.number_of_cells)
OPM_THROW(std::runtime_error, "Keyword EQLNUM: inconsistent array size.");
} }
else { else {
// No explicit equilibration region. // No explicit equilibration region.