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 <- Au$Ag_pfu / (Au$Fe_pfu + Au$Co_pfu + Au$Ni_pfu + Au$Cu_pfu + Au$Zn_pfu + Au$Mo_pfu + Au$Ag_pfu + Au$Cd_pfu + Au$Au$S_pfun_pfu + Au$Au_pfu + Au$Hg_pfu + Au$Pb_pfu + Au$As_pfu + Au$Au$S_pfub_pfu + Au$S_pfu);
    X <- Sph$S_pfuph$Fe_pfu / (Sph$S_pfuph$Fe_pfu + Sph$S_pfuph$Co_pfu + Sph$S_pfuph$Ni_pfu + Sph$S_pfuph$Cu_pfu + Sph$S_pfuph$Zn_pfu + Sph$S_pfuph$Mo_pfu + Sph$S_pfuph$Ag_pfu + Sph$S_pfuph$Cd_pfu + Sph$S_pfuph$Sph$S_pfun_pfu + Sph$S_pfuph$Au_pfu + Sph$S_pfuph$Hg_pfu + Sph$S_pfuph$Pb_pfu + Sph$S_pfuph$As_pfu + Sph$S_pfuph$Sph$S_pfub_pfu + Sph$S_pfu);
    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(){
#}