From de9dd08826d9ec80bae6df37ca531766e38ae873 Mon Sep 17 00:00:00 2001
From: Orestis Karapiperis <karapiperis@strw.leidenuniv.nl>
Date: Thu, 21 Nov 2024 13:10:52 +0100
Subject: [PATCH] Fixed mistake in free fall time calculation

---
 .../MinimalResMagnetisedCloudCollapse/plotSolution.py      | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/examples/MHDTests/MinimalResMagnetisedCloudCollapse/plotSolution.py b/examples/MHDTests/MinimalResMagnetisedCloudCollapse/plotSolution.py
index 3c99017a36..c1f024e83c 100644
--- a/examples/MHDTests/MinimalResMagnetisedCloudCollapse/plotSolution.py
+++ b/examples/MHDTests/MinimalResMagnetisedCloudCollapse/plotSolution.py
@@ -12,9 +12,12 @@ filename = sys.argv[1]
 data = load(filename)
 center = 0.5 * data.metadata.boxsize
 
-R0 = 0.015 * 3.086e18 * unyt.cm
+G = 6.6743015e-8 * unyt.cm**3 / unyt.g / unyt.s**2
+R0 = 0.015 * unyt.pc
+M = 1.0 * unyt.Msun
+
+tff = np.sqrt(2.0 * R0**3 / (G * M))
 
-tff = 3e4 * 3.156e7 * unyt.s
 tsim = data.metadata.time
 tplot = tsim / tff
 print("Showing results at %2f free fall times" % tplot)
-- 
GitLab