diff --git a/dumux/porousmediumflow/3pwateroil/model.hh b/dumux/porousmediumflow/3pwateroil/model.hh index 557061de497d4a510f044b1f6e187c66aed03d75..b3bdb846cb1742f38c2608e7164eb91de348e844 100644 --- a/dumux/porousmediumflow/3pwateroil/model.hh +++ b/dumux/porousmediumflow/3pwateroil/model.hh @@ -118,35 +118,47 @@ struct ThreePWaterOilModelTraits switch (state) { case Indices::threePhases: + { const std::vector<std::string> s1 = {"p_g", "S_w", "S_n"}; return s1[pvIdx]; + } case Indices::wPhaseOnly: + { const std::vector<std::string> s2 = {"p_w", "T", "x^" + FluidSystem::componentName(FluidSystem::nCompIdx) + "_" + FluidSystem::phaseName(FluidSystem::wPhaseIdx)}; return s2[pvIdx]; + } case Indices::gnPhaseOnly: + { const std::vector<std::string> s3 = {"p_g", "S_n", "x^" + FluidSystem::componentName(FluidSystem::wCompIdx) + "_" + FluidSystem::phaseName(FluidSystem::nPhaseIdx)}; return s3[pvIdx]; + } case Indices::wnPhaseOnly: + { const std::vector<std::string> s4 = {"p_w", "T", "S_n"}; return s4[pvIdx]; + } case Indices::gPhaseOnly: + { const std::vector<std::string> s5 = {"p_g", "T", "x^" + FluidSystem::componentName(FluidSystem::nCompIdx) + "_" + FluidSystem::phaseName(FluidSystem::gPhaseIdx)}; return s5[pvIdx]; + } case Indices::wgPhaseOnly: + { const std::vector<std::string> s6 = {"p_g", "S_w", "x^" + FluidSystem::componentName(FluidSystem::nCompIdx) + "_" + FluidSystem::phaseName(FluidSystem::gPhaseIdx)}; return s6[pvIdx]; + } } } };