Added routine to plot polarization as a function of intensity.
Hard coded CO(2-1) and SiO(5-4) frequencies.