Cpx: Putirka nomenative procedure, 2008

Putirka nomenative procedure.

Details

Reference K. D. Putirka Thermometers and Barometers for Volcanic Systems, 2008 [ref]
Recalc px_putirka
Results X_Jd, X_CaTs, X_CaTi, X_CrCaTs, X_DiHd, X_EnFs, X_En, X_Di

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

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

Код

source('recalc_px_putirka.r');

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

p <- recalc.px_putirka(wtp_data);
comps <- data.frame(Name = p$Name);

comps$X_Jd <- ifelse(p$Na_pfu<p$Al_c6_pfu, p$Na_pfu, p$Al_c6_pfu);
comps$X_CaTs <- ifelse(p$Al_c6_pfu > p$Na_pfu, p$Al_c6_pfu - p$Na_pfu, 0);
comps$X_CaTi <- ifelse(p$Al_c4_pfu > comps$X_CaTs, (p$Al_c4_pfu-comps$X_CaTs) / 2, 0);
comps$X_CrCaTs <- p$Cr_pfu / 2;
comps$X_DiHd <- p$Ca_pfu - comps$X_CaTs - comps$X_CaTi - comps$X_CrCaTs;
comps$X_EnFs <- (p$Fe_pfu + p$Mg_pfu - comps$X_DiHd) / 2;
comps$X_En <- comps$X_EnFs * (p$Mg_pfu / (p$Mg_pfu + p$Mn_pfu + p$Fe_pfu));
comps$X_Di <- comps$X_DiHd * (p$Mg_pfu / (p$Mg_pfu + p$Mn_pfu + p$Fe_pfu));

return(comps);
}

interpmin.cpx_putirka.test <- function (){
warning("No test routine!");
}