diff --git a/dumux/boxmodels/common/boxlocalresidual.hh b/dumux/boxmodels/common/boxlocalresidual.hh
index d64c6cb9de4c21b480f29b74b442a1f38a266ecb..fb48b16e5a8258c73b2c63ab990d653f0ee5eaf8 100644
--- a/dumux/boxmodels/common/boxlocalresidual.hh
+++ b/dumux/boxmodels/common/boxlocalresidual.hh
@@ -442,7 +442,13 @@ protected:
                 fvElemGeom_().boundaryFace[boundaryFaceIdx].area
                 * curVolVars_(scvIdx).extrusionFactor();
             Valgrind::CheckDefined(values);
-            residual_[scvIdx] += values;
+
+            // set the neumann conditions
+            for (int eqIdx = 0; eqIdx < numEq; ++eqIdx) {
+                if (!bcTypes.isNeumann(eqIdx))
+                    continue;
+                residual_[scvIdx] += values;
+            }
         }
     }