Это старая версия документа!


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

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

Код

  1. source('recalc_px_putirka.r');
  2.  
  3. interpmin.cpx_putirka <- function (wtp_data, ...) {
  4.  
  5. p <- recalc.px_putirka(wtp_data);
  6. comps <- data.frame(Name = p$Name);
  7.  
  8. comps$X_Jd <- ifelse(p$Na_pfu<p$Al_c6_pfu, p$Na_pfu, p$Al_c6_pfu);
  9. comps$X_CaTs <- ifelse(p$Al_c6_pfu > p$Na_pfu, p$Al_c6_pfu - p$Na_pfu, 0);
  10. comps$X_CaTi <- ifelse(p$Al_c4_pfu > comps$X_CaTs, (p$Al_c4_pfu-comps$X_CaTs) / 2, 0);
  11. comps$X_CrCaTs <- p$Cr_pfu / 2;
  12. comps$X_DiHd <- p$Ca_pfu - comps$X_CaTs - comps$X_CaTi - comps$X_CrCaTs;
  13. comps$X_EnFs <- (p$Fe_pfu + p$Mg_pfu - comps$X_DiHd) / 2;
  14. comps$X_En <- comps$X_EnFs * (p$Mg_pfu / (p$Mg_pfu + p$Mn_pfu + p$Fe_pfu));
  15. comps$X_Di <- comps$X_DiHd * (p$Mg_pfu / (p$Mg_pfu + p$Mn_pfu + p$Fe_pfu));
  16.  
  17. return(comps);
  18. }
  19.  
  20. interpmin.cpx_putirka.test <- function (){
  21. warning("No test routine!");
  22. }
  • interpmin/cpx_putirka.1458816305.txt.gz
  • Последние изменения: 2018/11/12 15:08
  • (внешнее изменение)