First version of advection with 6th order interpolation. What to do with
scalar advection near boundary? Using ghost cell data provides very
large gradients and problems with interpolation. For now the vertical
ghost cells are not used; interpolation stencil decreased to prevent
using them