====== Cpx: Nakamura, 2009 ======

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

=== Details ===

^ Reference | <ref:nakamura2009> |
^ Inputs | SiO2_wtp, TiO2_wtp, Al2O3_wtp, Cr2O3_wtp, FeO_wtp, MnO_wtp, MgO_wtp, CaO_wtp, Na2O_wtp, K2O_wtp |
^ Results | pfu, Fe_p2_pfu, Fe_p3_pfu |

<use:recalc.cpx_4cat_norm_jd>
<test:recalc.cpx_4cat_norm_jd>

<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_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);
}
</Code>