EulerEq: More follow ups to !155, !156

src/EulerEq/equation.jl

@@ -6,12 +6,13 @@
 @with_signature function flux(equation::EulerEquation)
   @accepts rho, q, E
   @unpack eos = equation
+  u = q/rho
   eps = E / rho - q^2 / (2 * rho^2)
   p = eos(Pressure, Density, rho, InternalEnergy, eps)
   flx_rho = q
   flx_q   = q^2 / rho + p
   flx_E   = q / rho * (E + p)
-  @returns flx_rho, flx_q, flx_E, p
+  @returns flx_rho, flx_q, flx_E, p, u
 end
 
 @with_signature function fv_bdry_flux(equation::EulerEquation)