Exporting separate movie for each of the quantities of interest (i.e. density, velocity, temperature, mass flux)
Separate sub-step method grid is updated after each substep
Adding mass source term
Calculate and store all necessary time derivs at the beginning of each substep
Changing variables to log_rho