diff --git a/src/hydro/PressureEnergy/hydro.h b/src/hydro/PressureEnergy/hydro.h index 68137d860f29ea529b94eb8d3026bea213e7703e..2671ac995409ea32981b535ca768c02f7fc0abce 100644 --- a/src/hydro/PressureEnergy/hydro.h +++ b/src/hydro/PressureEnergy/hydro.h @@ -153,7 +153,7 @@ hydro_get_comoving_soundspeed(const struct part *restrict p) { /* Compute the sound speed -- see theory section for justification */ /* IDEAL GAS ONLY -- P-U does not work with generic EoS. */ - const float square_rooted = sqrtf(p->rho / p->pressure_bar); + const float square_rooted = sqrtf(hydro_gamma * p->rho / p->pressure_bar); return p->u * hydro_gamma_minus_one * square_rooted; }