# Plot of core mass 

set style line 1 lc rgb '#0060ad' lt 1 lw 3 pt 5   
set style line 2 lc rgb '#dd181f' lt 1 lw 3 pt 7
set style line 3 lc rgb '#9922ab' lt 1 lw 3 pt 7
set style line 4 lc rgb '#44af22' lt 1 lw 3 pt 7
set style line 5 lc rgb '#ffaa00' lt 1 lw 3 pt 7
set style line 6 lc rgb '#00aaff' lt 1 lw 3 pt 7
set style line 7 lc rgb '#00abab' lt 1 lw 3 pt 7
set style line 8 lc rgb '#000000' lt 1 lw 3 pt 8

set st da l
set xlabel "Time (t_{cc})"
set ylabel "m_{core} (m_{c})"
set xrange [0:50]
set xtics 10

set bmargin 3.5
set tmargin 1.
set rmargin 1.5
set lmargin 9

#set key top left Left
set key top right
unset label
set label "a)" at graph 0.02, graph 0.95

set yrange [0:1.2]


plot "m10c1p1_sphere_2drz_r128_stats3.dat" using 2:($7/0.13323) ls 1 title 'c1shock'
replot "m10c1_wind1_2drz_r128_stats3.dat" using 2:($7/0.13269) ls 2 title 'c1wind1'
replot "m10c1_wind1A_2drz_r128_stats3.dat" using ($1/0.074):($7/0.13247) ls 3 title 'c1wind1a'
replot "m10c1_wind1B_2drz_r128_stats3.dat" using ($1/0.023):($7/0.13218) ls 7 title 'c1wind1b'
replot "m10c1_wind1C_2drz_r128_stats3.dat" using ($1/0.0074):($7/0.13168) ls 6 title 'c1wind1c'

#plot "m10c1b1l2o45r32_bx_stats3.dat" using (($4-3.13)/1.357):($8/0.2936) ls 1 title 'l2'
#time is less the amount before zero, divided by tcs/tcc - don’t need tcs here

set term postscript eps enhanced color
set output "mcore_wind.eps"
set size 0.6,0.6
replot