From ee06d4827bf6b7ddebc3c92aac7575c9f8ac970d Mon Sep 17 00:00:00 2001
From: Bert Vandenbroucke <bert.vandenbroucke@ugent.be>
Date: Fri, 19 Aug 2016 17:38:31 +0100
Subject: [PATCH] Added gravitational acceleration to primitive variable drift.
---
src/hydro/Gizmo/hydro.h | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/hydro/Gizmo/hydro.h b/src/hydro/Gizmo/hydro.h
index 69935d1a4e..7dd909351d 100644
--- a/src/hydro/Gizmo/hydro.h
+++ b/src/hydro/Gizmo/hydro.h
@@ -299,9 +299,9 @@ __attribute__((always_inline)) INLINE static void hydro_predict_extra(
} else {
p->primitives.rho *= expf(w);
}
- p->primitives.v[0] += p->a_hydro[0] * dt;
- p->primitives.v[1] += p->a_hydro[1] * dt;
- p->primitives.v[2] += p->a_hydro[2] * dt;
+ p->primitives.v[0] += (p->a_hydro[0] + p->gravity.old_a[0]) * dt;
+ p->primitives.v[1] += (p->a_hydro[1] + p->gravity.old_a[1]) * dt;
+ p->primitives.v[2] += (p->a_hydro[2] + p->gravity.old_a[2]) * dt;
float u = p->conserved.energy + p->du_dt * dt;
p->primitives.P =
hydro_gamma_minus_one * u * p->primitives.rho / p->conserved.mass;
--
GitLab