Grt-Hbl: LavrentEva, Perchuk, 1989 <A>

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

source("utilities_tpf.r");

sensors.t.grt_hbl_lavrenteva_perchuk_1989_a <- function(Grt, Hbl, ...){
    Grt <- tpf.recalc_minaral("GRT2",Grt);
    Hbl <- tpf.recalc_minaral("HBL",Hbl);

    C <- Grt$Ca_pfu;
    M <- Grt$Mg_pfu / (Grt$Fe_p2_pfu + Grt$Fe_p3_pfu + Grt$Mg_pfu + Grt$Mn_pfu);
    S <- Hbl$Si_pfu + Hbl$Ti_pfu + Hbl$Al_c4_pfu + Hbl$Al_c6_pfu + Hbl$Cr_pfu + Hbl$Fe_p2_pfu + Hbl$Fe_p3_pfu + Hbl$Mn_pfu + Hbl$Mg_pfu + Hbl$Ca_pfu + Hbl$Na_pfu;
    N <- Hbl$Na_pfu;
    A <- Hbl$Al_c4_pfu + Hbl$Al_c6_pfu - Hbl$Na_pfu;
    CH <- Hbl$Ca_pfu;
    MA <- Hbl$Mg_pfu / (Hbl$Fe_p2_pfu + Hbl$Fe_p3_pfu + Hbl$Mg_pfu + Hbl$Mn_pfu);
    K <- 15 / S;
    N <- N * K;
    A <- A * K;
    CH <- CH * K;
    LNK <- log(MA / (1 - MA) * (1 - M) / M);
    T <- (3366 - 1444 * N - 323 * A + 188 * CH + C * (728 * C - 148)) / (LNK + 2.417 - 1.393 * N -      0.266 * A + 0.051 * CH) - 273;
    warning(" CA,NA,AL CONTENT IN HBL IS TAKEN INTO ACCOUNT");
    return(data.frame(Name=paste0(Hbl$Name,"-", Grt$Name),T_calc=T));
}

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

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