diff --git a/Fig7_zones_below_lsmooth.py b/Fig7_zones_below_lsmooth.py index 51ea17dc3917e1b29115503cd30d67f64ec5cf6e..d0794ed55aef220cc11920bc1613ee861ed8fc47 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)