Here is the f90 code of the model described in Pivato et al, 2014. This version provides the animation described at the end of the article ( .plt file, use tecplot to post process)
https://drive.google.com/file/d/0ByVlxqWw8eVkaEYxNHlTbHlDMUk/view?usp=sharing
