diff --git a/dumux/implicit/common/implicitvelocityoutput.hh b/dumux/implicit/common/implicitvelocityoutput.hh index 5466be50cd78dc4f445da702b8b2926585f23941..5de58cd233a434cd50473eec2cd8174f9f0cd2e0 100644 --- a/dumux/implicit/common/implicitvelocityoutput.hh +++ b/dumux/implicit/common/implicitvelocityoutput.hh @@ -218,7 +218,7 @@ public: elemVolVars); Scalar flux = fluxVars.volumeFlux(phaseIdx); - scvVelocities[fIdx] = flux; + scvVelocities[fIdx] += flux; fIdxInner++; } @@ -239,7 +239,7 @@ public: for (int i = 0; i < dim; i++) refVelocity[i] = 0.5 * (scvVelocities[2*i + 1] - scvVelocities[2*i]); - Dune::FieldVector<Scalar, dim> scvVelocity(0); + Dune::FieldVector<Scalar, dimWorld> scvVelocity(0); jacobianT2.mtv(refVelocity, scvVelocity); scvVelocity /= element.geometry().integrationElement(localPos);