Merged requested to merge feature/improve-upwind-scheme into next
Now we just tell the flux variables which quantities should be receive upwind treatment. The flux variables will decide how to upwind. This way upwind gets more flexible. This is necessary to treat e.g. bifurcations differently where we don't have only two participating volvars but still need to know the physical quantities to upwind from the local residual.