linestyle='-',
color='r',
linewidth=1,
alpha=1,
)
axes[1, 0].set_xscale('log')
axes[1, 0].set_yscale('log')
axes[1, 0].set_ylabel(r'$f_{\rm gas}(<r) = M_{\rm gas} / (M_{500}\ f_b)$')
axes[1, 0].set_xlabel(r'$r/r_{500}$')
axes[1, 0].set_ylim([1e-5, 1])
axes[1, 0].set_xlim([0.01, 1])
axes[1, 1].plot(
radial_bin_centres,
entropy_profile * gas_fraction_enclosed**(2 / 3) *
cosmology.ez_function(gas_profile_obj.z)**(2 / 3),
linestyle='-',
color='r',
linewidth=1,
alpha=1,
)
axes[1, 1].set_xscale('log')
axes[1, 1].set_yscale('log')
axes[1, 1].set_ylabel(
r'$E(z) ^ {2/3}~(K/K_{500})~\times~f_{\rm gas}(<r) ^ {2/3}$')
axes[1, 1].set_xlabel(r'$r/r_{500}$')
axes[1, 1].set_ylim([5e-5, 3])
axes[1, 1].set_xlim([0.01, 1])
sun_observations = Sun2009()
r_r500, S_S500_50, S_S500_10, S_S500_90 = sun_observations.get_shortcut()
)
axes[1, 1].plot(
radial_bin_centres,
entropy_profile,
**line_style
)
axes[2, 0].plot(
radial_bin_centres,
gas_fraction_enclosed,
**line_style
)
axes[2, 1].plot(
radial_bin_centres,
entropy_profile * gas_fraction_enclosed ** (2 / 3) * cosmology.ez_function(redshift) ** (2 / 3),
**line_style
)
delete_last_line()
axes[0, 0].remove()
axes[0, 1].remove()
top_row_axes = fig.add_subplot(gs[0, :])
cax = top_row_axes.inset_axes([0.25, 0, 0.5, 1], transform=top_row_axes.transAxes)
top_row_axes.axis('off')
colorbar = mpl.colorbar.ColorbarBase(
cax,
cmap=mpl.cm.jet,
norm=mpl.colors.LogNorm(vmin=t_bounds[0], vmax=t_bounds[1]),
orientation='horizontal'