====== Ol-Spl-Opx: Sen, 1985 ======
 
=== Details ===

<WRAP center round info 60%>
This sensor from [[:about_tpf|TPF]] collection
</WRAP>

^ Article |  |
^ Collibration |  |
^ Precision |  |
^ Ranges |  |
^ Renew sensors |  |

<use:sensors.t.ol_spl_opx_sen_1985>
<test:sensors.t.ol_spl_opx_sen_1985>

<WRAP center round todo 60%>
  * <del>заголовок</del>
  * краткое описание
  * статья
  * ссылка
  * карточка
  * подробное описание
  * <del>код</del>
  * тест
</WRAP>

=== Source ===
<Code linenums lang-r>
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(){
#
#}
</Code>

{{tag> tpf Spl Ol Opx }}