OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-07-07 04:53:03 -05:00
Code Issues Packages Projects Releases Wiki Activity

Rewrite to properly loop over components

Browse Source
This commit is contained in:
Stein Krogstad 2024-06-25 11:35:25 +02:00
parent 1efcd2c936
commit fd199d2558
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
opm/simulators/wells/StandardWell_impl.hpp
View File

@ -1655,9 +1655,10 @@ namespace Opm
// fetch potentials (sign is updated on the outside).
well_potentials.clear();
well_potentials.resize(np, 0.0);
for (int phase = 0; phase < np; ++phase) {
const EvalWell rate = well_copy.primary_variables_.getQs(phase);
well_potentials[this->modelCompIdxToFlowCompIdx(phase)] = rate.value();
for (int comp_idx = 0; comp_idx < this->num_components_; ++comp_idx) {
if (has_solvent && comp_idx == Indices::contiSolventEqIdx) continue; // we do not store the solvent in the well_potentials
const EvalWell rate = well_copy.primary_variables_.getQs(comp_idx);
well_potentials[this->modelCompIdxToFlowCompIdx(comp_idx)] = rate.value();
}
// the solvent contribution is added to the gas potentials