====== Cpx: Ravna, 2000 ======

В расчет берутся: SiO2, TiO2, Al2O3, FeO, MnO, MgO, CaO, Na2O. Рассчитывается также содержание валентных форм железа и алюминия в IV и VI координациях.

=== Details ===

^ Reference | <ref:ravna2000> |
^ Input | SiO2_wtp, TiO2_wtp, Al2O3_wtp, FeO_wtp, MnO_wtp, MgO_wtp, CaO_wtp, Na2O_wtp |
^ Result | pfu, Fe_p2_pfu, Fe_p3_pfu, Al_c4_pfu, Al_c6_pfu |

<use:recalc.cpx_4cat_norm>
<test:recalc.cpx_4cat_norm>

<WRAP center round todo 60%>
  * заголовок
  * <del>краткое описание</del>
  * <del>статья</del>
  * <del>ссылка</del>
  * <del>карточка</del>
  * подробное описание
  * <del>код</del>
  * тест
</WRAP>

=== Source ===

<Code linenums lang-r>
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 <- selectNames(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');

anion_summ <- rowSums(rowApply('*',cat_normed[,pfu_group],c(2,2,1.5,1,1,1,1,0.5)));

cat_normed$Fe_p3_pfu <- (6-anion_summ)*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);
}
</Code>