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

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

source("utilities_tpf.r");

sensors.f.mag_ilm_powell_powell_1977_b <- 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);
    FMA <- (2 * log(ZMI) + 6 * log(YMI) - 6 * log(XMA) + 10.702 * KMA - 4.94) / 2.30258;
    FMI <- (2 * log(ZMA) + 6 * log(YMA) - 6 * log(XMI) + 10.702 * KMI - 4.94) / 2.30258;
    LGFO2 <- (FMA + FMI) / 2;
    return(data.frame(Name=paste0(Ilm$Name,"-", Mag$Name),LGFO2_calc=LGFO2));
}

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

#sensors.f.mag_ilm_powell_powell_1977_b.test <- function(){
#
#}