====== Grt: Интерпретация состава граната по TPF, 1991 ======

Рассчитывается Fe2/(Fe2+Mg) и 6 компонентов граната: Almandine, Spessartine, Pyrope, Grossular, Andradite, Ouvarovite.

=== Подробнее ===

^ Публикация | [[:about_tpf|TPF]] |
^ Пересчет | [[recalc:grt_tpf_realc0]] |
^ Результат | X_Alm, X_Sps, X_Prp, X_Grs, X_Ouv, X_Fe_p2_Fe_p2_Mg |

<use:interpmin.grt_tpf_recalc0>
<test:interpmin.grt_tpf_recalc0>

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

==== Код ====
<Code linenums lang-r>
source('utilities_tpf.r');

interpmin.grt_tpf_recalc0 <- function (wtp_data, ...) {

rdat <- tpf.recalc_minaral("GRT", wtp_data)

rdat <- selectNames(rdat, c('Fe_p2_pfu','Fe_p3_pfu','Mg_pfu','Ca_pfu','Mn_pfu','Cr_pfu'));

rdat$X_Fe_p2_Fe_p2_Mg <- rdat$Fe_p2_pfu / (rdat$Fe_p2_pfu + rdat$Mg_pfu);
so <- rdat$Fe_p2_pfu + rdat$Mg_pfu + rdat$Ca_pfu + rdat$Mn_pfu;

rdat$X_Alm <- rdat$Fe_p2_pfu / so;
rdat$X_Sps <- rdat$Mn_pfu / so;
rdat$X_Prp <- rdat$Mg_pfu / so;
rdat$X_Grs <- (rdat$Ca_pfu - (rdat$Fe_p3_pfu+rdat$Cr_pfu)*1.5) / so;
rdat$X_And <- rdat$Fe_p3_pfu*1.5 / so;
rdat$X_Ouv <- rdat$Cr_pfu*1.5 / so;


return(rdat);
}

interpmin.grt_tpf_recalc0.test <- function (){
warning("No test routine!");
}
</Code>

{{tag> order71}}