Hbl-Pl: Holland, Blundy, 1994

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

source("utilities_tpf.r");

sensors.t.hbl_pl_holland_blundy_1994 <- function(Hbl, Pl, 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);
    XAN <- Pl$Ca_pfu / (Pl$Ca_pfu + Pl$Na_pfu + Pl$K_pfu);
    if (XAB>0.5) {
    YABAN <- 3.0;
    }
    if (XAB<0.5) {
    YABAN <- 12.0 * (2 * XAB - 1) + 3.0;
    }
    if (XAB==0.5) {
    YABAN <- 3.0;
    }
    KR <- (27 * XNAM * XSI * XAN) / (64 * XCA * XALT * XAB);
    T <- (78.44 + YABAN - 33.6 * XNAM - (66.8 - 2.92 * P) * XALM + 78.5 * XALT + 9.4 * XNAA) /      (0.0721 - 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.hbl_pl_holland_blundy_1994.args <- function(){
    ars <- list();
    ars$minerals <- c("Hbl", "Pl");
    ars$ranges <- c("P");
    return(ars);
}

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