====== Px: Ca-Mg-Fe and Na pyroxenes with accepted names, Morimoto et al., 1988 ======

Ca-Mg-Fe and Na pyroxenes with accepted names. Quad represents the Ca-Mg-Fe pyroxene area [[classifydiag:cpx_ca_mg_fe]].

=== Details ===

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

<use:classifydiag.px_na_ca_mg_fe>
<test:classifydiag.px_na_ca_mg_fe>

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

=== Source ===

<Code linenums lang-r>
source('recalc_cpx_4cat_norm.r');

classifydiag.px_na_ca_mg_fe <- function (wtp_data, annotate=T, labels=F, ...)
{  
ds <- selectNames(recalc.cpx_4cat_norm(wtp_data),
                  c('Name','Ca_pfu','Mg_pfu','Fe_pfu','Fe_p2_pfu','Fe_p3_pfu','Na_pfu'));

//marks <- data.frame(Label=c('Diopside','Hedenbergite','Augite','Pigeonite','Clinoenstotite','Clinoferrosilite'),
//                     X=c(0.48,0.12,0.33,0.44,0.78,0.18),
//                     Y=c(0.55,0.55,0.33,0.12,0.03,0.03),
//                     Z=c(0.12,0.48,0.33,0.44,0.18,0.78));

pts <- data.frame (Name = ds$Name,
                   X = ds$Al_pfu,
                   Y = 2*ds$Na_pfu,
                   Z = ds$Fe_p3_pfu);

labels <- c("Jd", "Q", "Ae");

ret <- c();
 
ret$type <- "ternogram";
ret$points <- pts;
ret$marks <- marks;
ret$labels  <- labels;

return (ret);
}
</Code>

{{tag>classifydiag mineral}}