OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-03 21:16:54 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #1895 from joakim-hove/apply-actions

Browse Source 
Add applyActions() to EclProblem
This commit is contained in:
Joakim Hove 2019-08-29 08:58:06 +02:00 committed by GitHub
commit 8992476651
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 30 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
ebos/eclproblem.hh
View File

@ -87,6 +87,7 @@
#include <opm/parser/eclipse/EclipseState/EclipseState.hpp> #include <opm/parser/eclipse/EclipseState/EclipseState.hpp>
#include <opm/parser/eclipse/EclipseState/Tables/Eqldims.hpp> #include <opm/parser/eclipse/EclipseState/Tables/Eqldims.hpp>
#include <opm/parser/eclipse/EclipseState/Schedule/Schedule.hpp> #include <opm/parser/eclipse/EclipseState/Schedule/Schedule.hpp>
#include <opm/parser/eclipse/EclipseState/Schedule/Action/ActionContext.hpp>
#include <opm/parser/eclipse/EclipseState/SummaryConfig/SummaryConfig.hpp> #include <opm/parser/eclipse/EclipseState/SummaryConfig/SummaryConfig.hpp>
#include <opm/parser/eclipse/EclipseState/Tables/RockwnodTable.hpp> #include <opm/parser/eclipse/EclipseState/Tables/RockwnodTable.hpp>
#include <opm/parser/eclipse/EclipseState/Tables/OverburdTable.hpp> #include <opm/parser/eclipse/EclipseState/Tables/OverburdTable.hpp>
@ -966,10 +967,13 @@ public:
void endEpisode() void endEpisode()
{ {
auto& simulator = this->simulator(); auto& simulator = this->simulator();
const auto& schedule = simulator.vanguard().schedule(); auto& schedule = simulator.vanguard().schedule();
const auto& timeMap = schedule.getTimeMap(); const auto& timeMap = schedule.getTimeMap();
int episodeIdx = simulator.episodeIndex(); int episodeIdx = simulator.episodeIndex();
this->applyActions(episodeIdx + 1,
schedule,
simulator.vanguard().summaryState());
// check if we're finished ... // check if we're finished ...
if (episodeIdx + 1 >= static_cast<int>(timeMap.numTimesteps())) { if (episodeIdx + 1 >= static_cast<int>(timeMap.numTimesteps())) {
@ -1013,6 +1017,31 @@ public:
eclWriter_->writeOutput(isSubStep); eclWriter_->writeOutput(isSubStep);
} }
void applyActions(int reportStep,
Opm::Schedule& schedule,
const Opm::SummaryState& summaryState) {
const auto& actions = schedule.actions();
if (actions.empty())
return;
Opm::Action::Context context( summaryState );
auto simTime = schedule.simTime(reportStep);
for (const auto& action : actions.pending(simTime)) {
auto actionResult = action->eval(simTime, context);
if (actionResult) {
std::string msg = "The action: " + action->name() + " evaluated to true at report step: " + std::to_string(reportStep);
Opm::OpmLog::info(msg);
schedule.applyAction(reportStep, *action, actionResult);
} else {
std::string msg = "The action: " + action->name() + " evaluated to false sat report step: " + std::to_string(reportStep);
Opm::OpmLog::info(msg);
}
}
}
/*! /*!
* \copydoc FvBaseMultiPhaseProblem::intrinsicPermeability * \copydoc FvBaseMultiPhaseProblem::intrinsicPermeability
*/ */