diff --git a/dumux/multidomain/boundary/stokesdarcy/couplingdata.hh b/dumux/multidomain/boundary/stokesdarcy/couplingdata.hh index 259d7da18e1e0406372e99e487246246b2ae190d..ca88f85ccfd1a99abcb8271ec22ff33f8d2ee963 100644 --- a/dumux/multidomain/boundary/stokesdarcy/couplingdata.hh +++ b/dumux/multidomain/boundary/stokesdarcy/couplingdata.hh @@ -299,7 +299,7 @@ public: const Scalar mu = stokesContext.volVars.viscosity(darcyPhaseIdx); const Scalar rho = stokesContext.volVars.density(darcyPhaseIdx); const Scalar distance = (stokesContext.element.geometry().center() - scvf.center()).two_norm(); - const Scalar g = couplingManager_.problem(darcyIdx).gravity()[scvf.directionIndex()]; + const Scalar g = -scvf.directionSign() * couplingManager_.problem(darcyIdx).gravity()[scvf.directionIndex()]; const Scalar interfacePressure = ((scvf.directionSign() * velocity * (mu/darcyPermeability(scvf))) + rho * g) * distance + darcyPressure; momentumFlux = interfacePressure; }