Grt-L: Ellis, 1986

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

source("utilities_tpf.r");

sensors.t.grt_l_ellis_1986 <- function(Grt, L, P, ...){
    L <- tpf.recalc_minaral("L",L);
    Grt <- tpf.recalc_minaral("GRT",Grt);

    X <- (Grt$Fe_p2_pfu + Grt$Fe_p3_pfu) / Grt$Mg_pfu;
    Y <- (L$Fe_p2_pfu + L$Fe_p3_pfu) / L$Mg_pfu;
    K <- log(Y / X);
    T <- (4975 + 0.0166 * P * 1000) / (4.27 - K) - 273;
    return(data.frame(Name=paste0(L$Name,"-", Grt$Name),P=P,T_calc=T));
}

sensors.t.grt_l_ellis_1986.args <- function(){
    ars <- list();
    ars$minerals <- c("Grt", "L");
    ars$ranges <- c("P");
    return(ars);
}

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