Extension to fluid solver: Output interfacial advection speed

[ekeever1]

Accurate advection of passive fields which is self-consistent with Godunov methods will require that the Riemann solver output the interfacial flux velocities to be used in calculating the passive scalar's advection flux.