Q-J diagram Morimoto and Kitamura, 1983

Classifiaction of pyroxene in Q-J system.

require("ggplot2");
source("recalc_cpx_sturm.r");

classifydiag.cpx_morimoto_kitamura <- function (wtp_data, annotate=T, ...) {  

ds <- selectNames(recalc.cpx_sturm(wtp_data),
                  c("Name", "Ca_pfu", "Mg_pfu", "Fe_M2_pfu", "Na_pfu"));

pts <- data.frame(Name=ds$Name,
                  X=2*ds$Na_pfu,
                  Y=ds$Ca_pfu + ds$Fe_M2_pfu + ds$Mg_pfu);

marks  <- data.frame(Label=c("Quad", "Ca-Na", "Na", "Other",
                             "Jo,Ka,\nPe,Es", "En,Fs,Df,\nHd,Wo"),
                     X=c(0.1, 0.9, 1.5, 0.3, 0.0, 0.0),
                     Y=c(1.5, 0.9, 0.2, 0.5, 1.0, 2.0));

limits <- c(0,2,0,2);

lines <- data.frame(X1 = c(0.0, 0.0, 2.0, 1.5),
                    Y1 = c(0.0, 0.0, 0.0, 0.0),
                    X2 = c(2.0, 0.5, 0.0, 0.0),
                    Y2 = c(0.5, 2.0, 2.0, 1.5));

labels <- c("J", "Q");

ret <- c();

ret$points <- pts;
ret$limits <- limits;
ret$marks <- marks;
ret$lines <- lines;
ret$labels  <- labels;

return(ret);
}