Пересчет для клинопироксена на кислород с нормированием на катионы

В расчет берутся: SiO2 TiO2 Al2O3 FeO MnO MgO CaO Na2O. Расчет валентного железа - по кристалохимии.

По статье Ravna, 2001.

source('recalc_common_cation.r')

recalc.cpx_4cat_norm <- function (wtp_data, ...)
{

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

wtp_group <- wtp_data[,ox_group];

cat_normed <- recalc.common_cation (wtp_group, 4);

pfu_group <-c('Si_pfu','Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu', 'Na_pfu');

ox_calc <- rowSums(rowApply('*',cat_normed[,pfu_group],c(2,2,1.5,1,1,1,1,0.5)));
cat_normed$Fe_p3_pfu <- (6-ox_calc)*2;
cat_normed$Fe_p2_pfu <- cat_normed$Fe_pfu - cat_normed$Fe_p3_pfu;

cat_normed$Al_c4_pfu <- 2-cat_normed$Si_pfu;
cat_normed$Al_c6_pfu <- cat_normed$Al_pfu-cat_normed$Al_c4_pfu;

return (cat_normed);
}