Различия

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

--- recalc:cpx_4cat_norm [2015/09/08 16:37]
83.149.255.10 создано
+++ recalc:cpx_4cat_norm [2015/09/24 13:57]
admin
@@ Строка 1: / Строка 1: @@
-Пересчет для граната на кислород с нормированием на катионы
+====== Cpx: Пересчет для клинопироксена на кислород с нормированием на катионы ======
-В расчет берутся: SiO2 TiO2 Al2O3 FeO MnO MgO CaO. Расчет валентного железа - по кристалохимии.
+В расчет берутся: SiO2, TiO2, Al2O3, FeO, MnO, MgO, CaO, Na2O. Рассчитывается также содержание валентных форм железа и алюминия в IV и VI координациях.
-По статье Ravna, 2001.
+=== Подробнее ===
-<code>
+^ Информация о публикации | По статье Krogh Ravna, The garnet–clinopyroxene Fe2+–Mg geothermometer: an updated calibration, J. metamorphic Geol., 2000, 18, 211–219 |
-source('recalc_common_oxygen.r')
+^ Ссылка | {{:sources:ravna2000.pdf}} |
+^ Учитываемые компоненты | SiO2, TiO2, Al2O3, FeO, MnO, MgO, CaO, Na2O |
+^ Рассчитываемые компоненты | Fe_p2_pfu, Fe_p3_pfu, Al_c4_pfu, Al_c6_pfu |
-recalc.grt_12ox_norm <- function (wtp_data, ...)
+<Code linenums lang-r>
+source('recalc_common_cation.r')
+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 <- selectNames(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_p3_pfu <- (6-anion_summ)*2;
+cat_normed$Fe_p2_pfu <- cat_normed$Fe_pfu - cat_normed$Fe_p3_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>