fipy
This example shows the failure of advecting a square pulse with a first order explicit upwind scheme.