Basalts Geodinamic Settings Meschede

Basalts Geodinamic Settings Meschede.

require('ggtern');

classifydiag.basalts_meschede <- function (wtp_data, annotate=T, labels=F, ...)
{  
ds <- selectNames(wtp_data,c('Name','Nb_gpt','Zr_gpt','Y_gpt'));

pts <- data.frame(Name=ds$Name,
                  X = ds$Zr_gpt/4,
                  Y = ds$Nb_gpt*10,
                  Z = ds$Y_gpt);

labels <- c("Zr_gpt/4", "2*Nb_gpt", "Y_gpt");
 
ret <- c();
ret$type <- "ternogram";
ret$points <- pts;
ret$labels  <- labels;
 
return(ret);
}