# plot erodstats.out
set term pdf
set output 'erodstats.pdf'
set title 'Lansing - L, Saginawtri - S, Seed 0-4'
set xlabel 'Erosion (tn/ac)'
set ylabel 'Number of Events (#/yr)'
plot 'erodstats.out' index 0 using 3:2:1 with labels font "Helvetica,4" notitle, 'erodstats.out' index 0 using 3:2 with points title '15 year', \
     'erodstats.out' index 1 using 3:2:1 with labels font "Helvetica,4" notitle, 'erodstats.out' index 1 using 3:2 with points title '20 year', \
     'erodstats.out' index 2 using 3:2:1 with labels font "Helvetica,4" notitle, 'erodstats.out' index 2 using 3:2 with points title '30 year', \
     'erodstats.out' index 3 using 3:2:1 with labels font "Helvetica,4" notitle, 'erodstats.out' index 3 using 3:2 with points title '100 year', \
     'erodstats.out' index 4 using 3:2:1 with labels font "Helvetica,4" notitle, 'erodstats.out' index 4 using 3:2 with points title '1000 year'

set title
set logscale x
set xrange [:1500]
set xlabel 'Number of simulation years'

set ylabel 'Erosion events (#/yr)'
plot 'erodstats.sum' index 0 using 1:3:2:4 with errorbars lt 1 title 'Lansing', \
     'erodstats.sum' index 0 using 1:2 with lines lt 1 notitle, \
     'erodstats.sum' index 0 using 1:3 with lines lt 1 notitle, \
     'erodstats.sum' index 0 using 1:4 with lines lt 1 notitle, \
     'erodstats.sum' index 1 using 1:3:2:4 with errorbars lt 2 title 'Saginawtri', \
     'erodstats.sum' index 1 using 1:2 with lines lt 2 notitle, \
     'erodstats.sum' index 1 using 1:3 with lines lt 2 notitle, \
     'erodstats.sum' index 1 using 1:4 with lines lt 2 notitle

set ylabel 'Erosion Loss (Tn/ac/yr)'
plot 'erodstats.sum' index 0 using 1:6:5:7 with errorbars lt 1 title 'Lansing', \
     'erodstats.sum' index 0 using 1:5 with lines lt 1 notitle, \
     'erodstats.sum' index 0 using 1:6 with lines lt 1 notitle, \
     'erodstats.sum' index 0 using 1:7 with lines lt 1 notitle, \
     'erodstats.sum' index 1 using 1:6:5:7 with errorbars lt 2 title 'Saginawtri', \
     'erodstats.sum' index 1 using 1:5 with lines lt 2 notitle, \
     'erodstats.sum' index 1 using 1:6 with lines lt 2 notitle, \
     'erodstats.sum' index 1 using 1:7 with lines lt 2 notitle

set ylabel 'Precipitation (in/yr)'
plot 'erodstats.sum' index 0 using 1:9:8:10 with errorbars lt 1 title 'Lansing', \
     'erodstats.sum' index 0 using 1:8 with lines lt 1 notitle, \
     'erodstats.sum' index 0 using 1:9 with lines lt 1 notitle, \
     'erodstats.sum' index 0 using 1:10 with lines lt 1 notitle, \
     'erodstats.sum' index 1 using 1:9:8:10 with errorbars lt 2 title 'Saginawtri', \
     'erodstats.sum' index 1 using 1:8 with lines lt 2 notitle, \
     'erodstats.sum' index 1 using 1:9 with lines lt 2 notitle, \
     'erodstats.sum' index 1 using 1:10 with lines lt 2 notitle

set ylabel 'Yield (lbs/ac)'
plot 'erodstats.sum' index 0 using 1:12:11:13 with errorbars lt 1 title 'Lansing', \
     'erodstats.sum' index 0 using 1:11 with lines lt 1 notitle, \
     'erodstats.sum' index 0 using 1:12 with lines lt 1 notitle, \
     'erodstats.sum' index 0 using 1:13 with lines lt 1 notitle, \
     'erodstats.sum' index 1 using 1:12:11:13 with errorbars lt 2 title 'Saginawtri', \
     'erodstats.sum' index 1 using 1:11 with lines lt 2 notitle, \
     'erodstats.sum' index 1 using 1:12 with lines lt 2 notitle, \
     'erodstats.sum' index 1 using 1:13 with lines lt 2 notitle