Do not try to compute initial solution for inactive multi-segment wells

This commit is contained in:
Vegard Kippe
2024-11-19 15:57:49 +01:00
parent 5b002b2a77
commit c42f02e802

View File

@@ -724,6 +724,12 @@ void WellState<Scalar>::initWellStateMSWell(const std::vector<Well>& wells_ecl,
n_activeperf++; n_activeperf++;
} }
} }
// If perf_data does not contain all active perfs, the well is distributed or inactive
// For inactive wells there is no need for calculating an initial solution
// \todo{ Update the procedure below to work for actually distributed wells. }
if (static_cast<int>(ws.perf_data.size()) != n_activeperf)
continue;
std::vector<std::vector<int>> segment_inlets(well_nseg); std::vector<std::vector<int>> segment_inlets(well_nseg);
for (int seg = 0; seg < well_nseg; ++seg) { for (int seg = 0; seg < well_nseg; ++seg) {