====== Mag-Ilm: Anderson, Lindsley, 1985 <B> ======
 
=== Details ===

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

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

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

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

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

sensors.t.mag_ilm_anderson_lindsley_1985_b <- function(Mag, Ilm, ...){
    Ilm <- tpf.recalc_minaral("ILM",Ilm);
    Mag <- tpf.recalc_minaral("MAG",Mag);

    X1 <- Mag$Ti_pfu * Mag$Fe_p2_pfu / (Mag$Fe_p2_pfu + Mag$Mn_pfu + Mag$Mg_pfu + Mag$Ca_pfu);
    X2 <- (0.5 * Mag$Fe_p3_pfu * Mag$Fe_p3_pfu / (Mag$Fe_p3_pfu + Mag$Al_pfu + Mag$Cr_pfu) + Mag$Ti_pfu * Mag$Fe_p2_pfu / (Mag$Fe_p2_pfu + Mag$Mn_pfu + Mag$Mg_pfu + Mag$Ca_pfu));
    Y <- sqrt(Ilm$Fe_p2_pfu * Ilm$Ti_pfu) / (0.5 * Ilm$Fe_p3_pfu + sqrt(Ilm$Fe_p2_pfu * Ilm$Ti_pfu));
    X <- X1 / X2;
    A1 <-  - 3 * X^2 + 4 * X - 1;
    A2 <- 3 * X^2 - 2 * X;
    A3 <-  - 3 * Y^2 + 4 * Y - 1;
    A4 <- 3 * Y^2 - 2 * Y;
    K <- 8.314 * log(X * (1 - Y)^2 / Y^2 / (1 - X));
    T <- ( - A1 * 48050 - A2 * 32075 + A3 * 136900 + A4 * 33020 + 30782) /      ( - A1 * 41.4 - A2 * 30.3 + A3 * 111.965 + A4 * 2.5008 + 6.0655 - K) - 273;
    warning(" X(i)- after STORMER J.C.,JR.(1983)");
    return(data.frame(Name=paste0(Ilm$Name,"-", Mag$Name),T_calc=T));
}

sensors.t.mag_ilm_anderson_lindsley_1985_b.args <- function(){
    ars <- list();
    ars$minerals <- c("Mag", "Ilm");
    ars$ranges <- c();
    return(ars);
}

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

{{tag> tpf Ilm Mag }}