Currently, they use `defs.rho0`. This should be an optional parameter (defaulting to `defs.rho0`), just like the speed of sound.