It seems that the well-balanced scheme doesn't work anymore. If I run: ``` ./pyro.py compressible_rk hse inputs.hse vis.dovis=1 compressible.well_balanced=1 ``` it goes 16 steps, and has all NaNs in the velocity in the output