#
# GNUPlot start creation file for Geochemistry Topic: The Upper Crust
# - Eric Thern, Feb. 08, 2002
# usage:  # gnuplot < setup-plots
#

#
# RPD
#
reset
set label "SiO2" at 1, -4
set label "TiO2" at 2, -4
set label "Al2O3" at 3, -4
set label "Fe2O3" at 4, -4
set label "FeO" at 5, -4
set label "MnO" at 6, -4
set label "MgO" at 7, -4
set label "CaO" at 8, -4
set label "Na2O" at 9, -4
set label "K2O" at 10, -4
set label "P2O5" at 11, -4
set xlabel "Oxide"
set ylabel "Relative Percent Difference (%)"
set title "RPD % of Oxides in samples of upper crust"
set data style linespoints
set terminal postscript eps color
set output "RPD.eps"
plot 'PAAS-PLOT1.csv' using ($13):($12) title "PAAS", 'RIVERS-PLOT1.csv' using ($10):($9) title "Major Rivers", 'PLUTONICS-PLOT1.csv' using ($10):($9) title "Igneous Plutonics", 'VOLCANICS-PLOT1.csv' using ($10):($9) title "Igneous Volcanics", 'SHIELD-PLOT1.csv' using ($10):($9) title "Precambrian Shield"

#
# standard deviation
#
reset
set label "SiO2" at 1, -0.5
set label "TiO2" at 2, -0.5
set label "Al2O3" at 3, -0.5
set label "Fe2O3" at 4, -0.5
set label "FeO" at 5, -0.5
set label "MnO" at 6, -0.5
set label "MgO" at 7, -0.5
set label "CaO" at 8, -0.5
set label "Na2O" at 9, -0.5
set label "K2O" at 10, -0.5
set label "P2O5" at 11, -0.5
set xlabel "Oxide"
set ylabel "Standard Deviation"
set title "Standard Deviation of Oxides in samples of upper crust"
set data style linespoints
set terminal postscript eps color
set output "standard-deviation.eps"
plot 'PAAS-PLOT1.csv' using ($13):($11) title "PAAS", 'RIVERS-PLOT1.csv' using ($10):($8) title "Major Rivers", 'PLUTONICS-PLOT1.csv' using ($10):($8) title "Igneous Plutonics", 'VOLCANICS-PLOT1.csv' using ($10):($8) title "Igneous Volcanics", 'SHIELD-PLOT1.csv' using ($10):($8) title "Precambrian Shield"










#
# Mean values
#
reset
set label "SiO2" at 1, 0.05
set label "TiO2" at 2, 0.05
set label "Al2O3" at 3, 0.05
set label "Fe2O3" at 4, 0.05
set label "FeO" at 5, 0.05
set label "MnO" at 6, 0.05
set label "MgO" at 7, 0.05
set label "CaO" at 8, 0.05
set label "Na2O" at 9, 0.05
set label "K2O" at 10, 0.05
set label "P2O5" at 11, 0.05
set xlabel "Oxide"
set ylabel "Mean values (logarithmic scale)"
set title "Mean values of Oxides in samples of upper crust"
set data style linespoints
set terminal postscript eps color
set output "mean.eps"
set logscale y
plot 'PAAS-PLOT1.csv' using ($13):($10) title "PAAS", 'RIVERS-PLOT1.csv' using ($10):($7) title "Major Rivers", 'PLUTONICS-PLOT1.csv' using ($10):($7) title "Igneous Plutonics", 'VOLCANICS-PLOT1.csv' using ($10):($7) title "Igneous Volcanics", 'SHIELD-PLOT1.csv' using ($10):($7) title "Precambrian Shield"


