Ol-Spl-Opx: Sen, 1985

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

source("utilities_tpf.r");

sensors.t.ol_spl_opx_sen_1985 <- function(Ol, Spl, Opx, ...){
    Opx <- tpf.recalc_minaral("OPX",Opx);
    Spl <- tpf.recalc_minaral("SPL",Spl);
    Ol <- tpf.recalc_minaral("OL",Ol);

    FOL <- Ol$Mg_pfu / (Ol$Mg_pfu + Ol$Fe_p2_pfu + Ol$Fe_p3_pfu);
    AOP <- Opx$Al_c6_pfu;
    POP <- (Opx$Mg_pfu / (Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu)) * (1 - Opx$Ca_pfu - Opx$Na_pfu - Opx$Mn_pfu);
    MSP <- Spl$Mg_pfu / (Spl$Mg_pfu + Spl$Fe_p2_pfu);
    ASP <- (Spl$Al_c4_pfu + Spl$Al_c6_pfu) / (Spl$Al_c4_pfu + Spl$Al_c6_pfu + Spl$Cr_pfu);
    KB <- FOL * FOL * AOP / (POP * MSP * ASP * ASP);
    T <- 3666 / (0.8808 - log(KB)) - 273.15;
    return(data.frame(Name=paste0(Spl$Name,"-", Ol$Name,"-", Opx$Name),T_calc=T));
}

sensors.t.ol_spl_opx_sen_1985.args <- function(){
    ars <- list();
    ars$minerals <- c("Ol", "Spl", "Opx");
    ars$ranges <- c();
    return(ars);
}

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