Hbl-Pl: Spear, 1981 <B>

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

source('utilities_tpf.r');

#main code
sensors.t.hbl_pl_spear_1981_b <- function(Hbl, Pl, P, ...){
    Pl <- tpf.recalc_minaral('PL',Pl);
    Hbl <- tpf.recalc_minaral('HBL',Hbl);

    X <- Pl$Ca_pfu/Pl$Na_pfu;
    Y <- Hbl$Ca_pfu;
    Z <- Hbl$Na_pfu;
    K <- Hbl$Mg_pfu+Hbl$Fe_p2_pfu+Hbl$Fe_p3_pfu+Hbl$Ti_pfu+Hbl$Cr_pfu+Hbl$Al_c4_pfu+Hbl$Al_c6_pfu+Hbl$Si_pfu;
    N <- 13/K;
    Y <- Y*N;
    Z <- Z*N;
    if (Y<2.0) {
    if (Y+Z>2) {
    N <- 2-Y;
    T <- (0.24*P*1000-47000)/(2*1.9872*log(X*N/Y)-44.27)-273;
    }
    if (Y+Z<2) {
    T <- (0.24*P*1000-47000)/(2*1.9872*log(X*Z/Y)-44.27)-273;
    }
    }
    if (Y>2.0) {
    warning(' CANNOT BE USED FOR CA IN HBL > 2')
    }
    return(data.frame(Name=paste0(Pl$Name,'-', Hbl$Name),P=P,T_calc=T));
}

sensors.t.hbl_pl_spear_1981_b.args <- function(){
    ars <- list();
    ars$minerals = c('Hbl', 'Pl');
    ars$ranges = c('P');
    return(ars);
}

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