# fig, ax = pyplot.subplots(nrows=1, ncols=1, figsize=(10, 13))
fig, ax = plt.subplots(nrows=1, ncols=1, figsize=(10, 10))
# kxrange = [3.0, 130.]
kxrange = [3.0, 50.0]
# ann_xvals = [41.0, 50.0]
ann_xvals = [35.0, 42.0]
spec_name = "IRS"
norm_wave_range = [6.0, 10.0]
# ann_wave_range = [15.0, 18.0]
col_vals = ["b", "g", "r", "m", "c", "y"]
modnames = ["BC15000g400v2", "BC15000g175v10", "BC30000g400v2", "BC30000g300v10"]
path = "/home/kgordon/Python_git/extstar_data/Models/"
plot_mir_set(ax, modnames, ann_xvals=[9.0, 9.0], subpath="", path=path)
ax.plot([3.0, 40.0], [1.0, 1.0], "k-")
ax.set_yscale("linear")
ax.set_ylim(0.5, 9.0)
ax.set_xscale("log")
ax.set_xlim(kxrange)
ax.set_xlabel(r"$\lambda$ [$\mu m$]", fontsize=1.3 * fontsize)
ax.set_ylabel(
r"$\lambda^4 F(\lambda)/F(8 \mu m)$ + offset", fontsize=1.3 * fontsize
)
ax.tick_params("both", length=10, width=2, which="major")
ax.tick_params("both", length=5, width=1, which="minor")
"hd096042",
"hd197702",
"hd149404",
"hd169454",
"hd229059",
"hd166734",
"hd206773",
"vicyg5",
"hd034921",
"hd152408",
]
plot_mir_set(
ax,
starnames,
extra_off_val=0.0,
ann_xvals=[9.0, 9.0],
ann_wave_range=[6.0, 10.0] * u.micron,
ann_offset=0.1,
col_vals=col_vals,
)
ann_set(
ax,
fontsize,
[37.0, 40.0],
[9.4, 0.8],
[45.0, 42.0],
"Unusable for Extinction",
"(ordered by wind strength)",
)
# col_vals = ['k', 'k', 'k', 'k', 'k', 'k']
starnames = [
"hd281159",
"hd283809",
"hd147933",
"hd029309",
"hd147889",
"hd204827",
"cygob2-1",
]
plot_mir_set(
ax,
starnames,
ann_xvals=[6.0, 6.0] * u.micron,
ann_wave_range=[5.0, 6.0] * u.micron,
ann_rot=2.5,
ann_offset=0.1,
col_vals=col_vals,
)
ann_set(
ax,
fontsize,
[37.0, 40.0],
[4.4, 0.9],
[45.0, 42.0],
"Main Sequence",
"(ordered by spectral type)",
)