Это старая версия документа!
Cpx: Putirka nomenative procedure
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!");
- }