The interpolation method in a semi-Lagrangian scheme is decisive to its performance. Given the number of grid points one is considering to use for the interpolation, it does not necessarily follow that maximum formal accuracy should give the best results. For the advection equation, the driving force of this method is the method of the characteristics, which accounts for the flow of information in the model equation. This leads naturally to an interpolation problem since the foot point is not in general located on a grid point. We use another interpolation scheme that will allow achieving the high order for the box initial condition.
 Morton, K.W. (1980) Stability of Finite Difference Approximations to a Diffusion-Convection Equation. International Journal for Numerical Methods in Engineering, 15, 677-683.
 Hundsdorfer, W. and Koren, B. (1995) A Positive Finite-Difference Advection Scheme Applied on Locally Refined Grids. Journal of Computational Physics, 117, 35-36. http://dx.doi.org/10.1006/jcph.1995.1042
 Canuto, C. and Hussaini, M. (1988) Spectral Methods in Fluids Dynamics. Springer Series in Computational Physics, Springer-Verlag, Berlin. http://dx.doi.org/10.1007/978-3-642-84108-8
 Dehghan, M. (2005) On the Numerical Solution of the One-Dimensional Convection-Diffusion Equation. Mathematical Problems in Engineering, 1, 61-74. http://dx.doi.org/10.1155/MPE.2005.61