====== Mag-Ilm: Powell, Powell, 1977 <B> ======
 
=== Details ===

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

^ Article |  |
^ Collibration |  |
^ Presition |  |
^ Ranges |  |
^ Renew sensors |  |

<use:sensors.f.mag_ilm_powell_powell_1977_b>
<test:sensors.f.mag_ilm_powell_powell_1977_b>

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

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

sensors.f.mag_ilm_powell_powell_1977_b <- 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);
    FMA <- (2 * log(ZMI) + 6 * log(YMI) - 6 * log(XMA) + 10.702 * KMA - 4.94) / 2.30258;
    FMI <- (2 * log(ZMA) + 6 * log(YMA) - 6 * log(XMI) + 10.702 * KMI - 4.94) / 2.30258;
    LGFO2 <- (FMA + FMI) / 2;
    return(data.frame(Name=paste0(Ilm$Name,"-", Mag$Name),LGFO2_calc=LGFO2));
}

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

#sensors.f.mag_ilm_powell_powell_1977_b.test <- function(){
#
#}
</Code>

{{tag> tpf Ilm Mag }}