Open
Description
Describe the bug
The HLLD solver sets the flux for passive scalars (including mass scalars) to zero:
Line 329 in edc1dbc
(Within the Riemann solver, the passive scalars should be treated the same as the transverse momentum components.)
Additional context
Problems with passive scalars or mass scalars will not work correctly until this is added.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment