Ol-Spl: Roeder, 1979

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

source("utilities_tpf.r");

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

    X <- Spl$Cr_pfu / (Spl$Cr_pfu + Spl$Al_pfu + Spl$Fe_p3_pfu);
    Y <- Spl$Al_pfu / (Spl$Cr_pfu + Spl$Al_pfu + Spl$Fe_p3_pfu);
    Z <- Spl$Fe_p3_pfu / (Spl$Cr_pfu + Spl$Al_pfu + Spl$Fe_p3_pfu);
    N <- Spl$Fe_p2_pfu / Spl$Mg_pfu;
    M <- Ol$Mg_pfu / Ol$Fe_p2_pfu;
    T <- ((X * 3480 + Y * 1018 - Z * 1720 + 2400) / (X * 2.23 + Y * 2.56 -      Z * 3.08 - 1.47 + 1.987 * log(M * N))) - 273;
    return(data.frame(Name=paste0(Spl$Name,"-", Ol$Name),T_calc=T));
}

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

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