La-Y-Nb: Continental crust magmatic mixtures, Cabanis

Continental crust magmatic mixtures

<data plot> Reference [wiki]: Cabanis, B. Le diagramme La/10-Y/15-Nb/8 : unoutil pour la discrimination des series volcaniques et la mise en evidence des processus de melande et/ou de contamination crustale, 1989 [ref] Data: bulk Input: La_gpt, Y_gpt, Nb_gpt </data>

require('ggtern');

classifydiag.crust_pearce <- function (wtp_data, ...)
{  
ds <- selectNames(wtp_data,c('Name','La_gpt', 'Y_gpt', 'Nb_gpt'))

pts <- data.frame(Name=ds$Name,
                  X = ds$La_gpt/10,
                  Y = ds$Y_gpt/15,
                  Z = ds$Nb_gpt/8);

labels <- c("La_gpt/10", "Y_gpt/15", "Nb_gpt/8");
 
lines <- data.frame(X1 = c(0.00, 0.32, 0.43, 0.27, 0.12, 0.16, 0.22, 0.29),
                    Y1 = c(1.00, 0.68, 0.57, 0.56, 0.80, 0.68, 0.50, 0.28),
                    Z1 = c(0.00, 0.00, 0.00, 0.17, 0.08, 0.17, 0.38, 0.43),
                    X2 = c(0.62, 0.24, 0.34, 0.19, 0.38, 0.00, 0.00, 0.00),
                    Y2 = c(0.00, 0.62, 0.47, 0.60, 0.00, 0.76, 0.60, 0.42),
                    Z2 = c(0.38, 0.38, 0.20, 0.21, 0.62, 0.24, 0.40, 0.58));
 
marks  <- data.frame(Label=c("VAT", "NMORB", "*", "EMORB", "Calc-alkali", "Cont.", "Alkaline", "Intercont. rifts"),
                     X=c(0.24, 0.07, 0.20, 0.11, 0.70, 0.45, 0.17, 0.18),
                     Y=c(0.70, 0.80, 0.63, 0.45, 0.10, 0.10, 0.20, 0.10),
                     Z=c(0.07, 0.13, 0.17, 0.43, 0.20, 0.45, 0.63, 0.72));
ret <- c();
ret$type <- "ternogram";
ret$points <- pts;
ret$labels  <- labels;
ret$lines <- lines;
ret$marks <- marks;
 
return(ret);
}