====== Spr-Spl: Owen, Greenough, 1991 ======
 
=== Details ===

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

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

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

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

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

sensors.t.spr_spl_owen_greenough_1991 <- function(Spr, Spl, ...){
    Spr <- tpf.recalc_minaral("SPR",Spr);
    Spl <- tpf.recalc_minaral("SPL",Spl);

    A <- (Spr$Al_c4_pfu + Spr$Al_c6_pfu) * 2;
    S <- Spr$Si_pfu * 2;
    CT <- (Spr$Cr_pfu + 2 * Spr$Ti_pfu) * 2;
    F <- (Spr$Fe_p2_pfu + Spr$Fe_p3_pfu) * 2;
    M <- Spr$Mg_pfu * 2;
    Y <- Spl$Fe_p2_pfu / Spl$Mg_pfu;
    A4 <- 6 - S;
    A6 <- A - A4;
    if (A4>(A6 + CT)) {
    F3 <- A4 - A6 - CT;
    F2 <- F - F3;
    }
    if (A4<(A6 + CT)) {
    F2 <- F;
    }
    KD <- Y / (F2 / M);
    T <- 800 + 228 * log(KD) - 273.15;
    return(data.frame(Name=paste0(Spr$Name,"-", Spl$Name),T_calc=T));
}

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

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

{{tag> tpf Spr Spl }}