OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Include EQUIL-based inititalisation for 3-phase.

Browse Source
This commit is contained in:
osae
2014-06-26 15:11:34 +02:00
parent caa49bf3ce
commit 24506cc637
1 changed files with 6 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
examples/sim_fibo_ad.cpp
View File

@@ -26,6 +26,7 @@
#include <opm/core/wells/WellsManager.hpp> #include <opm/core/wells/WellsManager.hpp>
#include <opm/core/utility/ErrorMacros.hpp> #include <opm/core/utility/ErrorMacros.hpp>
#include <opm/core/simulator/initState.hpp> #include <opm/core/simulator/initState.hpp>
#include <opm/core/simulator/initStateEquil.hpp>
#include <opm/core/simulator/SimulatorReport.hpp> #include <opm/core/simulator/SimulatorReport.hpp>
#include <opm/core/simulator/SimulatorTimer.hpp> #include <opm/core/simulator/SimulatorTimer.hpp>
#include <opm/core/utility/miscUtilities.hpp> #include <opm/core/utility/miscUtilities.hpp>
@@ -135,6 +136,11 @@ try
/ state.surfacevol()[c*np + pu.phase_pos[Oil]]; / state.surfacevol()[c*np + pu.phase_pos[Oil]];
} }
} }
} else if (deck->hasKeyword("EQUIL") && props->numPhases() == 3) {
state.init(*grid->c_grid(), props->numPhases());
const double grav = param.getDefault("gravity", unit::gravity);
initStateEquil(*grid->c_grid(), *props, deck, eclipseState, grav, state);
state.faceflux().resize(grid->c_grid()->number_of_faces, 0.0);
} else { } else {
initBlackoilStateFromDeck(*grid->c_grid(), *props, deck, gravity[2], state); initBlackoilStateFromDeck(*grid->c_grid(), *props, deck, gravity[2], state);
} }