Различия
Здесь показаны различия между двумя версиями данной страницы.
Следующая версия | Предыдущая версия Следующая версия Следующая версия справа и слева | ||
recalc:cpx_4cat_norm [2015/09/08 16:37] 83.149.255.10 создано |
recalc:cpx_4cat_norm [2015/09/09 07:45] admin |
||
---|---|---|---|
Строка 1: | Строка 1: | ||
- | Пересчет для граната на кислород с нормированием на катионы | + | ====== Пересчет для клинопироксена на кислород с нормированием на катионы ====== |
- | В расчет берутся: SiO2 TiO2 Al2O3 FeO MnO MgO CaO. Расчет валентного железа - по кристалохимии. | + | В расчет берутся: SiO2 TiO2 Al2O3 FeO MnO MgO CaO Na2O. Расчет валентного железа - по кристалохимии. |
По статье Ravna, 2001. | По статье Ravna, 2001. | ||
<code> | <code> | ||
- | source('recalc_common_oxygen.r') | + | source('recalc_common_cation.r') |
- | recalc.grt_12ox_norm <- function (wtp_data, ...) | + | recalc.cpx_4cat_norm <- function (wtp_data, ...) |
{ | { | ||
- | ox_group <-c('SiO2_wtp','TiO2_wtp', 'Al2O3_wtp', 'FeO_wtp', 'MnO_wtp', 'MgO_wtp', 'CaO_wtp'); | + | ox_group <-c('Name','SiO2_wtp','TiO2_wtp', 'Al2O3_wtp', 'FeO_wtp', 'MnO_wtp', 'MgO_wtp', 'CaO_wtp', 'Na2O_wtp'); |
- | wtp_group <- wtp_data[,c('Name','SiO2_wtp','TiO2_wtp', 'Al2O3_wtp', 'FeO_wtp', 'MnO_wtp', 'MgO_wtp', 'CaO_wtp')]; | + | wtp_group <- wtp_data[,ox_group]; |
- | ox_normed <- recalc.common_oxygen (wtp_group, 12); | + | cat_normed <- recalc.common_cation (wtp_group, 4); |
- | ox_normed$unnormed_cations <- rowSums(ox_normed[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu')]); | + | pfu_group <-c('Si_pfu','Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu', 'Na_pfu'); |
- | normed_cations <- ox_normed[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu')] * (5/ox_normed$unnormed_cations); | + | anion_summ <- rowSums(rowApply('*',cat_normed[,pfu_group],c(2,2,1.5,1,1,1,1,0.5))); |
- | normed_cations$Si_pfu = 3.0; | + | cat_normed$Fe_p2_pfu <- (6-anion_summ)*2; |
+ | cat_normed$Fe_p3_pfu <- cat_normed$Fe_pfu - cat_normed$Fe_p2_pfu; | ||
- | normed_cations$Fe_p3_pfu <- 2-rowSums(normed_cations[,c('Ti_pfu', 'Al_pfu')]); | + | cat_normed$Al_c4_pfu <- 2-cat_normed$Si_pfu; |
- | normed_cations$Fe_p2_pfu <- 3-rowSums(normed_cations[,c('Mn_pfu', 'Mg_pfu', 'Ca_pfu')]); | + | cat_normed$Al_c6_pfu <- cat_normed$Al_pfu-cat_normed$Al_c4_pfu; |
- | return (normed_cations); | + | return (cat_normed); |
} | } | ||
</code> | </code> |