set output './stellar_prop_c.ps' 
set term postscript color enhanced 'Helvetica,18'  
set grid 
set logscale y 
set key top right 
set ylabel 'Solar units' 
set xlabel 'Spectral Type' 
set xtics ('O5' 1, 'O8' 2, 'B0' 3, 'B3' 4, 'B5' 5, 'A0' 6, 'A5' 7, 'G2' 8) 
set title 'Photons from stars of different spectral types ' 
set label 1 '100 percent ionizing' at 6,3.0e10 font 'Helvetica,14' tc rgb 'blue' 
plot [0:9][0.1:1.0e16] './compare_nuv.dat'    using 1:($11/4.44E44)  with linespoints ls 1 t 'Photon/s'  ,    './compare_nuv.dat'    using 1:(1e10*($12/$11))  with linespoints ls 3 pt 3 ps 1.5 t 'Fraction ionizing'  ,    './compare_nuv.dat'    using 1:($12/1.96E35) with linespoints ls 2 pt 6 ps 2 t 'Ionizing photon/s'  ,    1e10 with lines ls 3 dt 2 t ''  
set output 
quit 
