New Staggered with multicomponents needs a more sophisticated implementation in dumux!
To get the same boundary conditions as in the previous setup, we need to adapt the implementation in new staggered such that we can apply mixed boundary conditions in the mass problem. This was no issue when dealing with 1p1c systems so far. So this is a TODO for Dumux, I would maybe have to talk with Matthis...