====== Mag-Ilm: Powell, Powell, 1977 <A> ======
 
=== 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_powell_powell_1977_a>
<test:sensors.t.mag_ilm_powell_powell_1977_a>

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

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

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

    XMA <- Mag$Ti_pfu;
    XMI <- Mag$Ti_pfu - 0.5 * (Mag$Mg_pfu + Mag$Mn_pfu + Mag$Ca_pfu);
    YMI <- Ilm$Ti_pfu - (Ilm$Mg_pfu + Ilm$Mn_pfu + Ilm$Ca_pfu);
    YMA <- Ilm$Ti_pfu - (Ilm$Mg_pfu + Ilm$Mn_pfu + Ilm$Ca_pfu);
    ZMI <- Mag$Fe_p3_pfu / 2 - (Mag$Mg_pfu + Mag$Mn_pfu + Mag$Ca_pfu);
    ZMA <- Mag$Fe_p3_pfu / 2;
    HMI <- 0.5 * (Ilm$Fe_p3_pfu - Ilm$Al_pfu - Ilm$Cr_pfu);
    HMA <- Ilm$Fe_p3_pfu / 2;
    KMA <- log(XMA * HMA / ZMI / YMI);
    KMI <- log(XMI * HMI / ZMA / YMA);
    TMA <- 8155 / (4.59 - KMA) - 273;
    TMI <- 8155 / (4.59 - KMI) - 273;
    T <- (TMI + TMA) / 2;
    return(data.frame(Name=paste0(Ilm$Name,"-", Mag$Name),T_calc=T));
}

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

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

{{tag> tpf Ilm Mag }}