Commit d33b4473 authored by Lukas Riedel's avatar Lukas Riedel

Merge branch '159-follow-up-from-resolve-add-upwinding-option-to-richards-fv-solver' into 'master'

Resolve "Follow-up from "Resolve "Add upwinding option to Richards FV solver"""

Closes #159

See merge request !162
parents 373885e6 a49788e6
......@@ -173,7 +173,8 @@ public:
// Upwinding conductivity factor
RangeU cond_upwind_factor(0.0);
if (dudn < 0.0) {
// NOTE: Sign seems inverted due to definition of 'water_flux_n' below.
if (dudn > 0.0) {
// inward flux
cond_upwind_factor = cond_o / cond_sat_o;
}
......@@ -188,7 +189,7 @@ public:
else // RichardsFVUpwinding::fullUpwind
{
// No average: Use upwind conductivity
if (dudn < 0.0) {
if (dudn > 0.0) {
// inward flux
cond_f = cond_o;
}
......
Markdown is supported
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