From 08f308724eb1dbddb34bf8aaf830c11046b7b231 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 15:39:31 +0200
Subject: [PATCH] Minor updates to Fig 7

---
 Fig7_zones_below_lsmooth.py | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/Fig7_zones_below_lsmooth.py b/Fig7_zones_below_lsmooth.py
index 51ea17d..d0794ed 100644
--- a/Fig7_zones_below_lsmooth.py
+++ b/Fig7_zones_below_lsmooth.py
@@ -4,7 +4,6 @@ import matplotlib.pyplot as plt
 from unyt import msun, pc
 from unyt import proton_mass_cgs as mH
 import matplotlib.pylab as pl
-from astropy.cosmology import Planck18
 import matplotlib.patches as mpatches
 from matplotlib.collections import PatchCollection
 
@@ -194,7 +193,6 @@ def get_title(eps, hmin):
     return title
 
 def make_one_plot(outputname):
-    mB   = 1.e5
     fig, ax = plt.subplots(nrows=1, ncols=2, figsize = (3.5*1.3*1.8,3.5))
     fig.suptitle("Gravitational stability below the hydro spatial resolution limit l$_{\mathrm{smooth}}$")
     fig.subplots_adjust(bottom = 0.18, top = 0.8, left = 0.07, right = 0.98)
@@ -204,9 +202,13 @@ def make_one_plot(outputname):
 
     hmin =  0.2 / kernel_support_over_smoothing_length
 
+    ######### Plot nr 1 ###############
+    mB   = 1.e5
     epsarr = [100.]
     im = make_subplot(ax[0], hmin, mB, epsarr)
 
+    ######### Plot nr 2 ###############
+    mB   = 1.e5
     epsarr = [10., 20., 50., 100., 200., 500.]
     im = make_subplot(ax[1], hmin, mB, epsarr)
 
-- 
GitLab