set terminal postscript eps color enhanced solid font "Helvetica,20" size 13cm,17cm
set colorsequence classic
set output 'mollow.eps'
set multiplot layout 4, 1

infile0='data_spectrum/80amp'
infile0_s='sim_spectrum/80amp.dat'
infile1='data_spectrum/195amp'
infile1_s='sim_spectrum/195amp.dat'
infile2='data_spectrum/265amp'
infile2_s='sim_spectrum/265amp.dat'
infile3='data_spectrum/350amp'
infile3_s='sim_spectrum/350amp.dat'
infile4='data_spectrum/455amp'
infile4_s='sim_spectrum/455amp.dat'

unset ylabel
unset xlabel
unset key

#
set yrange [0:1.1]
set xrange [-90:90]
set ytics 0.25 nomirror
set xtics 20 nomirror
set format x ""
set lmargin 8
set rmargin 1
set tmargin at screen 0.98
set bmargin at screen 0.82
set label 1 at -87, 0.98
set label 1 "(a)" tc lt 7
plot infile0 u ($1-85):($2/0.03968):(sqrt($3)/($6*0.03968)) w yerr lc 1 ps 0.9 pt 7 t '195amp_data',\
infile0_s u ($1-85):($2) w lines lw 3 lc 3
#
unset label 1
set label 1 at -87, 0.98
set label 1 "(b)" tc lt 7
set label 2 at -80, 0.6
set label 2 "{/Symbol W} = 2{/Symbol p} x 16.6(3) MHz" tc lt 7
set tmargin at screen 0.80
set bmargin at screen 0.63
plot infile1 u ($1-134):($2/0.04561):(sqrt($3)/($6*0.04561)) w yerr lc 1 ps 0.9 pt 7 t '195amp_data',\
infile1_s u ($1-134):($2) w lines lw 3 lc 3
#
unset label 1
unset label 2
set label 1 at -87, 0.98
set label 1 "(c)" tc lt 7
set label 2 at -80, 0.6
set label 2 "{/Symbol W} = 2{/Symbol p} x 28.3(3) MHz" tc lt 7
set tmargin at screen 0.61
set bmargin at screen 0.44
set ylabel "Normalized Intensity" offset 1,0
plot infile2 u ($1-136):($2/0.03514):(sqrt($3)/($6*0.03514)) w yerr lc 1 ps 0.9 pt 7 t '5 pi',\
infile2_s u ($1-136):($2) w lines lw 3 lc 3
#
unset label 1
unset label 2
set label 1 at -87, 0.98
set label 1 "(d)" tc lt 7
set label 2 at -80, 0.6
set label 2 "{/Symbol W} = 2{/Symbol p} x 45.7(3) MHz" tc lt 7
set tmargin at screen 0.42
set bmargin at screen 0.25
unset ylabel
plot infile3 u ($1-134):($2/0.03078):(sqrt($3)/($6*0.03078)) w yerr lc 1 ps 0.9 pt 7 t '15 pi',\
infile3_s u ($1-134):($2) w lines lw 3 lc 3
#
unset label 1
unset label 2
set label 1 at -87, 0.98
set label 1 "(e)" tc lt 7
set label 2 at -80, 0.6
set label 2 "{/Symbol W} = 2{/Symbol p} x 64.1(2) MHz" tc lt 7
set format x "%.0f"
set tmargin at screen 0.23
set bmargin at screen 0.06
set xlabel "Frequency {/Symbol w}/2{/Symbol p} (MHz)" offset 0,0.4
plot infile4 u ($1-134):($2/0.03300):(sqrt($3)/($6*0.03300)) w yerr lc 1 ps 0.9 pt 7 t '15 pi',\
infile4_s u ($1-134):($2) w lines lw 3 lc 3
unset multiplot
pause -1