Merge pull request #92 from AviFriedlander/unitFix

Fix velocity_units in freefall.py
grackle-project · Feb 16, 2022 · a4e05a4 · a4e05a4
2 parents a45c3e2 + dba43cb
commit a4e05a4
Showing 1 changed file with 1 addition and 3 deletions.
diff --git a/src/python/examples/freefall.py b/src/python/examples/freefall.py
@@ -66,9 +66,7 @@
     my_chemistry.density_units  = mass_hydrogen_cgs # rho = 1.0 is 1.67e-24 g
     my_chemistry.length_units   = cm_per_mpc        # 1 Mpc in cm
     my_chemistry.time_units     = sec_per_Myr       # 1 Myr in s
-    my_chemistry.velocity_units = my_chemistry.a_units * \
-        (my_chemistry.length_units / my_chemistry.a_value) / \
-        my_chemistry.time_units
+    my_chemistry.velocity_units = my_chemistry.length_units / my_chemistry.time_units
 
     # set initial density and temperature
     initial_temperature = 50000. # start the gas at this temperature