Merge pull request #3341 from joakim-hove/explicitly-clear-well-potentials

Explicitly zero well potentials at the start of report step
This commit is contained in:
Bård Skaflestad 2021-06-02 18:14:41 +02:00 committed by GitHub
commit a03d7589b2
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23

View File

@ -337,17 +337,13 @@ void WellState::init(const std::vector<double>& cellPressures,
}
}
perfRateSolvent_.clear();
perfRateSolvent_.resize(nperf, 0.0);
productivity_index_.resize(nw * this->numPhases(), 0.0);
conn_productivity_index_.resize(nperf * this->numPhases(), 0.0);
well_potentials_.resize(nw * this->numPhases(), 0.0);
perfRateSolvent_.assign(nperf, 0.0);
productivity_index_.assign(nw * this->numPhases(), 0.0);
conn_productivity_index_.assign(nperf * this->numPhases(), 0.0);
well_potentials_.assign(nw * this->numPhases(), 0.0);
perfRatePolymer_.clear();
perfRatePolymer_.resize(nperf, 0.0);
perfRateBrine_.clear();
perfRateBrine_.resize(nperf, 0.0);
perfRatePolymer_.assign(nperf, 0.0);
perfRateBrine_.assign(nperf, 0.0);
for (int w = 0; w < nw; ++w) {
switch (wells_ecl[w].getStatus()) {