Iterating over grid, compute boundary fluxes

Add some slides/exercise on how to do that