...
 
Commits (3)
...@@ -175,7 +175,8 @@ public: ...@@ -175,7 +175,8 @@ public:
// Upwinding conductivity factor // Upwinding conductivity factor
RangeU cond_upwind_factor(0.0); 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 // inward flux
cond_upwind_factor = cond_o / cond_sat_o; cond_upwind_factor = cond_o / cond_sat_o;
} }
...@@ -190,7 +191,7 @@ public: ...@@ -190,7 +191,7 @@ public:
else // RichardsFVUpwinding::fullUpwind else // RichardsFVUpwinding::fullUpwind
{ {
// No average: Use upwind conductivity // No average: Use upwind conductivity
if (dudn < 0.0) { if (dudn > 0.0) {
// inward flux // inward flux
cond_f = cond_o; cond_f = cond_o;
} }
......