Broadcast for `var_on_centers` and `var_on_vertices`

[bmxam]

We should allow

values_u, values_v = var_on_centers(mesh, (u, v))

because up to now we have to call var_on_centers twice in a row for this (so we loop twice on the mesh), which is a bit stupid.