diff --git a/dumux/porousmediumflow/2p2c/primaryvariableswitch.hh b/dumux/porousmediumflow/2p2c/primaryvariableswitch.hh index 3b9282e91932f9e0be28ddab02857d5d81d00afc..703da6cb14871b652ab1f359cdeea665d17bb058 100644 --- a/dumux/porousmediumflow/2p2c/primaryvariableswitch.hh +++ b/dumux/porousmediumflow/2p2c/primaryvariableswitch.hh @@ -107,9 +107,9 @@ protected: << xww + xwn << std::endl; newPhasePresence = bothPhases; if (formulation == pnsw) - priVars[switchIdx] = 0.0; + priVars[switchIdx] = 0.0001; else if (formulation == pwsn) - priVars[switchIdx] = 1.0; + priVars[switchIdx] = 0.9999; } } else if (phasePresence == wPhaseOnly) @@ -135,9 +135,9 @@ protected: << xnw + xnn << std::endl; newPhasePresence = bothPhases; if (formulation == pnsw) - priVars[switchIdx] = 0.999; + priVars[switchIdx] = 0.9999; else if (formulation == pwsn) - priVars[switchIdx] = 0.001; + priVars[switchIdx] = 0.0001; } } else if (phasePresence == bothPhases)