recalc:grt_12ox_norm [pin.ifz.ru]

Эта страница только для чтения. Вы можете посмотреть её исходный текст, но не можете его изменить. Сообщите администратору, если считаете, что это неправильно.

====== Grt: Ravna, 2000 ======

Расчет на 23 атомов кислорода: SiO2 TiO2 Al2O3 FeO MnO MgO CaO. Расчет валентного железа - по кристалохимии.

=== Details ===

^ Reference | <ref:ravna2000> |
^ Input | SiO2_wtp, Al2O3_wtp, TiO2_wtp, FeO_wtp, MnO_wtp, MgO_wtp, CaO_wtp |
^ Results | Fe_p3_pfu, Fe_p2_pfu |
<use:recalc.grt_12ox_norm>
<test:recalc.grt_12ox_norm>

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

=== Source ===

<Code linenums lang-r>
source('recalc_common_oxygen.r')

recalc.grt_12ox_norm <- function (wtp_data, ...)
{
ox_group <-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);

ox_normed$unnormed_cations <- rowSums(ox_normed[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu')]);

normed_cations <- ox_normed[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu')] * (5/ox_normed$unnormed_cations);

normed_cations$Si_pfu = 3.0;

normed_cations$Fe_p3_pfu <- 2-rowSums(normed_cations[,c('Ti_pfu', 'Al_pfu')]);
normed_cations$Fe_p2_pfu <- 3-rowSums(normed_cations[,c('Mn_pfu', 'Mg_pfu', 'Ca_pfu')]);

ret <- data.frame(Name=wtp_data$Name,Si_pfu=ox_normed$Si_pfu);

return (cbind(ret,normed_cations[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu','Fe_p2_pfu','Fe_p3_pfu')]));

}
</Code>