set term qt 1 persist
set format y "%.5f"
set xlabel "Pool Height [cm]"
set ylabel "k_{eff}"
set key outside top center horizontal
set grid
set style fill transparent solid 0.25 noborder

set arrow from first 10.16, graph 0 to first 10.16, graph 1 nohead dt 2 lw 2 lc rgb "grey"
set arrow from first 35.16, graph 0 to first 35.16, graph 1 nohead dt 2 lw 2 lc rgb "grey"

plot "results/keff_45s_table" using 1:($3-$4):($3+$4) with filledcurves lc rgb "blue" notitle, \
    "results/keff_45s_table" using 1:3 with linespoints pt 7 ps 0.8 lw 2 lc rgb "blue" title "K_{eff} (4/4 Saturated Wick)", \
    "results/keff_35s_table" using 1:($3-$4):($3+$4) with filledcurves lc rgb "red" notitle, \
    "results/keff_35s_table" using 1:3 with linespoints pt 7 ps 0.8 lw 2 lc rgb "red" title "K_{eff} (3/4 Saturated Wick)", \
    "results/keff_25s_table" using 1:($3-$4):($3+$4) with filledcurves lc rgb "green" notitle, \
    "results/keff_25s_table" using 1:3 with linespoints pt 7 ps 0.8 lw 2 lc rgb "green" title "K_{eff} (2/4 Saturated Wick)", \
    "results/keff_15s_table" using 1:($3-$4):($3+$4) with filledcurves lc rgb "orange" notitle, \
    "results/keff_15s_table" using 1:3 with linespoints pt 7 ps 0.8 lw 2 lc rgb "orange" title "K_{eff} (1/4 Saturated Wick)", \

set term qt 2 persist

plot "results/range_table" using 1:2:3 with filledcurves lc rgb "blue" notitle, \
    "results/range_table" using 1:3 with linespoints pt 7 ps 0.8 lw 2 lc rgb "red" title "Maximum K_{eff} (Fully Saturated Wick)", \
    "results/range_table" using 1:2 with linespoints pt 7 ps 0.8 lw 2 lc rgb "green" title "Minimum K_{eff} (Unsaturated Wick)"