Au-Sph-Py-Arg: Shikazono, 1985 <A>

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

source("utilities_tpf.r");

#main code
sensors.t.au_sph_py_arg_shikazono_1985_a <- function(Au, Sph, Py, Arg, ...){
    Au <- tpf.recalc_minaral("AU",Au);
    Sph <- tpf.recalc_minaral("SPH",Sph);

    N <- AG / (Au$Fe_pfu + CO + NI + CU + ZN + MO + AG + CD + SN + AU + HG + PB + AS + SB + S);
    X <- Sph$Fe_pfu / (Sph$Fe_pfu + CO + NI + CU + ZN + MO + AG + CD + SN + AU + HG + PB + AS + SB + S);
    T <- (28765 + 22600 * (1 - N)^2 - 6400 * (1 - N)^3) /      (49.008 - 9.152 * log(X) + 18.2961 * log(N) + 5.5 * (1 - N)^2) - 273.15;
    warning(" EQUATION (6)");
    return(data.frame(Name=paste0(Au$Name,"-", Sph$Name),T_calc=T));
}

sensors.t.au_sph_py_arg_shikazono_1985_a.args <- function(){
    ars <- list();
    ars$minerals <- c("Au", "Sph", "Py", "Arg");
    ars$ranges <- c();
    return(ars);
}

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