Adding some extra const-ness where it was appropriate to grid
Advection core now works, with explicit euler and runge-kutta 2 working as examples.