correcting a bug in computeWellFlux in BlackoilMultiSegmentModel

it should be segment rates instead of well rates.
This commit is contained in:
Kai Bao 2015-09-24 16:46:08 +02:00
parent 02b187e14c
commit 5aa258fdd6

View File

@ -681,7 +681,7 @@ namespace Opm {
// const ADB& q_ps = well->wellOps().p2s * cq_ps[phase];
const ADB& q_ps = well->wellOps().p2s_gather * cq_ps[phase];
const int n_total_segments = state.segp.size();
const ADB& q_s = subset(state.qs, Span(nseg, 1, phase * n_total_segments + start_segment));
const ADB& q_s = subset(state.segqs, Span(nseg, 1, phase * n_total_segments + start_segment));
Selector<double> injectingPhase_selector(q_s.value(), Selector<double>::GreaterZero);
const int pos = pu.phase_pos[phase];