Abstract
Errors in the numerical analog to a simple convection model can lead to an instability, even though existing methods of analysis predict that the scheme is stable.
The solutions obtained from a number of finite difference schemes and different resolutions are compared. This shows that increasing the resolution does not significantly improve the results, also that increasing the order of the finite difference scheme can lead to greater errors. The mechanism of the instability is discussed.