====== Ol-Spl: Hill, Sack, 1987 ======
 
=== 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_hill_sack_1987>
<test:sensors.t.ol_spl_hill_sack_1987>

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

=== Source ===
<Code linenums lang-r>
source("utilities_tpf.r");

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

    Y <- (Ol$Fe_p2_pfu + Ol$Fe_p3_pfu) / (Ol$Fe_p2_pfu + Ol$Fe_p3_pfu + Ol$Mg_pfu);
    Y1 <- Ol$Mg_pfu / (Ol$Fe_p2_pfu + Ol$Fe_p3_pfu);
    X1 <- Spl$Mg_pfu / Spl$Fe_p2_pfu;
    X3 <- Spl$Cr_pfu / 2;
    X4 <- Spl$Ti_pfu;
    X5 <- Spl$Fe_p3_pfu / 2;
    A <- Spl$Al_pfu;
    if (A>0.20) {
    warning("CANNOT BE USED FOR AL IN SPL >0.2");
    }
    if (A<0.20) {
    K <- 0.0019872 * log(Y1 / X1);
    T <- (1.33795 + 4 * 80 * X3 + 3.5561 * X4 + 6.39 * X5 + 0.84 * (1 - 2 * Y)) /        (K + 0.0009539 - 0.0015698 * X4) - 273;
    }
    return(data.frame(Name=paste0(Spl$Name,"-", Ol$Name),T_calc=T));
}

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

#sensors.t.ol_spl_hill_sack_1987.test <- function(){
#
#}
</Code>

{{tag> tpf Spl Ol }}