Putirka nomenative procedure.
| 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!");
}