Amph: Пересчет для амфибола - кислородный метод

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

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 <- selectNames(wtp_data,c('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/70.94 *13 /  cat_sub_NaCa;
pfu_group$Mg_pfu/40.31 *13 /  cat_sub_NaCa;
pfu_group$Ca_pfu/56.08 *13 /  cat_sub_NaCa;
pfu_group$Na_pfu/61.98 *2 *13 /  cat_sub_NaCa;
pfu_group$K_pfu/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);

return (pfu_group);
}