====== Px: Q-J Diagram, Morimoto, 1988 ======

Classifiaction of pyroxene in Q-J system.

=== Details ===

| Reference | <ref:morimoto1988> |
| Data | mineral |
| Input | wtp |

<use:classifydiag.cpx_morimoto_kitamura>
<test:classifydiag.cpx_morimoto_kitamura>

<WRAP center round todo 60%>
  * <del>заголовок</del>
  * <del>краткое описание</del>
  * <del>статья</del>
  * <del>ссылка</del>
  * <del>карточка</del>
  * подробное описание
  * <del>код</del>
  * тест
</WRAP>

=== Код метода ===

<Code linenums lang-r>

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$type <- "diagram";
ret$points <- pts;
ret$limits <- limits;
ret$marks <- marks;
ret$lines <- lines;
ret$labels  <- labels;

return(ret);
}
</Code>

{{tag>classifydiag mineral}}