Bugfix: do not compute MDU terms for noflow faces.

This commit is contained in:
Atgeirr Flø Rasmussen
2013-04-17 10:38:08 +02:00
parent 40fe2abf04
commit f8db9cad66
+1 -1
View File
@@ -219,7 +219,7 @@ namespace Opm
// Add flux to upwind_term or downwind_term_[face|cell_factor].
if (flux < 0.0) {
upwind_term += flux*face_tof_[f];
} else {
} else if (flux > 0.0) {
double fterm, cterm_factor;
multidimUpwindTerms(f, cell, fterm, cterm_factor);
downwind_term_face += fterm*flux;