#
# GNUPlot start creation file
# Eric Thern
# usage:  # gnuplot < setupplots
#

#
# major-oxides mean plot
#
reset
set label "SiO2" at 1, 0.04
set label "Al2O3" at 2, 0.04
set label "Fe2O3" at 3, 0.04
set label "MgO" at 4, 0.04
set label "CaO" at 5, 0.04
set label "Na2O" at 6, 0.04
set label "K2O" at 7, 0.04
set label "TiO2" at 8, 0.04
set label "P2O5" at 9, 0.04
set label "MnO" at 10, 0.04
set xlabel "Major Oxides in Peach Bottom Slate"
set ylabel "Mean %"
set title "Mean of the major oxides"
set grid
set data style linespoints
set terminal postscript eps color
set output "major-oxides-pbs.eps"
set logscale y
plot 'major-oxides' using 1:2 title "Mean of Major Oxides in Peach Bottom Slate"

#
# PAAS mean plot
#
reset
set label "SiO2" at 1, 0.04
set label "Al2O3" at 2, 0.04
set label "Fe2O3" at 3, 0.04
set label "MgO" at 4, 0.04
set label "CaO" at 5, 0.04
set label "Na2O" at 6, 0.04
set label "K2O" at 7, 0.04
set label "TiO2" at 8, 0.04
set xlabel "Major Oxides in PAAS"
set ylabel "Concentration %"
set title "Major Oxides in PAAS"
set grid
set data style linespoints
set terminal postscript eps color
set output "major-oxides-paas.eps"
set logscale y
plot 'major-oxides' using 1:3 title "Major Oxides in PAAS"


#
# PAAS and PB plot
#
reset
set label "SiO2" at 1, 0.04
set label "Al2O3" at 2, 0.04
set label "Fe2O3" at 3, 0.04
set label "MgO" at 4, 0.04
set label "CaO" at 5, 0.04
set label "Na2O" at 6, 0.04
set label "K2O" at 7, 0.04
set label "TiO2" at 8, 0.04
set label "P2O5" at 9, 0.04
set label "MnO" at 10, 0.04
set xlabel "Major Oxides in PAAS and PBS"
set ylabel "Concentration %"
set title "Major Oxides in PAAS and PBS"
set grid
set data style linespoints
set terminal postscript eps color
set output "major-oxides.eps"
set logscale y
plot 'major-oxides' using 1:3 title "Major Oxides in PAAS",\
 'major-oxides' using 1:2 title "Major Oxides in PBS"



#
# PAAS and PB plot with errorbars
#
reset
set label "SiO2" at 1, 0.0025
set label "Al2O3" at 2, 0.0025
set label "Fe2O3" at 3, 0.0025
set label "MgO" at 4, 0.0025
set label "CaO" at 5, 0.0025
set label "Na2O" at 6, 0.0025
set label "K2O" at 7, 0.0025
set label "TiO2" at 8, 0.0025
set label "P2O5" at 9, 0.0025
set label "MnO" at 10, 0.0025
set xlabel "Major Oxides in PAAS and PBS"
set ylabel "Concentration %"
set title "Major Oxides in PAAS and PBS including errorbars"
set grid
set data style linespoints
set terminal postscript eps color
set output "major-oxides-errorbars.eps"
set logscale y
plot 'major-oxides' using 1:3 title "Major Oxides in PAAS",\
 'major-oxides' using 1:2 title "Major Oxides in PBS",\
 'major-oxides' using 1:2:4 title "Major Oxides in PBS Errorbars" with errorbars




#
# REE normalization to chondritic values & PAAS
#
reset
set label "SiO2" at 1, 0.0025
set label "Al2O3" at 2, 0.0025
set label "Fe2O3" at 3, 0.0025
set label "MgO" at 4, 0.0025
set label "CaO" at 5, 0.0025
set label "Na2O" at 6, 0.0025
set label "K2O" at 7, 0.0025
set label "TiO2" at 8, 0.0025
set label "P2O5" at 9, 0.0025
set label "MnO" at 10, 0.0025
set xlabel "REE Elements"
set ylabel "Normalized Concentration of REE's to Chondritic values"
set title "Normalized REE and PAAS values"
set grid
set data style linespoints
set terminal postscript eps color
set output "lichen-normalized.eps"
set logscale y
plot 'lichen' using 1:2 title "Gabbro",\
 'lichen' using 1:3 title "Quartzite",\
 'lichen' using 1:4 title "Granodiorite",\
 'lichen' using 1:5 title "Dolostone",\
 'lichen' using 1:6 title "Slate",\
 'lichen' using 1:7 title "Sand",\
 'lichen' using 1:8 title "Andesite",\
 'lichen' using 1:9 title "PAAS"




#
# REE's normalization to chondritic values & PAAS -- REGRESSION LINE
#
reset
set xlabel "Chondrite Normalized PAAS REE's"
set ylabel "Chondrite Normalized Lichen REE's"
set title "Normalized REE and PAAS values with regression"
set grid
set data style linespoints
set terminal table
set output 'lichen-regression.dat'
plot 'lichen' using 9:10
set terminal postscript eps color
set output "lichen-normalized-regression.eps"
f(x) = a*x +b
fit f(x) 'lichen-regression.dat' via a,b
plot f(x), 'lichen-regression.dat' title "REE lichen vs. REE PAAS"