diff --git a/src/fluid/drag.hpp b/src/fluid/drag.hpp index 29eca309..63b5e1bd 100644 --- a/src/fluid/drag.hpp +++ b/src/fluid/drag.hpp @@ -37,8 +37,8 @@ class GammaDrag { // Assume a fixed size, hence for both Epstein or Stokes, gamma~1/rho_g/cs // Get the sound speed #if HAVE_ENERGY == 1 - cs = std::sqrt(eos.GetGamma(VcGas(PRS,k,j,i),VcGas(RHO,k,j,i) - *VcGas(PRS,k,j,i)/VcGas(RHO,k,j,i))); + cs = std::sqrt( eos.GetGamma(VcGas(PRS,k,j,i),VcGas(RHO,k,j,i)) + *VcGas(PRS,k,j,i)/VcGas(RHO,k,j,i)); #else cs = eos.GetWaveSpeed(k,j,i); #endif