Skip to content
Snippets Groups Projects
Commit c224d57b authored by Bernd Flemisch's avatar Bernd Flemisch
Browse files

Merge branch 'fix/priVar-iniSat-2p2c' into 'master' 

[2p2c] unify 2p2c and 2pnc init saturation guess

Closes #407

See merge request !881
parents ea9649be cc28be8d
No related branches found
No related tags found
1 merge request!881[2p2c] unify 2p2c and 2pnc init saturation guess
Hide whitespace changes
Inline Side-by-side
dumux/porousmediumflow/2p2c/primaryvariableswitch.hh
+ 4
− 4
View file @ c224d57b
...@@ -107,9 +107,9 @@ protected: ...@@ -107,9 +107,9 @@ protected:
<< xww + xwn << std::endl; << xww + xwn << std::endl;
newPhasePresence = bothPhases; newPhasePresence = bothPhases;
if (formulation == pnsw) if (formulation == pnsw)
priVars[switchIdx] = 0.0; priVars[switchIdx] = 0.0001;
else if (formulation == pwsn) else if (formulation == pwsn)
priVars[switchIdx] = 1.0; priVars[switchIdx] = 0.9999;
} }
} }
else if (phasePresence == wPhaseOnly) else if (phasePresence == wPhaseOnly)
...@@ -135,9 +135,9 @@ protected: ...@@ -135,9 +135,9 @@ protected:
<< xnw + xnn << std::endl; << xnw + xnn << std::endl;
newPhasePresence = bothPhases; newPhasePresence = bothPhases;
if (formulation == pnsw) if (formulation == pnsw)
priVars[switchIdx] = 0.999; priVars[switchIdx] = 0.9999;
else if (formulation == pwsn) else if (formulation == pwsn)
priVars[switchIdx] = 0.001; priVars[switchIdx] = 0.0001;
} }
} }
else if (phasePresence == bothPhases) else if (phasePresence == bothPhases)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment