diff --git a/opm/models/blackoil/blackoilintensivequantities.hh b/opm/models/blackoil/blackoilintensivequantities.hh index e1223b4a6..bd59dbec3 100644 --- a/opm/models/blackoil/blackoilintensivequantities.hh +++ b/opm/models/blackoil/blackoilintensivequantities.hh @@ -348,7 +348,11 @@ public: if(FluidSystem::phaseIsActive(oilPhaseIdx)){ x = rockCompressibility*(fluidState_.pressure(oilPhaseIdx) - rockRefPressure); }else{ - x = rockCompressibility*(fluidState_.pressure(waterPhaseIdx) - rockRefPressure); + if(FluidSystem::phaseIsActive(waterPhaseIdx)){ + x = rockCompressibility*(fluidState_.pressure(waterPhaseIdx) - rockRefPressure); + }else{ + x = rockCompressibility*(fluidState_.pressure(gasPhaseIdx) - rockRefPressure); + } } porosity_ *= 1.0 + x + 0.5*x*x; }