diff --git a/README.md b/README.md
index 018da2f9365eaaa4362e34bc08f90161ea5ac9f5..9cf1f721f94854503c4f34ec941f80580247fdc3 100644
--- a/README.md
+++ b/README.md
@@ -2,15 +2,18 @@
 
 The python routines to reproduce each figure from Ploeckinger et al. ([arXiv](https://arxiv.org/abs/2310.10721)) can be found here. For the figures showing results from the isolated galaxy simulations, the simulations have to be run first by following the instructions [here](https://gitlab.phaidra.org/softenedjeanscriteria/setup-galaxy-simulations). 
 
+**Required python packages:**
+- matplotlib
+- numpy
+- unyt
+
 ### Fig. 1: Comparison between Newtonian and softened Jeans mass
 Reproduce Fig. 1 with:
 
 ```
 python3 Fig1_Jeans_mass_comparison.py
 ```
-**Required python packages:**
-- matplotlib
-- numpy
+
 
 **Comment:**
 Easily adapt this figure to different constant Plummer-equivalent softening lengths by changing the line `eps = 100. # in pc`.
@@ -21,10 +24,6 @@ Reproduce Fig. 2 with:
 ```
 python3 Fig2_SPH_density_sketch.py
 ```
-**Required python packages:**
-- matplotlib
-- numpy
-- unyt
 
 **Comment:**
 Re-running the figures can lead to slightly different results because the particle positions are random. The minimum smoothing length can be varied by changing the line `hmin_in_pc = 1.55 * 10. # minimum smoothing length in pc`.
@@ -35,10 +34,6 @@ Reproduce Fig. 3 with:
 ```
 python3 Fig3_runaway_collapse_zone.py
 ```
-**Required python packages:**
-- matplotlib
-- numpy
-- unyt
 
 **Comment:**
 The plotted resolution parameters can be updated by changing the lines that set `epsarr`, `mBarr`, and `h_min_ratio_arr` (see comments in the python script).