def main():
"""
Produces the figure and saves it as a PDF.
"""
plt.clf()
axes = setup_axes()
visuals.plot_output_4axes(axes,
"../../simulations/sudden_5Gyr_5e9Msun_schmidt",
"crimson", "O")
visuals.plot_output_4axes(
axes, "../../simulations/sudden_5Gyr_5e9Msun_ts1p0_schmidt",
"deepskyblue", "O")
visuals.plot_output_4axes(axes,
"../../simulations/sudden_5Gyr_5e9Msun",
"black",
"O",
second_linestyle=':')
visuals.plot_track_points_intervals(
axes[2], vice.history("../../simulations/sudden_5Gyr_5e9Msun"))
visuals.sfr_ifr_legend(axes[0])
visuals.legend(axes[2], ["black", "crimson", "deepskyblue"], [
r"$\tau_*\propto M_\text{g}^0$ \qquad$\tau_\text{s}$ = 0",
r"$\tau_*\propto M_\text{g}^{-1/2}$\quad$\tau_\text{s}$ = 0",
r"$\tau_*\propto M_\text{g}^{-1/2}$\quad$\tau_\text{s}$ = 1 Gyr"
])
plot_ifr(axes[0], "../../simulations/sudden_5Gyr_5e9Msun_schmidt",
"crimson")
plot_ifr(axes[0], "../../simulations/sudden_5Gyr_5e9Msun_ts1p0_schmidt",
"deepskyblue")
plot_ifr(axes[0], "../../simulations/sudden_5Gyr_5e9Msun", "black")
plt.tight_layout()
visuals.yticklabel_formatter(axes[3])
plt.savefig(sys.argv[1])
plt.clf()
def legend(ax):
"""
Draws the legend denoting the oscillatory infall and efficiency models
Parameters
==========
ax :: subplot
The matplotlib subplot to put the legend on
"""
visuals.legend(ax, ["crimson", "deepskyblue"],
[r"oscillatory $\dot{M}_\text{in}$", r"oscillatory $\tau_*$"],
loc = "upper right", bbox_to_anchor = (0.99, 0.99))
def main():
"""
Produces the figure and save it as a PDF.
"""
plt.clf()
axes = setup_axes()
plot_gas_driven_models(axes[0])
plot_eff_driven_models(axes[1])
visuals.legend(axes[0][1], ["black", "crimson", "deepskyblue"],
["No Burst", "2 Gyr", "5 Gyr"])
visuals.legend(axes[1][2], ["black", "crimson", "deepskyblue"],
["No Burst", "2 Gyr", "5 Gyr"])
visuals.sfr_ifr_legend(axes[0][0])
for i in range(2):
visuals.yticklabel_formatter(axes[i][-1])
plt.tight_layout()
plt.savefig(sys.argv[1])
plt.clf()
def main():
"""
Produces the figure and saves it as a PDF.
"""
plt.clf()
axes = setup_axes()
plot_gas_driven_models(axes[0])
plot_eff_driven_models(axes[1])
visuals.sfr_ifr_legend(axes[0][0])
visuals.legend(axes[0][-1], ["black", "crimson", "deepskyblue"], [
r"$\tau_\text{s}$ = 0", r"$\tau_\text{s}$ = 0.5 Gyr",
r"$\tau_\text{s}$ = 1 Gyr"
],
loc="upper right",
bbox_to_anchor=(0.99, 0.99))
for i in range(2):
visuals.yticklabel_formatter(axes[i][-1])
plt.tight_layout()
plt.subplots_adjust(right=0.985)
plt.savefig(sys.argv[1])
plt.clf()
def main():
"""
Produces the figure and saves it as a PDF.
"""
plt.clf()
axes = setup_axes()
plot_gas_driven_sudden_onset_models(axes[0])
plot_gas_driven_prolonged_models(axes[1])
plot_eff_driven_models(axes[2])
visuals.legend(axes[0][1], ["black", "crimson", "deepskyblue"],
["No Burst", "2 Gyr", "5 Gyr"])
visuals.legend(axes[1][1], ["black", "crimson", "deepskyblue"],
["Sudden", "0.5 Gyr", "1 Gyr"])
visuals.legend(axes[2][2], ["black", "crimson", "deepskyblue"],
["No Burst", "2 Gyr", "5 Gyr"])
visuals.sfr_ifr_legend(axes[0][0])
visuals.sfr_ifr_legend(axes[1][0])
for i in range(len(axes)):
visuals.yticklabel_formatter(axes[i][-1])
plt.tight_layout()
plt.subplots_adjust(right=0.985)
plt.savefig(sys.argv[1])
plt.clf()