Grt: Ravna, 2000

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

Details

Reference Krogh Ravna The garnet-clinopyroxene Fe2+-Mg geothermometer: an updated calibration, 2000 [ref]
Input SiO2_wtp, Al2O3_wtp, TiO2_wtp, FeO_wtp, MnO_wtp, MgO_wtp, CaO_wtp
Results Fe_p3_pfu, Fe_p2_pfu

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

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

Source

  1. source('recalc_common_oxygen.r')
  2.  
  3. recalc.grt_12ox_norm <- function (wtp_data, ...)
  4. {
  5. ox_group <-c('Name', 'SiO2_wtp','TiO2_wtp', 'Al2O3_wtp', 'FeO_wtp', 'MnO_wtp', 'MgO_wtp', 'CaO_wtp');
  6.  
  7. wtp_group <- selectNames(wtp_data,ox_group);
  8.  
  9. ox_normed <- recalc.common_oxygen (wtp_group, 12);
  10.  
  11. ox_normed$unnormed_cations <- rowSums(ox_normed[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu')]);
  12.  
  13. normed_cations <- ox_normed[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu')] * (5/ox_normed$unnormed_cations);
  14.  
  15. normed_cations$Si_pfu = 3.0;
  16.  
  17. normed_cations$Fe_p3_pfu <- 2-rowSums(normed_cations[,c('Ti_pfu', 'Al_pfu')]);
  18. normed_cations$Fe_p2_pfu <- 3-rowSums(normed_cations[,c('Mn_pfu', 'Mg_pfu', 'Ca_pfu')]);
  19.  
  20. ret <- data.frame(Name=wtp_data$Name,Si_pfu=ox_normed$Si_pfu);
  21.  
  22. return (cbind(ret,normed_cations[,c('Ti_pfu', 'Al_pfu', 'Fe_pfu', 'Mn_pfu', 'Mg_pfu', 'Ca_pfu','Fe_p2_pfu','Fe_p3_pfu')]));
  23.  
  24. }
  • recalc/grt_12ox_norm.txt
  • Последние изменения: 2020/02/19 14:29
  • — f0ma