Ol-Spl: Hill, Sack, 1987

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

source("utilities_tpf.r");

sensors.t.ol_spl_hill_sack_1987 <- function(Ol, Spl, ...){
    Spl <- tpf.recalc_minaral("SPL",Spl);
    Ol <- tpf.recalc_minaral("OL",Ol);

    Y <- (Ol$Fe_p2_pfu + Ol$Fe_p3_pfu) / (Ol$Fe_p2_pfu + Ol$Fe_p3_pfu + Ol$Mg_pfu);
    Y1 <- Ol$Mg_pfu / (Ol$Fe_p2_pfu + Ol$Fe_p3_pfu);
    X1 <- Spl$Mg_pfu / Spl$Fe_p2_pfu;
    X3 <- Spl$Cr_pfu / 2;
    X4 <- Spl$Ti_pfu;
    X5 <- Spl$Fe_p3_pfu / 2;
    A <- Spl$Al_pfu;
    if (A>0.20) {
    warning("CANNOT BE USED FOR AL IN SPL >0.2");
    }
    if (A<0.20) {
    K <- 0.0019872 * log(Y1 / X1);
    T <- (1.33795 + 4 * 80 * X3 + 3.5561 * X4 + 6.39 * X5 + 0.84 * (1 - 2 * Y)) /        (K + 0.0009539 - 0.0015698 * X4) - 273;
    }
    return(data.frame(Name=paste0(Spl$Name,"-", Ol$Name),T_calc=T));
}

sensors.t.ol_spl_hill_sack_1987.args <- function(){
    ars <- list();
    ars$minerals <- c("Ol", "Spl");
    ars$ranges <- c();
    return(ars);
}

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