changed d[energy]/d[rho] flux entry to be based on total energy instead of total enthalpy