From 34ad6b199266703f10013b510aa09350c7dbe41d Mon Sep 17 00:00:00 2001
From: Geoffroy Lesur <geoffroy.lesur@univ-grenoble-alpes.fr>
Date: Fri, 17 Oct 2025 22:25:52 +0200
Subject: [PATCH] fix "size" drag force with adiabatic equation of state (sound
 speed was not properly computed)

---
 src/fluid/drag.hpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

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