====== Amp: Nakamura, 2009 ======

Метод пересчета для граната на 23 атома кислорода. В расчет берутся оксиды: SiO2, TiO2, Al2O3, Cr2O3, FeO, MnO, MgO, CaO, Na2O, K2O.

=== Details ===

^ Reference | <ref:nakamura2009> |
^ Input | SiO2_wtp, TiO2_wtp, Al2O3_wtp, Cr2O3_wtp, FeO_wtp, MnO_wtp, MgO_wtp, CaO_wtp, Na2O_wtp, K2O_wtp |
^ Output | pfu Si_T1_pfu, Ti_M2_pfu, Al_M2_pfu, Al_T1_pfu, F_pfu |
<use:recalc.amph_nakamura>
<test:recalc.amph_nakamura>

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

=== Source ===

<Code linenums lang-r>
recalc.amph_nakamura <- function (wtp_data, ...)
{

ox_group <-c('Name','SiO2_wtp','TiO2_wtp', 'Al2O3_wtp', 'Cr2O3_wtp', 'FeO_wtp', 'MnO_wtp', 'MgO_wtp', 'CaO_wtp', 'Na2O_wtp', 'K2O_wtp');

wtp_group <- selectNames(wtp_data,ox_group);

cat_sub_NaCa <- wtp_group$SiO2_wtp/60.09 +
            wtp_group$TiO2_wtp/79.88 +
            wtp_group$Al2O3_wtp/101.96*2 +
            wtp_group$Cr2O3/152*2 +
            wtp_group$FeO/71.85 +
            wtp_group$MnO/70.94 +
            wtp_group$MgO/40.31;


pfu_group <- data.frame(Name = wtp_data$Name)

pfu_group$Si_pfu <- wtp_group$SiO2_wtp/60.09 * 13 / cat_sub_NaCa;
pfu_group$Si_T1_pfu <- pfu_group$Si_pfu - 4;
pfu_group$Ti_pfu <- wtp_group$TiO2_wtp/79.88 * 13 /  cat_sub_NaCa;
pfu_group$Al_pfu <- wtp_group$Al2O3_wtp/101.96 * 2 *13 /  cat_sub_NaCa;
pfu_group$Al_T1_pfu <- 8 - pfu_group$Si_pfu;
pfu_group$Al_M2_pfu <- pfu_group$Si_pfu + pfu_group$Al_pfu -8;
pfu_group$Cr_pfu <- wtp_group$Cr2O3/152 * 2 * 13/cat_sub_NaCa;
pfu_group$Fe_pfu <- wtp_group$FeO_wtp/71.85 * 13 /  cat_sub_NaCa;
pfu_group$Mn_pfu <- wtp_group$MnO_wtp/70.94 *13 /  cat_sub_NaCa;
pfu_group$Mg_pfu <- wtp_group$MgO_wtp/40.31 *13 /  cat_sub_NaCa;
pfu_group$Ca_pfu <- wtp_group$CaO_wtp/56.08 *13 /  cat_sub_NaCa;
pfu_group$Na_pfu <- wtp_group$Na2O_wtp/61.98 *2 *13 /  cat_sub_NaCa;
pfu_group$K_pfu <- wtp_group$K2O_wtp/94.2 *2 *13 /  cat_sub_NaCa;

pfu_group$F_pfu <- (pfu_group$Si_pfu*2 +
                    pfu_group$Ti_pfu*2 +
                    pfu_group$Al_pfu*1.5 +
                    pfu_group$Cr_pfu*1.5 + 
                    pfu_group$Fe_pfu*1.5 +
                    pfu_group$Mn_pfu + 
                    pfu_group$Mg_pfu +
                    pfu_group$Ca_pfu +
                    pfu_group$Na_pfu*0.5 +
                    pfu_group$K_pfu*0.5-23)/(0.5*pfu_group$Fe_pfu);

return (pfu_group);
}
</Code>