Mag-Ilm: Powell, Powell, 1977 <A>

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

source("utilities_tpf.r");

sensors.t.mag_ilm_powell_powell_1977_a <- function(Mag, Ilm, ...){
    Ilm <- tpf.recalc_minaral("ILM",Ilm);
    Mag <- tpf.recalc_minaral("MAG",Mag);

    XMA <- Mag$Ti_pfu;
    XMI <- Mag$Ti_pfu - 0.5 * (Mag$Mg_pfu + Mag$Mn_pfu + Mag$Ca_pfu);
    YMI <- Ilm$Ti_pfu - (Ilm$Mg_pfu + Ilm$Mn_pfu + Ilm$Ca_pfu);
    YMA <- Ilm$Ti_pfu - (Ilm$Mg_pfu + Ilm$Mn_pfu + Ilm$Ca_pfu);
    ZMI <- Mag$Fe_p3_pfu / 2 - (Mag$Mg_pfu + Mag$Mn_pfu + Mag$Ca_pfu);
    ZMA <- Mag$Fe_p3_pfu / 2;
    HMI <- 0.5 * (Ilm$Fe_p3_pfu - Ilm$Al_pfu - Ilm$Cr_pfu);
    HMA <- Ilm$Fe_p3_pfu / 2;
    KMA <- log(XMA * HMA / ZMI / YMI);
    KMI <- log(XMI * HMI / ZMA / YMA);
    TMA <- 8155 / (4.59 - KMA) - 273;
    TMI <- 8155 / (4.59 - KMI) - 273;
    T <- (TMI + TMA) / 2;
    return(data.frame(Name=paste0(Ilm$Name,"-", Mag$Name),T_calc=T));
}

sensors.t.mag_ilm_powell_powell_1977_a.args <- function(){
    ars <- list();
    ars$minerals <- c("Mag", "Ilm");
    ars$ranges <- c();
    return(ars);
}

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