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

making the single phase RESV injection work.

Browse Source
This commit is contained in:
Kai Bao
2017-11-24 11:13:02 +01:00
parent 796f6cd13c
commit 59722e552c
1 changed files with 11 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
opm/autodiff/BlackoilWellModel_impl.hpp
View File

@@ -1151,7 +1151,6 @@ namespace Opm {
std::vector<double> distr (np);
std::vector<double> hrates(np);
std::vector<double> prates(np);
for (std::vector<int>::const_iterator
rp = resv_wells.begin(), e = resv_wells.end();
@@ -1169,21 +1168,20 @@ namespace Opm {
if (0 <= rctrl) {
const std::vector<double>::size_type off = (*rp) * np;
if (is_producer) {
// Convert to positive rates to avoid issues
// in coefficient calculations.
std::transform(well_state_.wellRates().begin() + (off + 0*np),
well_state_.wellRates().begin() + (off + 1*np),
prates.begin(), std::negate<double>());
} else {
std::copy(well_state_.wellRates().begin() + (off + 0*np),
well_state_.wellRates().begin() + (off + 1*np),
prates.begin());
}
const int fipreg = 0; // Hack. Ignore FIP regions.
rateConverter_->calcCoeff(fipreg, pvtreg, distr);
if (!is_producer) { // injectors
well_controls_assert_number_of_phases(ctrl, np);
// original distr contains 0 and 1 to indicate phases under control
const double* old_distr = well_controls_get_current_distr(ctrl);
for (int p = 0; p < np; ++p) {
distr[p] *= old_distr[p];
}
}
well_controls_iset_distr(ctrl, rctrl, & distr[0]);
}
}