#
# REE for three samples
#
reset
set label "La" at 1, 0.05
set label "Ce" at 2, 0.05
set label "Pr" at 3, 0.05
set label "Nd" at 4, 0.05
set label "Sm" at 5, 0.05
set label "Eu" at 6, 0.05
set label "Gd" at 7, 0.05
set label "Tb" at 8, 0.05
set label "Dy" at 9, 0.05
set label "Ho" at 10, 0.05
set label "Er" at 11, 0.05
set label "Yb" at 12, 0.05
set xlabel "Element"
set ylabel "REE Content"
set title "REE's for three samples (2.5, 0.85, 0.27)"
set data style linespoints
set terminal postscript eps color
set output "REE.eps"
set logscale y
plot 'REE.csv' using ($13):($2) title "sample 2.5", 'REE.csv' using ($13):($5) title "sample 0.85", 'REE.csv' using ($13):($8) title "sample 0.27", 'REE.csv' using ($13):($12) title "RPD Values"












#
# RIVER and PAAS La/Th ratio vs. Th/Sc ratio
# 
reset
set label "2.5" at 5, -0.3
set label "1.7" at 10, -0.3
set label "1.7" at 15, -0.3
set label "0.85" at 20, -0.3
set label "0.44" at 25, -0.3
set label "0.32" at 30, -0.3
set label "0.27" at 35, -0.3
set label "0.22" at 40, -0.3
set label "Amazon" at 8, -0.5
set label "Congo" at 16, -0.5
set label "Ganges" at 24, -0.5
set label "MacKenzie" at 32, -0.5
set label "Mekong" at 40, -0.5
set xlabel "Sample #"
set ylabel "Ratio"
set title "La/Th ratio vs. Th/Sc ratio in PAAS and river deposits"
set data style linespoints
set terminal postscript eps color
set output "ratios.eps"
plot 'REE-RIVER-RATIO.csv' using ($3):($1) title "RIVER La/Th Ratio", 'REE-RIVER-RATIO.csv' using ($3):($2) title "RIVER Th/Sc Ratio", 'REE-PAAS-RATIO.csv' using ($3):($1) title "PAAS La/Th Ratio", 'REE-PAAS-RATIO.csv' using ($3):($2) title "PAAS Th/Sc Ratio"

#
# Plutonics vs. Silica
#
reset
set label "Granite" at 1, 30
set label "Granodiorite" at 2, 30
set label "Tonalite" at 3, 30
set label "Diorite" at 4, 30
set label "Gabbro" at 5, 30
set xlabel "Igneous Rock Name"
set ylabel "Percentage"
set title "Plutonic Igneous Rocks"
set data style linespoints
set terminal postscript eps color
set output "plutonics1.eps"
set logscale y
plot 'PLUTONICS-vs-silica.csv' using ($7):($1) title "Si02", 'PLUTONICS-vs-silica.csv' using ($7):($2) title "Fe2O3", 'PLUTONICS-vs-silica.csv' using ($7):($3) title "FeO", 'PLUTONICS-vs-silica.csv' using ($7):($4) title "MgO", 'PLUTONICS-vs-silica.csv' using ($7):($5) title "CaO", 'PLUTONICS-vs-silica.csv' using ($7):($6) title "K2O"

#
# Volcanics vs. Silica
#
reset
set label "Rhyolite" at 1, 30
set label "Rhyodacite" at 2, 30
set label "Dacite" at 3, 30
set label "Andesite" at 4, 30
set label "Basalt" at 5, 30
set xlabel "Igneous Rock Name"
set ylabel "Percentage"
set title "Volcanic Igneous Rocks"
set data style linespoints
set terminal postscript eps color
set output "volcanics1.eps"
set logscale y
plot 'VOLCANICS-vs-silica.csv' using ($7):($1) title "Si02", 'VOLCANICS-vs-silica.csv' using ($7):($2) title "Fe2O3", 'VOLCANICS-vs-silica.csv' using ($7):($3) title "FeO", 'VOLCANICS-vs-silica.csv' using ($7):($4) title "MgO", 'VOLCANICS-vs-silica.csv' using ($7):($5) title "CaO", 'VOLCANICS-vs-silica.csv' using ($7):($6) title "K2O"