From c4bf296bae3f66a52880089f9135da42e7994399 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <schaller@strw.leidenuniv.nl>
Date: Wed, 18 Dec 2019 13:35:44 +0100
Subject: [PATCH] Fix typo in the velociraptor interface when computing the
 linking length.

---
 src/velociraptor_interface.c | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/src/velociraptor_interface.c b/src/velociraptor_interface.c
index 7599d4ca5e..5cd4566eca 100644
--- a/src/velociraptor_interface.c
+++ b/src/velociraptor_interface.c
@@ -483,11 +483,13 @@ void velociraptor_invoke(struct engine *e, const int linked_with_snap) {
                   MPI_COMM_WORLD);
 #endif
 
+    const double Omega_m = e->cosmology->Omega_m;
+    const double Omega_b = e->cosmology->Omega_b;
+
     /* Linking length based on the mean DM inter-particle separation
      * in the zoom region and assuming the mean density of the Universe
      * is used in the zoom region. */
-    double mean_matter_density =
-        e->cosmology->Omega_m * e->cosmology->critical_density_0;
+    double mean_matter_density;
     if (s->with_hydro)
       mean_matter_density = (Omega_m - Omega_b) * critical_density_0;
     else
-- 
GitLab