Reduced number of loops in assemble_explicit()
Adding gravitational force term. You can use "set gravity = 1.0" in input file to enable the force.
Gravity source term in energy equation was wrong. W[dim-1] = rho * velocity but we were adding another factor of density in the source term.