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

Merge pull request #2096 from totto82/avoid_nan_density

Browse Source 
avoid nan in the density calculation
This commit is contained in:
Atgeirr Flø Rasmussen 2019-10-22 14:27:45 +02:00 committed by GitHub
commit 5d1e125ab9
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -2344,10 +2344,10 @@ namespace Opm
const unsigned oilpos = Indices::canonicalToActiveComponentIndex(FluidSystem::oilCompIdx); const unsigned oilpos = Indices::canonicalToActiveComponentIndex(FluidSystem::oilCompIdx);
double rs = 0.0; double rs = 0.0;
double rv = 0.0; double rv = 0.0;
if (!rsmax_perf.empty() && mix[oilpos] > 0.0) { if (!rsmax_perf.empty() && mix[oilpos] > 1e-12) {
rs = std::min(mix[gaspos]/mix[oilpos], rsmax_perf[perf]); rs = std::min(mix[gaspos]/mix[oilpos], rsmax_perf[perf]);
} }
if (!rvmax_perf.empty() && mix[gaspos] > 0.0) { if (!rvmax_perf.empty() && mix[gaspos] > 1e-12) {
rv = std::min(mix[oilpos]/mix[gaspos], rvmax_perf[perf]); rv = std::min(mix[oilpos]/mix[gaspos], rvmax_perf[perf]);
} }
if (rs != 0.0) { if (rs != 0.0) {