diff --git a/dumux/freeflow/stokes/localresidual.hh b/dumux/freeflow/stokes/localresidual.hh index 93b9e4d833c07c703c5ac799a184a87974a2942d..7951945fc9ad49b28c10fef6210f5e3faee809a4 100644 --- a/dumux/freeflow/stokes/localresidual.hh +++ b/dumux/freeflow/stokes/localresidual.hh @@ -519,6 +519,7 @@ protected: PrimaryVariables values(0.0); asImp_()->computeFlux(values, boundaryFaceIdx, /*onBoundary=*/true); + values *= this->curVolVars_(scvIdx).extrusionFactor(); Valgrind::CheckDefined(values); for (int eqIdx = 0; eqIdx < numEq; ++eqIdx) diff --git a/dumux/implicit/box/localresidual.hh b/dumux/implicit/box/localresidual.hh index 4b69761d836ed384918e3209c2c968a4abfbfdc5..1eacceb3b1711de1024d04f5ab68b0717caa6627 100644 --- a/dumux/implicit/box/localresidual.hh +++ b/dumux/implicit/box/localresidual.hh @@ -204,6 +204,7 @@ protected: //calculate outflow fluxes PrimaryVariables values(0.0); this->asImp_().computeFlux(values, boundaryFaceIdx, true); + values *= this->curVolVars_(scvIdx).extrusionFactor(); Valgrind::CheckDefined(values); for (int equationIdx = 0; equationIdx < numEq; ++equationIdx)