OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-27 09:40:59 -06:00
Code Issues Packages Projects Releases Wiki Activity

Initialize pressure in top segment for non MSW well to zero

Browse Source
This commit is contained in:
Joakim Hove 2021-06-01 07:29:03 +02:00
parent df7d96494b
commit fc185afc02
1 changed files with 1 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/simulators/wells/WellState.cpp
View File

@ -531,7 +531,6 @@ void WellState::init(const std::vector<double>& cellPressures,
for (int w = 0; w < nw; ++w) {
top_segment_index_[w] = w;
seg_number_[w] = 1; // Top segment is segment #1
this->seg_press_[w] = this->bhp(w);
}
//seg_rates_ = wellRates();
seg_rates_.assign(nw*np, 0);
@ -841,7 +840,7 @@ void WellState::initWellStateMSWell(const std::vector<Well>& wells_ecl,
if ( !well_ecl.isMultiSegment() ) { // not multi-segment well
nseg_ += 1;
seg_number_.push_back(1); // Assign single segment (top) as number 1.
seg_press_.push_back(bhp(w));
seg_press_.push_back(0);
for (int p = 0; p < np; ++p) {
seg_rates_.push_back(rates[p]);
}