updating pressure and phase rates for top segments in updateWellControls

This commit is contained in:
Kai Bao 2015-10-07 17:19:52 +02:00
parent 70d6d253fe
commit 2b4922a986

View File

@ -994,6 +994,7 @@ namespace Opm {
switch (well_controls_iget_type(wc, current)) {
case BHP:
xw.bhp()[w] = target;
xw.segPress()[xw.topSegmentLoc()[w]] = target;
break;
case THP: {
@ -1051,7 +1052,8 @@ namespace Opm {
case SURFACE_RATE:
for (int phase = 0; phase < np; ++phase) {
if (distr[phase] > 0.0) {
xw.wellRates()[np*w + phase] = target * distr[phase];
xw.wellRates()[np * w + phase] = target * distr[phase];
xw.segPhaseRates()[np * xw.topSegmentLoc()[w] + phase] = target * distr[phase];
}
}
break;