set terminal svg enhanced size 1250 400 fname "Nimbus Sans L" fsize 10
set output "volcanic_forcing-reconstructed.svg"
set title "Global volcanic aerosol forcing 500 BCE - 2000 CE"  font "Nimbus Sans L, 20"
set datafile missing '"#N/A"'
set key top right
set yrange [0:-40]
set ytics  0,-10,-40
set mytics 2
set xrange [-500:2000]
set xtics -500,200,2000
set mxtics 4
set xlabel 'Year'
set ylabel 'Forcing in W / m^2'
set boxwidth 1
set style fill solid
set grid mytics lt 0 lw 0.5 lc rgb "grey"
set grid mxtics lt 0 lw 0.5 lc rgb "grey"
set grid ytics lt 0 lw 0.5 lc rgb "grey"
set grid xtics lt 0 lw 0.5 lc rgb "grey"
set style line 1 lc rgb "#228B22"
set style line 2 lc rgb "red"
set style line 3 lc rgb "blue"
plot 'volcanic_forcing_t.csv' using 1:2 with boxes ls 1 title "tropical", \
     'volcanic_forcing_n.csv' using 1:2 with boxes ls 2 title "northern hemisphere", \
     'volcanic_forcing_s.csv' using 1:2 with boxes ls 3 title "southern hemisphere", \
     'volcanic_forcing_m.csv' using 1:($2-1):1 with labels title "(eruptions < -10 W/m^2 labeled with year)"