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, labels=F, ...)
{  
font_size=3;
line_height=0.7;
ds <- selectNames(recalc.cpx_sturm(wtp_data),c('Name','Ca_pfu','Mg_pfu','Fe_M2_pfu','Na_pfu'));

ds$J <- 2*ds$Na_pfu;
ds$Q <- ds$Ca_pfu + ds$Fe_M2_pfu + ds$Mg_pfu;

txts <- data.frame(Label=c('Quad','Ca-Na','Na','Other','Jo,Ka,\nPe,Es','En,Fs,Df,\nHd,Wo'),
                     J=c(0.1,0.9,1.5,0.3,0,0.0),
                     Q=c(1.5,0.9,0.2,0.5,1,2));

p <- ggplot(ds, aes(x=J,y=Q)) +
  geom_point() +
  theme_bw() +
  xlim(c(0,2)) +
  ylim(c(0,2)) + 
  geom_abline(intercept = 0, slope = 0.25) +
  geom_abline(intercept = 0, slope = 4) +
  geom_abline(intercept = 1.5, slope = -1) +
  geom_abline(intercept = 2.0, slope = -1) +
  xlab("J")+
  ylab("Q")+
  geom_text(data=txts, aes(label=Label), size=font_size, lineheight=line_height);
  
if (labels)
{
    p <- p + geom_text(aes(label=Name),hjust=0, vjust=0, size=font_size, lineheight=line_height);
}
  
  
return(p);
}