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

Merge pull request #6010 from vkip/allow_empty_non_group_leaf_node

Browse Source 
Allow empty non-group leaf nodes
This commit is contained in:
Kai Bao authored and GitHub committed 2025-02-19 11:13:11 +01:00
2 parents f40d88a284 0a36678709
commit 86cc12bece
1 file changed
+6 -5
Download Patch File Download Diff File Expand all files Collapse all files
opm/simulators/wells/WellGroupHelpers.cpp
+6 -5
View File
@@ -755,12 +755,13 @@ updateNetworkLeafNodeProductionRates(const Schedule& schedule,
if (network.active()) {
const int np = wellState.numPhases();
for (const auto& group_name : network.leaf_nodes()) {
assert(schedule[reportStepIdx].groups.has(group_name));
const auto& group = schedule[reportStepIdx].groups.get(group_name);
std::vector<Scalar> network_rates(np, 0.0);
if (group.numWells() > 0) {
for (int phase = 0; phase < np; ++phase) {
network_rates[phase] = sumWellPhaseRates(false, group, schedule, wellState, reportStepIdx, phase, /*isInjector*/ false, /*network*/ true);
if (schedule[reportStepIdx].groups.has(group_name)) { // Allow empty leaf nodes that are not groups
const auto& group = schedule[reportStepIdx].groups.get(group_name);
if (group.numWells() > 0) {
for (int phase = 0; phase < np; ++phase) {
network_rates[phase] = sumWellPhaseRates(false, group, schedule, wellState, reportStepIdx, phase, /*isInjector*/ false, /*network*/ true);
}
}
}
group_state.update_network_leaf_node_production_rates(group_name, network_rates);