This code is the modified version of Jos Stam's code.
http://www.dgp.toronto.edu/people/stam/reality/Research/zip/CDROM_GDC03.zip

In this code, "Combined Lagrangian-Eulerian Approach for Accurate Advection" is 
used for advection of density field. 