Ol-Spl: Fabries, 1979 <B>

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

source("utilities_tpf.r");

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

    X <- Spl$Cr_pfu / (Spl$Al_c4_pfu + Spl$Al_c6_pfu + Spl$Cr_pfu + Spl$Fe_p3_pfu);
    Y <- Spl$Fe_p3_pfu / (Spl$Al_c4_pfu + Spl$Al_c6_pfu + Spl$Cr_pfu + Spl$Fe_p3_pfu);
    Z <- Spl$Mg_pfu / Spl$Fe_p2_pfu;
    A <- Ol$Mg_pfu / (Ol$Fe_p2_pfu + Ol$Fe_p3_pfu);
    T <- (4250 * X + 1343) / (log(A / Z) - 1.8 * Y + 1.825 * X + 0.571) - 273.15;
    warning(" FOR EQ.(5) LN(K)=1.8");
    return(data.frame(Name=paste0(Spl$Name,"-", Ol$Name),T_calc=T));
}

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

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