This models an unstable shear flow. The model is initialized with a strong shear, then is given a small perturbation, resulting in the flow going unstable and the creation of Kelvin-Helmholtz billows.

The code outputs the fluid velocity, vorticity and particle trajectories to a NetCDF file. A Matlab script is included to recreate the movie shown above.

The code is available for download from GitHub