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

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

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

source('recalc_common_oxygen.r')

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

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

wtp_group <- wtp_data[,c('Name','SiO2_wtp','TiO2_wtp', 'Al2O3_wtp', 'FeO_wtp', 'MnO_wtp', 'MgO_wtp', 'CaO_wtp')];

ox_normed <- recalc.common_oxygen (wtp_group, 12);

ox_normed$unnormed_cations <- rowSums(ox_normed[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu')]);

normed_cations <- ox_normed[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu')] * (5/ox_normed$unnormed_cations);

normed_cations$Si_pfu = 3.0;

normed_cations$Fe_p3_pfu <- 2-rowSums(normed_cations[,c('Ti_pfu', 'Al_pfu')]);
normed_cations$Fe_p2_pfu <- 3-rowSums(normed_cations[,c('Mn_pfu', 'Mg_pfu', 'Ca_pfu')]);

return (normed_cations);
}