Различия

Здесь показаны различия между двумя версиями данной страницы.

--- recalc:amph_nakamura [2015/09/23 13:38]
admin
+++ recalc:amph_nakamura [2018/11/12 15:08] (текущий)
@@ Строка 1: / Строка 1: @@
-====== Amph: Пересчет для амфибола - кислородный метод ======
+====== Amp: Nakamura, 2009 ======
 Метод пересчета для граната на 23 атома кислорода. В расчет берутся оксиды: SiO2, TiO2, Al2O3, Cr2O3, FeO, MnO, MgO, CaO, Na2O, K2O.
-=== Подробнее ===
+=== Details ===
-^ Информация о публикации | |
+^ Reference | <ref:nakamura2009> |
-^ Ссылка |  |
+^ Input | SiO2_wtp, TiO2_wtp, Al2O3_wtp, Cr2O3_wtp, FeO_wtp, MnO_wtp, MgO_wtp, CaO_wtp, Na2O_wtp, K2O_wtp |
-^ Учитываемые компоненты | SiO2, TiO2, Al2O3, Cr2O3, FeO, MnO, MgO, CaO, Na2O, K2O |
+^ Output | pfu Si_T1_pfu, Ti_M2_pfu, Al_M2_pfu, Al_T1_pfu, F_pfu |
-^ Рассчитываемые компоненты | Si_T1_pfu, Ti_M2_pfu, F_pfu |
+<use:recalc.amph_nakamura>
+<test:recalc.amph_nakamura>
-<code>
+<WRAP center round todo 60%>
+  * заголовок
+  * краткое описание
+  * <del>статья</del>
+  * <del>ссылка</del>
+  * <del>карточка</del>
+  * подробное описание
+  * <del>код</del>
+  * тест
+</WRAP>
+=== Source ===
+<Code linenums lang-r>
 recalc.amph_nakamura <- function (wtp_data, ...)
 {
@@ Строка 26: / Строка 39: @@
             wtp_group$MnO/70.94 +
             wtp_group$MgO/40.31;
-pfu_group <- selectNames(wtp_data,c('Name'));
+pfu_group <- data.frame(Name = wtp_data$Name)
+pfu_group$Si_pfu <- wtp_group$SiO2_wtp/60.09 * 13 / cat_sub_NaCa;
+pfu_group$Si_T1_pfu <- pfu_group$Si_pfu - 4;
+pfu_group$Ti_pfu <- wtp_group$TiO2_wtp/79.88 * 13 /  cat_sub_NaCa;
+pfu_group$Al_pfu <- wtp_group$Al2O3_wtp/101.96 * 2 *13 /  cat_sub_NaCa;
+pfu_group$Al_T1_pfu <- 8 - pfu_group$Si_pfu;
+pfu_group$Al_M2_pfu <- pfu_group$Si_pfu + pfu_group$Al_pfu -8;
+pfu_group$Cr_pfu <- wtp_group$Cr2O3/152 * 2 * 13/cat_sub_NaCa;
+pfu_group$Fe_pfu <- wtp_group$FeO_wtp/71.85 * 13 /  cat_sub_NaCa;
+pfu_group$Mn_pfu <- wtp_group$MnO_wtp/70.94 *13 /  cat_sub_NaCa;
+pfu_group$Mg_pfu <- wtp_group$MgO_wtp/40.31 *13 /  cat_sub_NaCa;
+pfu_group$Ca_pfu <- wtp_group$CaO_wtp/56.08 *13 /  cat_sub_NaCa;
+pfu_group$Na_pfu <- wtp_group$Na2O_wtp/61.98 *2 *13 /  cat_sub_NaCa;
+pfu_group$K_pfu <- wtp_group$K2O_wtp/94.2 *2 *13 /  cat_sub_NaCa;
+pfu_group$F_pfu <- (pfu_group$Si_pfu*2 +
+                    pfu_group$Ti_pfu*2 +
+                    pfu_group$Al_pfu*1.5 +
+                    pfu_group$Cr_pfu*1.5 +
+                    pfu_group$Fe_pfu*1.5 +
+                    pfu_group$Mn_pfu +
+                    pfu_group$Mg_pfu +
+                    pfu_group$Ca_pfu +
+                    pfu_group$Na_pfu*0.5 +
+                    pfu_group$K_pfu*0.5-23)/(0.5*pfu_group$Fe_pfu);
 return (pfu_group);
 }
-</code>
+</Code>