Pl-Hbl-Qtz: Holland, Blundy, 1994

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

source("utilities_tpf.r");

sensors.t.pl_hbl_qtz_holland_blundy_1994 <- function(Pl, Hbl, Qtz, P, ...){
    Pl <- tpf.recalc_minaral("PL",Pl);
    Hbl <- tpf.recalc_minaral("HBL",Hbl);

    CM <- (Si + Al + Ti + Hbl$Fe_p3_pfu + Hbl$Fe_p2_pfu + Hbl$Mg_pfu + Hbl$Mn_pfu - 13.0);
    XSI <- (Hbl$Si_pfu - 4) / 4;
    XALT <- (8 - Hbl$Si_pfu) / 4;
    XALM <- (Hbl$Al_pfu + Hbl$Si_pfu - 8) / 2;
    XK <- Hbl$K_pfu;
    XV <- 3 - Hbl$Ca_pfu - Hbl$Na_pfu - Hbl$K_pfu - (Si + Al + Ti + Hbl$Fe_p3_pfu + Hbl$Fe_p2_pfu + Hbl$Mg_pfu + Hbl$Mn_pfu - 13.0);
    XNAA <- Hbl$Ca_pfu + Hbl$Na_pfu + (Si + Al + Ti + Hbl$Fe_p3_pfu + Hbl$Fe_p2_pfu + Hbl$Mg_pfu + Hbl$Mn_pfu - 13.0) - 2;
    XNAM <- (2 - Hbl$Ca_pfu - (Si + Al + Ti + Hbl$Fe_p3_pfu + Hbl$Fe_p2_pfu + Hbl$Mg_pfu + Hbl$Mn_pfu - 13.0)) / 2;
    XCA <- Hbl$Ca_pfu / 2;
    XAB <- Pl$Na_pfu / (Pl$Ca_pfu + Pl$Na_pfu + Pl$K_pfu);
    if (XAB>0.5) {
    YAB <- 0.0;
    }
    if (XAB<0.5) {
    YAB <- 12.0 * (1 - XAB)^2 - 3.0;
    }
    if (XAB==0.5) {
    YAB <- 0.0;
    }
    KR <- (27 * XV * XSI * XAB) / (256 * XNAA * XALT);
    T <- ( - 76.95 + 0.79 * P + YAB + 39.4 * XNAA + 22.4 * XK + (41.5 - 2.89 * P) * XALM) /      ( - 0.0650 - 0.0083144 * log(KR)) - 287;
    warning(" Edenite-tremolite thermometer, NO CORRECTION FOR FE3 IS MADE ! ");
    return(data.frame(Name=paste0(Hbl$Name,"-", Pl$Name),P=P,T_calc=T));
}

sensors.t.pl_hbl_qtz_holland_blundy_1994.args <- function(){
    ars <- list();
    ars$minerals <- c("Pl", "Hbl", "Qtz");
    ars$ranges <- c("P");
    return(ars);
}

#sensors.t.pl_hbl_qtz_holland_blundy_1994.test <- function(){
#
#}