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;

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.2) +
  geom_abline(intercept = 0, slope = 5) +
  geom_abline(intercept = 1.5, slope = -1) +
  geom_abline(intercept = 2.0, slope = -1) +
  xlab("J")+
  ylab("Q");
  
if (labels)
{
    p <- p + geom_text(aes(label=Name),hjust=0, vjust=0, size=font_size, lineheight=line_height);
}
  
  
return(p);
}