From 3d93710213d4524698b9a8c86f42d0bd3919de10 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:46:02 +0000
Subject: [PATCH] Update README.md

---
 README.md | 30 +++++++++++++++++++++++++++++-
 1 file changed, 29 insertions(+), 1 deletion(-)

diff --git a/README.md b/README.md
index 6f99379..1c224e5 100644
--- a/README.md
+++ b/README.md
@@ -51,7 +51,7 @@ 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. 
 
-### Fig. 5: 
+### Fig. 5: Gravitational stability at the length scale of a smoothing kernel (constant softening)
 Reproduce Fig. 5 with:
 
 ```
@@ -66,6 +66,34 @@ mB_default = 1
 lsoft = [10., 20., 50., 100., 200., 500.] # pc
 lsoft_default = 3
 ```
+Here, `mB_default` and `lsoft_default` are the indices for the default values for each panel. 
+
+### Fig. 6: Gravitational stability at the length scale of a smoothing kernel (adaptive softening)
+Reproduce Fig. 6 with:
+
+```
+python3 Fig6_zones_at_lsmooth_adap.py
+```
+
+**Comment:**
+Change the plot for different resolution parameters by modifying the lines:
+```
+######### Plot nr 1. ########
+
+lsoftmin = 2.25
+lsmoothmin = 2.25
+mB = 6.4e4 
+```
+to update the left panel and 
+```
+######### Plot nr 2. ########
+lsoftmin = 108.
+lsmoothmin = 0.
+mB = 8.5e4
+```
+to update the right panel.
+
+### Fig. 7: Gravitational stability within a smoothing kernel 
 
 # Authors and acknowledgment
 S. Ploeckinger, University of Vienna
-- 
GitLab