diff --git a/dumux/freeflow/stokes/stokeslocalresidual.hh b/dumux/freeflow/stokes/stokeslocalresidual.hh
index 5a983e145ec9081d398a96a5e6f4a4166abb8a99..fdc3c85761763259828dfaa5957e5814457b1fcd 100644
--- a/dumux/freeflow/stokes/stokeslocalresidual.hh
+++ b/dumux/freeflow/stokes/stokeslocalresidual.hh
@@ -227,11 +227,7 @@ protected:
             // TODO: dilatation term has to be accounted for in outflow, coupling, neumann
             //            velGradComp[velIdx] += 2./3*fluxVars.velocityDivAtIP;
 
-            if (massUpwindWeight_ == 1.0) // fully upwind
-                velGradComp *= up.viscosity();
-            else
-                velGradComp *= massUpwindWeight_ * up.viscosity() +
-                    (1.0 - massUpwindWeight_) * dn.viscosity();
+            velGradComp *= fluxVars.viscosityAtIP();
 
             flux[momentumXIdx + velIdx] -=
                 velGradComp*fluxVars.face().normal;