pro W51_plotSpectraNoAvg ;plots ammiona lines on common velocity axis ;region 1 (19:23:38.4,+14.31.34.9) file1='/users/kchristi/scratch/W51spectra/11spectraNoAvg1.txt' file2='/users/kchristi/scratch/W51spectra/22spectraNoAvg1.txt' file4='/users/kchristi/scratch/W51spectra/44spectraAvg2_1.txt' file5='/users/kchristi/scratch/W51spectra/55spectraAvg2_1.txt' file6='/users/kchristi/scratch/W51spectra/66spectraAvg2_1.txt' readcol, file1, velocity1, intensity1 readcol, file2, velocity2, intensity2 readcol, file4, velocity4, intensity4 readcol, file5, velocity5, intensity5 readcol, file6, velocity6, intensity6 DEVICE, DECOMPOSED=0 loadct, 5 !P.COLOR=0 set_plot,'ps' device,/inches, xsize=4, ysize=6,$ filename='/users/kchristi/scratch/W51spectra/W51_spectralPlotNoAvg.ps' plot, velocity1, intensity1,$ xrange = [0,120],$ yrange = [-0.5,4],$ position = [0.10, 0.05, 0.9, 0.27],$ background = 255,$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ xtitle = 'Velocity km/s',$ ytitle = '(1,1)',$ ystyle = 8,$ xstyle = 8,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2] plot, velocity2, intensity2,$ xrange = [0,120],$ yrange = [-0.5,4],$ position = [0.10, 0.22, 0.9, 0.44],$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ ytitle = '(2,2)',$ ystyle = 8,$ xstyle = 4,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2],$ /NoErase plot, velocity4, intensity4,$ xrange = [0,120],$ yrange = [-0.5,4],$ position = [0.10, 0.39, 0.9, 0.61],$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ ytitle = 'K/Beam!C!C(4,4)',$ ystyle = 8,$ xstyle = 4,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2],$ /NoErase plot, velocity5, intensity5,$ xrange = [0,120],$ yrange = [-0.5,4],$ position = [0.10, 0.56, 0.9, 0.78],$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ ytitle = '(5,5)',$ ystyle = 8,$ xstyle = 4,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2],$ /NoErase plot, velocity6, intensity6,$ xrange = [0,120],$ yrange = [-0.5,4],$ position = [0.10, 0.73, 0.9, 0.95],$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ ytitle = '(6,6)',$ ystyle = 8,$ xstyle = 4,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2],$ /NoErase plots, [0.10, 0.05], /Normal ;fix left axis plots, [0.90, 0.98], /Normal ;fix right axis device, /close set_plot, 'X' ;region 2 (19:23:40.2,+14.31.11.4) -maser file11='/users/kchristi/scratch/W51spectra/11spectraNoAvg2.txt' file22='/users/kchristi/scratch/W51spectra/22spectraNoAvg2.txt' file44='/users/kchristi/scratch/W51spectra/44spectraAvg2_2.txt' file55='/users/kchristi/scratch/W51spectra/55spectraAvg2_2.txt' file66='/users/kchristi/scratch/W51spectra/66spectraAvg2_2.txt' readcol, file11, velocity11, intensity11 readcol, file22, velocity22, intensity22 readcol, file44, velocity44, intensity44 readcol, file55, velocity55, intensity55 readcol, file66, velocity66, intensity66 DEVICE, DECOMPOSED=0 loadct, 5 !P.COLOR=0 set_plot,'ps' device, filename='/users/kchristi/scratch/W51spectra/W51_spectralPlotNoAvg2.ps' plot, velocity11, intensity11,$ xrange = [0,120],$ yrange = [-0.5,5],$ position = [0.10, 0.05, 0.9, 0.27],$ background = 255,$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ xtitle = 'Velocity km/s',$ ytitle = '(1,1)',$ ystyle = 8,$ xstyle = 8,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2] plot, velocity22, intensity22,$ xrange = [0,120],$ yrange = [-0.5,4],$ position = [0.10, 0.22, 0.9, 0.44],$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ ytitle = '(2,2)',$ ystyle = 8,$ xstyle = 4,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2],$ /NoErase plot, velocity44, intensity44,$ xrange = [0,120],$ yrange = [-0.5,4],$ position = [0.10, 0.39, 0.9, 0.61],$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ ytitle = 'K/Beam!C!C(4,4)',$ ystyle = 8,$ xstyle = 4,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2],$ /NoErase plot, velocity55, intensity55,$ xrange = [0,120],$ yrange = [-0.5,4],$ position = [0.10, 0.56, 0.9, 0.78],$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ ytitle = '(5,5)',$ ystyle = 8,$ xstyle = 4,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2],$ /NoErase plot, velocity66, intensity66,$ xrange = [0,120],$ yrange = [-0.5,4],$ position = [0.10, 0.73, 0.9, 0.95],$ charthick = 4,$ thick = 4,$ xthick = 4,$ ythick = 4,$ color = 0,$ ytitle = '(6,6)',$ ystyle = 8,$ xstyle = 4,$ yticks = 3,$ yminor = 0,$ ytickv = [0,1,2],$ /NoErase plots, [0.10, 0.05], /Normal ;fix left axis plots, [0.90, 0.98], /Normal ;fix right axis device, /close set_plot, 'X' end