diff --git a/examples/MHDTests/MinimalResMagnetisedCloudCollapse/plotSolution.py b/examples/MHDTests/MinimalResMagnetisedCloudCollapse/plotSolution.py
index 3c99017a361c1fba79841b987978b515a1be4d66..c1f024e83c1a8c6ef98b0bf075e6a3e488aef5b2 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)