From 830ae5ed80214a465c4eafeaf21e5da944fc565a Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Sylvia=20Pl=C3=B6ckinger?= <sylvia.ploeckinger@univie.ac.at>
Date: Wed, 18 Oct 2023 08:39:39 +0000
Subject: [PATCH] Update README.md

---
 README.md | 20 +++++++++++++++++++-
 1 file changed, 19 insertions(+), 1 deletion(-)

diff --git a/README.md b/README.md
index 9cf1f72..a6402da 100644
--- a/README.md
+++ b/README.md
@@ -7,6 +7,9 @@ The python routines to reproduce each figure from Ploeckinger et al. ([arXiv](ht
 - numpy
 - unyt
 
+**Additional python packages required for individual plots**
+- astropy (Fig. 4) 
+
 ### Fig. 1: Comparison between Newtonian and softened Jeans mass
 Reproduce Fig. 1 with:
 
@@ -39,8 +42,23 @@ python3 Fig3_runaway_collapse_zone.py
 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).  
 
 ### Fig. 4: Runaway collapse zones examples from the literature
+Reproduce Fig. 4 with:
+
+```
+python3 Fig4_runaway_collapse_zone_examples.py 
+```
+
+**Comment:**
+The astropy python package is used to convert the over-density into a physical density for the Springel & Hernquist (2003) effective equation of state. If the relevant lines that add this line to the subplot are commeted out, the astropy package is not necessary. The parameters for some other simulation projects can be found in the python script and can be used for any of the subplots by modifying the call to `make_subplot`. Some examples are commented out under "other options" towards the end of the python script. 
 
-TBC..
+### Fig. 5: 
+Reproduce Fig. 5 with:
+
+```
+python3 Fig5_zones_at_lsmooth.py 
+```
+
+**Comment:**
 
 
 # Authors and acknowledgment
-- 
GitLab