[upwind] Define the upwing quantities instead of the rule in local residual

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.