# python 导入数据及作图的实现

```filename="/home/res/user/csluo/test.txt"          #将文件名赋值为变量

X,Y,Z=[ ],[ ],[ ]                      #给三个空列表

with open(filename, "r") as f:               #打开文档

for i in lines:

value=[ float(s) for s in line.split( ) ]   #将 行数据 以空格分开，并以浮点型赋给value

X=append(value[0])              #将value 的第一个数值 添加到 列表X

Y=append(value[1])

Z=append(value[2])```

```a = numpy.loadtxt("a.txt")

a[:,0]```

```>>> import pandas as pd

```import matplotlib
import matplotlib.pyplot as plt
import numpy as np

##################################################################
######################################################################################### ######## 1
##################################################################

#data_a
X10,Y10= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]     #以空格分开 行 并输出 浮点数
X10.append(value[0])
Y10.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A1/xclass_spectrum_output.dat"
X11,Y11 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X11.append(value[0])
Y11.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A1/xclass_spectrum_output.dat"
X12,Y12 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X12.append(value[0])
Y12.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A1/xclass_spectrum_output.dat"
X13,Y13 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X13.append(value[0])
Y13.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A1/xclass_spectrum_output.dat"
X14,Y14 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X14.append(value[0])
Y14.append(value[1])

###################################################  data_b

filename = "/home/luoluo/Documents/python_scripts/Bdata_cut/B1_145122-145580.txt"
X15,Y15= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X15.append(value[0])
Y15.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B1/xclass_spectrum_output.dat"
X16,Y16 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X16.append(value[0])
Y16.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B1/xclass_spectrum_output.dat"
X17,Y17 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X17.append(value[0])
Y17.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B1/xclass_spectrum_output.dat"
X18,Y18 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X18.append(value[0])
Y18.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B1/xclass_spectrum_output.dat"
X19,Y19 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X19.append(value[0])
Y19.append(value[1])

###################################################################  #2 plot

plt.figure(1)
plt.subplot(211)
plt.axis([145122,145580, -1, 15])

new_ticks = np.linspace(145122,145580,10)
plt.xticks(new_ticks)

