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

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

Подробнее

Публикация TPF
Пересчет grt_tpf_realc0
Результат X_Alm, X_Sps, X_Prp, X_Grs, X_Ouv, X_Fe_p2_Fe_p2_Mg

Использовать Проверить

  • заголовок
  • краткое описание
  • статья
  • ссылка
  • карточка
  • подробное описание
  • код
  • тест

Код

  1. source('utilities_tpf.r');
  2.  
  3. interpmin.grt_tpf_recalc0 <- function (wtp_data, ...) {
  4.  
  5. rdat <- tpf.recalc_minaral("GRT", wtp_data)
  6.  
  7. rdat <- selectNames(rdat, c('Fe_p2_pfu','Fe_p3_pfu','Mg_pfu','Ca_pfu','Mn_pfu','Cr_pfu'));
  8.  
  9. rdat$X_Fe_p2_Fe_p2_Mg <- rdat$Fe_p2_pfu / (rdat$Fe_p2_pfu + rdat$Mg_pfu);
  10. so <- rdat$Fe_p2_pfu + rdat$Mg_pfu + rdat$Ca_pfu + rdat$Mn_pfu;
  11.  
  12. rdat$X_Alm <- rdat$Fe_p2_pfu / so;
  13. rdat$X_Sps <- rdat$Mn_pfu / so;
  14. rdat$X_Prp <- rdat$Mg_pfu / so;
  15. rdat$X_Grs <- (rdat$Ca_pfu - (rdat$Fe_p3_pfu+rdat$Cr_pfu)*1.5) / so;
  16. rdat$X_And <- rdat$Fe_p3_pfu*1.5 / so;
  17. rdat$X_Ouv <- rdat$Cr_pfu*1.5 / so;
  18.  
  19.  
  20. return(rdat);
  21. }
  22.  
  23. interpmin.grt_tpf_recalc0.test <- function (){
  24. warning("No test routine!");
  25. }
  • interpmin/grt_tpf_recalc0.txt
  • Последние изменения: 2018/11/12 15:07
  • (внешнее изменение)