Avoid segfault for two-phase case.

This commit is contained in:
Atgeirr Flø Rasmussen 2020-05-29 10:20:58 +02:00
parent 61eb74b1d6
commit b31f33f567

View File

@ -229,6 +229,9 @@ public:
}
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
if (!FluidSystem::phaseIsActive(phaseIdx)) {
continue;
}
if (pressureOutput_())
pressure_[phaseIdx][I] = Opm::getValue(fs.pressure(phaseIdx));
if (densityOutput_())