matplotlib.rcParams["xtick.direction"] = "in"
matplotlib.rcParams["ytick.direction"] = "in"
plt.plot(X10,Y10,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X11,Y11,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X12,Y12,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X13,Y13,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X14,Y14,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
#plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(145535,13,"IRAS16293A")

#plt.show()
#plt.savefig("8.png")

plt.subplot(212)
plt.axis([145122,145580, -1, 15])

new_ticks = np.linspace(145122,145580,10)
plt.xticks(new_ticks)

plt.plot(X15,Y15,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X16,Y16,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X17,Y17,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X18,Y18,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X19,Y19,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(145535,13,"IRAS16293B")

plt.plot([145263,145263],[1,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("C2H5OH", xy = (145263,1), xytext = (145280,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################
######################################################################################### ######## 2
##################################################################

#data_a
X20,Y20= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X20.append(value[0])
Y20.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A2/xclass_spectrum_output.dat"
X21,Y21 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X21.append(value[0])
Y21.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A2/xclass_spectrum_output.dat"
X22,Y22 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X22.append(value[0])
Y22.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A2/xclass_spectrum_output.dat"
X23,Y23 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X23.append(value[0])
Y23.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A2/xclass_spectrum_output.dat"
X24,Y24 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X24.append(value[0])
Y24.append(value[1])

###################################################  data_b

filename = "/home/luoluo/Documents/python_scripts/Bdata_cut/B2.txt"
X25,Y25= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X25.append(value[0])
Y25.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B2/xclass_spectrum_output.dat"
X26,Y26 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X26.append(value[0])
Y26.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B2/xclass_spectrum_output.dat"
X27,Y27 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X27.append(value[0])
Y27.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B2/xclass_spectrum_output.dat"
X28,Y28 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X28.append(value[0])
Y28.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B2/xclass_spectrum_output.dat"
X29,Y29 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X29.append(value[0])
Y29.append(value[1])

###################################################################  #2 plot

plt.figure(2)
plt.subplot(211)
plt.axis([145580,146050, -1, 15])

new_ticks = np.linspace(145580,146050,10)
plt.xticks(new_ticks)

matplotlib.rcParams["xtick.direction"] = "in"
matplotlib.rcParams["ytick.direction"] = "in"
plt.plot(X20,Y20,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X21,Y21,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X22,Y22,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X23,Y23,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X24,Y24,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
#plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(146000,13,"IRAS16293A")

plt.plot([145618.5,145618.5],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (145618,7), xytext = (145650,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146020,146020],[7,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146020,6), xytext = (146005,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()
#plt.savefig("8.png")

plt.subplot(212)
plt.axis([145580,146050, -1, 15])

new_ticks = np.linspace(145580,146050,10)
plt.xticks(new_ticks)

plt.plot(X25,Y25,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X26,Y26,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X27,Y27,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X28,Y28,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X29,Y29,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(146000,13,"IRAS16293B")

plt.plot([145618.5,145618.5],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (145618,7), xytext = (145650,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146020,146020],[7,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146020,6), xytext = (145990,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################
######################################################################################### ######## 3
##################################################################

#data_a
X30,Y30= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X30.append(value[0])
Y30.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A3/xclass_spectrum_output.dat"
X31,Y31 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X31.append(value[0])
Y31.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A3/xclass_spectrum_output.dat"
X32,Y32 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X32.append(value[0])
Y32.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A3/xclass_spectrum_output.dat"
X33,Y33 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X33.append(value[0])
Y33.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A3/xclass_spectrum_output.dat"
X34,Y34 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X34.append(value[0])
Y34.append(value[1])

###################################################  data_b

filename = "/home/luoluo/Documents/python_scripts/Bdata_cut/B3_146050-146507.txt"
X35,Y35= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X35.append(value[0])
Y35.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B3/xclass_spectrum_output.dat"
X36,Y36 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X36.append(value[0])
Y36.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B3/xclass_spectrum_output.dat"
X37,Y37 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X37.append(value[0])
Y37.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B3/xclass_spectrum_output.dat"
X38,Y38 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X38.append(value[0])
Y38.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B3/xclass_spectrum_output.dat"
X39,Y39 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X39.append(value[0])
Y39.append(value[1])

###################################################################  #2 plot

plt.figure(3)
plt.subplot(211)
plt.axis([146050,146507, -1, 15])

new_ticks = np.linspace(146050,146507,10)
plt.xticks(new_ticks)

matplotlib.rcParams["xtick.direction"] = "in"
matplotlib.rcParams["ytick.direction"] = "in"
plt.plot(X30,Y30,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X31,Y31,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X32,Y32,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X33,Y33,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X34,Y34,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
#plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(146465,13,"IRAS16293A")

plt.plot([146172,146172],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146172,7), xytext = (146190,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146232,146232],[3,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146232,10), xytext = (146210,5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146343,146343],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146343,7), xytext = (146320,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146355,146355],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146355,6), xytext = (146380,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()
#plt.savefig("8.png")

plt.subplot(212)
plt.axis([146050,146507, -1, 15])

new_ticks = np.linspace(146050,146507,10)
plt.xticks(new_ticks)

plt.plot(X35,Y35,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X36,Y36,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X37,Y37,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X38,Y38,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X39,Y39,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(146465,13,"IRAS16293B")

plt.plot([146343,146343],[5.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146343,7), xytext = (146320,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146355,146355],[5.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146355,6), xytext = (146380,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.plot([157324,157324],[2.5,15], color ="grey", linewidth=1, linestyle="--")
#plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################
######################################################################################### ######## 4
##################################################################

#data_a
X40,Y40= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X40.append(value[0])
Y40.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A4/xclass_spectrum_output.dat"
X41,Y41 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X41.append(value[0])
Y41.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A4/xclass_spectrum_output.dat"
X42,Y42 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X42.append(value[0])
Y42.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A4/xclass_spectrum_output.dat"
X43,Y43 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X43.append(value[0])
Y43.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A4/xclass_spectrum_output.dat"
X44,Y44 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X44.append(value[0])
Y44.append(value[1])

###################################################  data_b

filename = "/home/luoluo/Documents/python_scripts/Bdata_cut/B4_146507-146976.txt"
X45,Y45= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X45.append(value[0])
Y45.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B4/xclass_spectrum_output.dat"
X46,Y46 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X46.append(value[0])
Y46.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B4/xclass_spectrum_output.dat"
X47,Y47 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X47.append(value[0])
Y47.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B4/xclass_spectrum_output.dat"
X48,Y48 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X48.append(value[0])
Y48.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B4/xclass_spectrum_output.dat"
X49,Y49 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X49.append(value[0])
Y49.append(value[1])

###################################################################  #2 plot

plt.figure(4)
plt.subplot(211)
plt.axis([146507,146976, -1, 15])

new_ticks = np.linspace(146507,146976,10)
plt.xticks(new_ticks)

matplotlib.rcParams["xtick.direction"] = "in"
matplotlib.rcParams["ytick.direction"] = "in"
plt.plot(X40,Y40,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X41,Y41,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X42,Y42,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X43,Y43,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X44,Y44,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
#plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(146930,13,"IRAS16293A")

#plt.plot([157323,157323],[2.5,15], color ="grey", linewidth=1, linestyle="--")
#plt.annotate("CH3COCH3", xy = (157323,3), xytext = (157345,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()
#plt.savefig("8.png")

plt.subplot(212)
plt.axis([146507,146976, -1, 15])

new_ticks = np.linspace(146507,146976,10)
plt.xticks(new_ticks)

plt.plot(X45,Y45,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X46,Y46,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X47,Y47,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X48,Y48,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X49,Y49,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(146930,13,"IRAS16293B")

#plt.plot([157324,157324],[2.5,15], color ="grey", linewidth=1, linestyle="--")
#plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################
######################################################################################### ######## 5
##################################################################

#data_a
X50,Y50= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X50.append(value[0])
Y50.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A5/xclass_spectrum_output.dat"
X51,Y51 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X51.append(value[0])
Y51.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A5/xclass_spectrum_output.dat"
X52,Y52 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X52.append(value[0])
Y52.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A5/xclass_spectrum_output.dat"
X53,Y53 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X53.append(value[0])
Y53.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A5/xclass_spectrum_output.dat"
X54,Y54 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X54.append(value[0])
Y54.append(value[1])

###################################################  data_b

filename = "/home/luoluo/Documents/python_scripts/Bdata_cut/B5_156840-157303.txt"
X55,Y55= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X55.append(value[0])
Y55.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B5/xclass_spectrum_output.dat"
X56,Y56 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X56.append(value[0])
Y56.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B5/xclass_spectrum_output.dat"
X57,Y57 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X57.append(value[0])
Y57.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B5/xclass_spectrum_output.dat"
X58,Y58 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X58.append(value[0])
Y58.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B5/xclass_spectrum_output.dat"
X59,Y59 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X59.append(value[0])
Y59.append(value[1])

###################################################################  #2 plot

plt.figure(5)
plt.subplot(211)
plt.axis([156840, 157303, -1, 15])

new_ticks = np.linspace(156840,157303,10)
plt.xticks(new_ticks)

matplotlib.rcParams["xtick.direction"] = "in"
matplotlib.rcParams["ytick.direction"] = "in"
plt.plot(X50,Y50,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X51,Y51,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X52,Y52,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X53,Y53,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X54,Y54,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
#plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(157256,13,"IRAS16293A")

#plt.show()
#plt.savefig("8.png")

plt.subplot(212)
plt.axis([156820,157303, -1, 15])

new_ticks = np.linspace(156820,157303,10)
plt.xticks(new_ticks)

plt.plot(X55,Y55,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X56,Y56,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X57,Y57,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X58,Y58,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X59,Y59,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(157255,13,"IRAS16293B")

#plt.plot([157324,157324],[2.5,15], color ="grey", linewidth=1, linestyle="--")
#plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################
######################################################################################### ######## 6
##################################################################

#data_a
X60,Y60= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X60.append(value[0])
Y60.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A6/xclass_spectrum_output.dat"
X61,Y61 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X61.append(value[0])
Y61.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A6/xclass_spectrum_output.dat"
X62,Y62 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X62.append(value[0])
Y62.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A6/xclass_spectrum_output.dat"
X63,Y63 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X63.append(value[0])
Y63.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A6/xclass_spectrum_output.dat"
X64,Y64 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X64.append(value[0])
Y64.append(value[1])

###################################################  data_b

filename = "/home/luoluo/Documents/python_scripts/Bdata_cut/B6_157303-157820.txt"
X65,Y65= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X65.append(value[0])
Y65.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B6/xclass_spectrum_output.dat"
X66,Y66 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X66.append(value[0])
Y66.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B6/xclass_spectrum_output.dat"
X67,Y67 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X67.append(value[0])
Y67.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B6/xclass_spectrum_output.dat"
X68,Y68 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X68.append(value[0])
Y68.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B6/xclass_spectrum_output.dat"
X69,Y69 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X69.append(value[0])
Y69.append(value[1])

###################################################################  #2 plot

plt.figure(6)
plt.subplot(211)
plt.axis([157303, 157820, -1, 15])

new_ticks = np.linspace(157303, 157820,10)
plt.xticks(new_ticks)

matplotlib.rcParams["xtick.direction"] = "in"
matplotlib.rcParams["ytick.direction"] = "in"
plt.plot(X60,Y60,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X61,Y61,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X62,Y62,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X63,Y63,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X64,Y64,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
#plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(157765,13,"IRAS16293A")

plt.plot([157323,157323],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (157323,3), xytext = (157345,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()
#plt.savefig("8.png")

plt.subplot(212)
plt.axis([157303, 157820, -1, 15])

new_ticks = np.linspace(157303, 157820,10)
plt.xticks(new_ticks)

plt.plot(X65,Y65,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X66,Y66,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X67,Y67,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X68,Y68,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X69,Y69,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(157765,13,"IRAS16293B")

plt.plot([157324,157324],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([157505.5,157505.5],[5.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("C2H5OH", xy = (157505,6), xytext = (157450,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

###############################################################################################
################################################################################################################ 7
###############################################################################################

#data_a
X70,Y70= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X70.append(value[0])
Y70.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A7/xclass_spectrum_output.dat"
X71,Y71 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X71.append(value[0])
Y71.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A7/xclass_spectrum_output.dat"
X72,Y72 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X72.append(value[0])
Y72.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A7/xclass_spectrum_output.dat"
X73,Y73 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X73.append(value[0])
Y73.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A7/xclass_spectrum_output.dat"
X74,Y74 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X74.append(value[0])
Y74.append(value[1])
#########################################################
#data_b
filename = "/home/luoluo/Documents/python_scripts/Bdata_cut/B7_157820-158228.txt"
X75,Y75= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X75.append(value[0])
Y75.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B7/xclass_spectrum_output.dat"
X76,Y76 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X76.append(value[0])
Y76.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B7/xclass_spectrum_output.dat"
X77,Y77 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X77.append(value[0])
Y77.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B7/xclass_spectrum_output.dat"
X78,Y78 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X78.append(value[0])
Y78.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B7/xclass_spectrum_output.dat"
X79,Y79 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X79.append(value[0])
Y79.append(value[1])
#############################################################################################################################################
#2 plot
plt.figure(7)
plt.subplot(211)
plt.axis([157820, 158228, -1, 15])

new_ticks = np.linspace(157820,158228,10)
plt.xticks(new_ticks)

matplotlib.rcParams["xtick.direction"] = "in"
matplotlib.rcParams["ytick.direction"] = "in"
plt.plot(X70,Y70,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X71,Y71,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X72,Y72,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X73,Y73,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X74,Y74,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
#plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(158183,13,"IRAS16293A")

plt.plot([157894,157894],[5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (157895,7), xytext = (157880,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158120,158120],[5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158120,5), xytext = (158110,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158172,158172],[7,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158171,8), xytext = (158150,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158225,158225],[4,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158225,5), xytext = (158200,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()
#plt.savefig("8.png")

plt.subplot(212)
plt.axis([157820, 158228, -1, 15])

new_ticks = np.linspace(157820, 158228,10)
plt.xticks(new_ticks)

plt.plot(X75,Y75,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X76,Y76,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X77,Y77,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X78,Y78,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X79,Y79,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(158183,13,"IRAS16293B")

plt.plot([157895,157895],[5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (157896,7), xytext = (157880,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158121,158121],[5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158121,7), xytext = (158110,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158173,158173],[7,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158172,8), xytext = (158150,5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158226,158226],[4,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158225,5), xytext = (158205,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.plot([158389,158389],[3.5,15], color ="grey", linewidth=1, linestyle="--")
#plt.annotate("CH3OCHO v18=1", xy = (158390,5), xytext = (158400,13),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()
######################################################################################
###################################################################################### ###### 8
######################################################################################

#
#1 import data

#data_a
X80,Y80= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X80.append(value[0])
Y80.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A8/xclass_spectrum_output.dat"
X81,Y81 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X81.append(value[0])
Y81.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A8/xclass_spectrum_output.dat"
X82,Y82 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X82.append(value[0])
Y82.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A8/xclass_spectrum_output.dat"
X83,Y83 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X83.append(value[0])
Y83.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A8/xclass_spectrum_output.dat"
X84,Y84 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X84.append(value[0])
Y84.append(value[1])
#########################################################
#data_b
filename = "Bdata_cut/B8_158228-158680.txt"
X85,Y85= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X85.append(value[0])
Y85.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B8/xclass_spectrum_output.dat"
X86,Y86 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X86.append(value[0])
Y86.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B8/xclass_spectrum_output.dat"
X87,Y87 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X87.append(value[0])
Y87.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B8/xclass_spectrum_output.dat"
X88,Y88 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X88.append(value[0])
Y88.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B8/xclass_spectrum_output.dat"
X89,Y89 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X89.append(value[0])
Y89.append(value[1])
#############################################################
#2 plot
plt.figure(8)
plt.subplot(211)
plt.axis([158228, 158680, -1, 15])

new_ticks = np.linspace(158228,158680,10)
plt.xticks(new_ticks)

matplotlib.rcParams["xtick.direction"] = "in"
matplotlib.rcParams["ytick.direction"] = "in"
plt.plot(X80,Y80,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X81,Y81,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X82,Y82,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X83,Y83,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X84,Y84,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
#plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(158635,13,"IRAS16293A")

plt.plot([158378,158378],[1.9,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158378, 2.0), xytext = (158300,5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158388,158388],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158388, 4.5), xytext = (158370,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158430,158430],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158430, 3.5), xytext = (158400,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158484,158484],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158485, 4), xytext = (158450,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158519,158519],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158518,4), xytext = (158490,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158540,158540],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158541,4), xytext = (158550,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158675,158675],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158674,3), xytext = (158610,8),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()
#plt.savefig("8.png")

plt.subplot(212)
plt.axis([158228, 158680, -1, 15])

new_ticks = np.linspace(158228,158680,10)
plt.xticks(new_ticks)

plt.plot(X85,Y85,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X86,Y86,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X87,Y87,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X88,Y88,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X89,Y89,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(158635,13,"IRAS16293B")
plt.plot([158389,158389],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158390,5), xytext = (158400,13),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158484,158484],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158484,5), xytext = (158430,4),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158520,158520],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158520,4), xytext = (158500,6),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158540,158540],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158540,4), xytext = (158570,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158675,158675],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158675,4), xytext = (158640,8),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

####################################################################################################
#################################################################################################### ######## 9
####################################################################################################

#data_a

X90,Y90= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X90.append(value[0])
Y90.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A9/xclass_spectrum_output.dat"
X91,Y91 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X91.append(value[0])
Y91.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A9/xclass_spectrum_output.dat"
X92,Y92 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X92.append(value[0])
Y92.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A9/xclass_spectrum_output.dat"
X93,Y93 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X93.append(value[0])
Y93.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A9/xclass_spectrum_output.dat"
X94,Y94 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X94.append(value[0])
Y94.append(value[1])
#########################################################
#data_b
filename = "/home/luoluo/Documents/python_scripts/Bdata_cut/B9_158680-159155.txt"
X95,Y95= [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X95.append(value[0])
Y95.append(value[1])

#c2h5oh
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B9/xclass_spectrum_output.dat"
X96,Y96 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X96.append(value[0])
Y96.append(value[1])

#ch3coch3
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B9/xclass_spectrum_output.dat"
X97,Y97 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X97.append(value[0])
Y97.append(value[1])

#ch3ocho_v18=1
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B9/xclass_spectrum_output.dat"
X98,Y98 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X98.append(value[0])
Y98.append(value[1])

#ch3ocho_v=0
filename = "/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B9/xclass_spectrum_output.dat"
X99,Y99 = [],[]
with open(filename, "r") as f:
for line in lines:
value = [float(s) for s in line.split()]
X99.append(value[0])
Y99.append(value[1])

###2 plot
plt.figure(9)
#A
plt.subplot(211)
plt.axis([158680, 159155, -1, 15])

new_ticks = np.linspace(158680,159155,10)
plt.xticks(new_ticks)

matplotlib.rcParams["xtick.direction"] = "in"
matplotlib.rcParams["ytick.direction"] = "in"
plt.plot(X90,Y90,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X91,Y91,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X92,Y92,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X93,Y93,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X94,Y94,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
#plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(159110,13,"IRAS16293A")

plt.plot([158691,158691],[8,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158691, 5), xytext = (158730,3),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158701.9,158701.9],[8,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158701, 10), xytext = (158770,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([159003.5,159003.5],[3.4,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (159005, 3), xytext = (159020,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([159106,159106],[3.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (159107, 3.5), xytext = (159110,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#B
plt.subplot(212)
plt.axis([158680, 159155, -1, 15])

new_ticks = np.linspace(158680,159122,10)
plt.xticks(new_ticks)

plt.plot(X95,Y95,color="k",lw=1.0,ls="-",label="DATA")
plt.plot(X96,Y96,color="r",lw=1.0,ls="-",label="C2H5OH")
plt.plot(X97,Y97,color="g",lw=1.0,ls="-",label="CH3COCH3")
plt.plot(X98,Y98,color="cyan",lw=1.0,ls="-",label="CH3OCHO v18=1")
plt.plot(X99,Y99,color="b",lw=1.0,ls="-",label="CH3OCHO v=0")
plt.legend(loc="upper left")
plt.xlabel("Freq(MHZ)")
plt.ylabel("T(K)")
plt.text(159110,13,"IRAS16293B")

plt.plot([158689,158689],[7,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158689, 3.5), xytext = (158700,4),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.plot([158963,158963],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v=0", xy = (158694,6), xytext = (158760,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.plot([158963,158963],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v=0", xy = (158704,9.5), xytext = (158770,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158963,158963],[2.5,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158963, 3.5), xytext = (158980,6),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([159030,159030],[2.4,15], color ="grey", linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (159030,4), xytext = (159000,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

```