diff --git a/applications/ebos/ecloutputblackoilmodule.hh b/applications/ebos/ecloutputblackoilmodule.hh index 36a812e48..876f54121 100644 --- a/applications/ebos/ecloutputblackoilmodule.hh +++ b/applications/ebos/ecloutputblackoilmodule.hh @@ -167,11 +167,7 @@ public: const auto &fs = elemCtx.intensiveQuantities(dofIdx, /*timeIdx=*/0).fluidState(); typedef typename std::remove_const::type>::type FluidState; unsigned globalDofIdx = elemCtx.globalSpaceIndex(dofIdx, /*timeIdx=*/0); - unsigned regionIdx = elemCtx.primaryVars(dofIdx, /*timeIdx=*/0).pvtRegionIndex(); - Scalar po = Toolbox::value(fs.pressure(oilPhaseIdx)); - Scalar To = Toolbox::value(fs.temperature(oilPhaseIdx)); - Scalar XoG = Toolbox::value(fs.massFraction(oilPhaseIdx, gasCompIdx)); - Scalar XgO = Toolbox::value(fs.massFraction(gasPhaseIdx, oilCompIdx)); + unsigned pvtRegionIdx = elemCtx.primaryVars(dofIdx, /*timeIdx=*/0).pvtRegionIndex(); if (saturationsOutput_()) { for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) { @@ -187,22 +183,22 @@ public: } if (gasDissolutionFactorOutput_()) { gasDissolutionFactor_[globalDofIdx] = - FluidSystem::template saturatedDissolutionFactor(fs, gasPhaseIdx, regionIdx); + FluidSystem::template saturatedDissolutionFactor(fs, gasPhaseIdx, pvtRegionIdx); Valgrind::CheckDefined(gasDissolutionFactor_[globalDofIdx]); } if (gasFormationVolumeFactorOutput_()) { gasFormationVolumeFactor_[globalDofIdx] = - FluidSystem::template formationVolumeFactor(fs, gasPhaseIdx, regionIdx); + FluidSystem::template formationVolumeFactor(fs, gasPhaseIdx, pvtRegionIdx); Valgrind::CheckDefined(gasFormationVolumeFactor_[globalDofIdx]); } if (saturatedOilFormationVolumeFactorOutput_()) { saturatedOilFormationVolumeFactor_[globalDofIdx] = - FluidSystem::template saturatedFormationVolumeFactor(fs, oilPhaseIdx, regionIdx); + FluidSystem::template saturatedFormationVolumeFactor(fs, oilPhaseIdx, pvtRegionIdx); Valgrind::CheckDefined(saturatedOilFormationVolumeFactor_[globalDofIdx]); } if (oilSaturationPressureOutput_()) { oilSaturationPressure_[globalDofIdx] = - FluidSystem::template saturationPressure(fs, oilPhaseIdx, regionIdx); + FluidSystem::template saturationPressure(fs, oilPhaseIdx, pvtRegionIdx); Valgrind::CheckDefined(oilSaturationPressure_[globalDofIdx]); } }