Cpx: Nakamura, 2009

В расчет берутся: SiO2, TiO2, Al2O3, Cr2O3, FeO, MnO, MgO, CaO, Na2O, K2O.

Использовать Проверить

source('recalc_common_cation.r')

recalc.cpx_4cat_norm_jd <- 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_normed <- recalc.common_cation (wtp_group, 4);

cat_sx2 <- 2*(2 - cat_normed$Si_pfu - cat_normed$Ti_pfu);

cats <- wtp_group$SiO2_wtp/60.09 +
        wtp_group$TiO2_wtp/79.88 +
        wtp_group$Al2O3_wtp/101.96*2 +
        wtp_group$Cr2O3_wtp/152*2 + 
        wtp_group$FeO_wtp/71.85 + 
        wtp_group$MnO_wtp/70.94 + 
        wtp_group$MgO_wtp/40.31 + 
        wtp_group$CaO_wtp/56.08 + 
        wtp_group$Na2O_wtp/61.98*2 + 
        wtp_group$K2O_wtp/94.2*2;

cat_normed$Al_pfu <- wtp_group$Al2O3_wtp/101.96*2*4/cats;

jd <- cat_normed$Al_pfu-cat_sx2;

cat_normed$Fe_p3_pfu <- ifelse(cat_normed$Na_pfu+cat_normed$K_pfu-jd-cat_normed$Cr_pfu > 0, cat_normed$Na_pfu+cat_normed$K_pfu-jd-cat_normed$Cr_pfu, 0);

cat_normed$Fe_p2_pfu <- cat_normed$Fe_pfu - cat_normed$Fe_p3_pfu;

return (cat_normed);
}