Implement non-conforming output for box with interface solver
As soon as we have an interface solver (!907 (merged)), we should think about how to visualize correctly the results. The original idea was element-averaged saturations, but an alternative could be nodal, but nonconforming vtk output.