Riemann invariants for one-dimensional isentropic flow (source code)

= Riemann invariants for one-dimensional isentropic flow
{c}

For a barotropic pressure law with $c^2=dp/d\rho$, the invariants
$$
R_\pm=u\pm\int^\rho\frac{c(s)}s\,ds
$$
obey $[\partial_t+(u\pm c)\partial_x]R_\pm=0$.