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

removing one unnecessary if condition

Browse Source
This commit is contained in:
Kai Bao
2017-04-06 14:53:44 +02:00
parent e627f35176
commit 8fa3c4b324
1 changed files with 17 additions and 19 deletions
Download Patch File Download Diff File Expand all files Collapse all files

36
opm/autodiff/StandardWellsDense_impl.hpp
View File

@@ -1835,26 +1835,24 @@ namespace Opm {
// update/setup guide rates for each well based on the well_potentials // update/setup guide rates for each well based on the well_potentials
well_collection_->setGuideRatesWithPotentials(wellsPointer(), phase_usage_, well_potentials); well_collection_->setGuideRatesWithPotentials(wellsPointer(), phase_usage_, well_potentials);
// handling the situation that wells does not have a valid control
// it happens the well specified with GRUP and restarting due to non-convergencing
// putting the well under group control for this situation
if (wellCollection()->groupControlActive()) {
for (int w = 0; w < nw; ++w) {
WellControls* wc = wells().ctrls[w];
const int ctrl_index = well_controls_get_current(wc);
WellNode& well_node = well_collection_->findWellNode(std::string(wells().name[w]));
const int group_control_index = well_node.groupControlIndex();
if (group_control_index >= 0 && ctrl_index < 0) {
well_controls_set_current(wc, group_control_index);
well_state.currentControls()[w] = group_control_index;
well_node.setIndividualControl(false);
}
}
}
} }
// handling the situation that wells do not have a valid control
// it happens the well specified with GRUP and restarting due to non-convergencing
// putting the well under group control for this situation
for (int w = 0; w < nw; ++w) {
WellControls* wc = wells().ctrls[w];
const int ctrl_index = well_controls_get_current(wc);
WellNode& well_node = well_collection_->findWellNode(std::string(wells().name[w]));
const int group_control_index = well_node.groupControlIndex();
if (group_control_index >= 0 && ctrl_index < 0) {
well_controls_set_current(wc, group_control_index);
well_state.currentControls()[w] = group_control_index;
well_node.setIndividualControl(false);
}
}
applyVREPGroupControl(well_state); applyVREPGroupControl(well_state);
if (!wellCollection()->groupControlApplied()) { if (!wellCollection()->groupControlApplied()) {