====== Grt-Bt: Nikitina, 1991 ======
 
=== Details ===

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

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

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

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

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

sensors.t.grt_bt_nikitina_1991 <- function(Grt, Bt, ...){
    Bt <- tpf.recalc_minaral("BT",Bt);
    Grt <- tpf.recalc_minaral("GRT",Grt);

    CAG <- Grt$Ca_pfu / 3;
    FMG <- (Grt$Fe_p2_pfu + Grt$Mn_pfu) / 3;
    FMB <- (Bt$Fe_p2_pfu + Bt$Mn_pfu) / 3;
    ALB <- Bt$Al_c6_pfu / 3;
    TIB <- Bt$Ti_pfu / 3;
    D <-  - 4.08 * (ALB - 0.15) - 21.6 * (ALB - 0.15)^2 + 634.5 * (ALB - 0.15)^3 + 3628.8 * (ALB - 0.15)^4;
    E <- 1.8 * (TIB - 0.033) - 1.8 * (TIB - 0.033)^2;
    A <- 1 - FMG;
    C <-  - 2.55 * (CAG - 0.1) - 2.7 * (CAG - 0.1)^2 + 35.4 * (CAG - 0.1)^3;
    T <- ( - 14811.25 * A + 60034.08 * A^2 - 97177.13 * A^3 + 46500 * A^4) /      (log(FMB / FMG) - 10.946 * A + 48.238 * A^2 - 81.784 * A^3 + 43.802 * A^4 + C + D + E) - 273;
    return(data.frame(Name=paste0(Grt$Name,"-", Bt$Name),T_calc=T));
}

sensors.t.grt_bt_nikitina_1991.args <- function(){
    ars <- list();
    ars$minerals <- c("Grt", "Bt");
    ars$ranges <- c();
    return(ars);
}

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

{{tag> tpf Grt Bt }}