def plot_oscillatory_infall_model(axes, inset):
"""
Plots the models in which the infall rate oscillates with time
Parameters
==========
axes :: list
The 1-D list of matplotlib axis objects to plot on
inset :: subplot
The matplotlib axis to plot the tracks inset on
"""
visuals.plot_output_3axes(axes, "../../simulations/SFRoscil_amp0p3_per4",
"crimson", "O")
visuals.plot_output_3axes(axes, "../../simulations/SFRoscil_amp0p6_per2",
"deepskyblue", "O")
visuals.plot_output_3axes(axes, "../../simulations/SFRoscil_amp0p3_per2",
"black", "O")
visuals.plot_track_points_intervals(
axes[1], vice.history("../../simulations/SFRoscil_amp0p3_per2"))
# visuals.plot_output_3axes(axes, "../../simulations/slow", "lime", "O")
visuals.plot_inset(inset, "../../simulations/SFRoscil_amp0p3_per4",
"crimson")
visuals.plot_inset(inset, "../../simulations/SFRoscil_amp0p6_per2",
"deepskyblue")
visuals.plot_inset(inset, "../../simulations/SFRoscil_amp0p3_per2",
"black")
visuals.plot_inset(inset,
"../../simulations/default",
"black",
linestyle=":")
visuals.plot_track_points_intervals(
inset, vice.history("../../simulations/SFRoscil_amp0p3_per2"))
# visuals.plot_inset(inset, "../../simulations/slow", "lime")
visuals.plot_reference(axes)
visuals.sfr_ifr_legend(axes[0], ncol=2)
def plot_oscillatory_eff_model(axes, inset):
"""
Plots the models in which the SFE timescale oscillates with time
Parameters
==========
axes :: list
The 1-D list of matplotlib axis objects to plot on
inset :: subplot
The matplotlib axis to plot the tracks inset on
"""
visuals.plot_output_4axes(axes, "../../simulations/SFEoscil_amp0p2_per4",
"crimson", "O")
visuals.plot_output_4axes(axes, "../../simulations/SFEoscil_amp0p4_per2",
"deepskyblue", "O")
visuals.plot_output_4axes(axes, "../../simulations/SFEoscil_amp0p2_per2",
"black", "O")
visuals.plot_track_points_intervals(
axes[2], vice.history("../../simulations/SFEoscil_amp0p2_per2"))
# visuals.plot_output_4axes(axes,
# "../../simulations/slow", "lime", "O")
visuals.plot_inset(inset, "../../simulations/SFEoscil_amp0p2_per4",
"crimson")
visuals.plot_inset(inset, "../../simulations/SFEoscil_amp0p4_per2",
"deepskyblue")
visuals.plot_inset(inset, "../../simulations/SFEoscil_amp0p2_per2",
"black")
visuals.plot_inset(inset,
"../../simulations/default",
"black",
linestyle=':')
visuals.plot_track_points_intervals(
inset, vice.history("../../simulations/SFEoscil_amp0p2_per2"))
# visuals.plot_inset(inset,
# "../../simulations/slow", "lime")
visuals.plot_reference(axes)
axes[0].set_xlim([-1, 11])
axes[0].set_ylim([-1, 17])
axes[1].set_xlim([-1.7, 0.2])
axes[1].set_ylim([-0.1, 0.5])
axes[2].set_xlim([-0.1, 0.5])
axes[2].set_ylim([0.2, 50])
visuals.set_labels_3axes(axes, "O")
inset_xlim = [-0.26, -0.06]
inset_ylim = [0.06, 0.16]
visuals.draw_box(axes[1], inset_xlim, inset_ylim)
inset = visuals.zoom_box(axes[1], inset_xlim, inset_ylim, zoom=3.8)
plot_output_3axes(axes, "../../simulations/SFRoscil_amp0p3_per4",
"crimson")
plot_output_3axes(axes, "../../simulations/SFRoscil_amp0p6_per2", "blue")
plot_output_3axes(axes, "../../simulations/SFRoscil_amp0p3_per2", "black")
visuals.plot_inset(inset, "../../simulations/SFRoscil_amp0p3_per4",
"crimson")
visuals.plot_inset(inset, "../../simulations/SFRoscil_amp0p6_per2", "blue")
visuals.plot_inset(inset, "../../simulations/SFRoscil_amp0p3_per2",
"black")
visuals.plot_inset(inset,
"../../simulations/default",
"black",
linestyle=":")
visuals.plot_reference(axes)
visuals.yticklabel_formatter(axes[-1])
visuals.sfr_ifr_legend(axes[0], ncol=2)
plt.tight_layout()
plt.savefig(sys.argv[1])
plt.clf()