Mag-Ilm: Anderson, Lindsley, 1985 <C>

Использовать Проверить

source("utilities_tpf.r");

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

    X <- (Mag$Ti_pfu - Mag$Mn_pfu / 2) / (Mag$Ti_pfu - Mag$Mn_pfu / 2 + Mag$Fe_p3_pfu / 2);
    Y <- (Ilm$Ti_pfu - Ilm$Mn_pfu + Ilm$Al_pfu / 2) / (Ilm$Ti_pfu - Ilm$Mn_pfu + Ilm$Al_pfu / 2 + Ilm$Fe_p3_pfu / 2);
    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;
    TT <- T + 273;
    LgFO2 <-  - 23847.6 / TT + 13.48 + (12 * log(1 - Y) - 4 * log(1 - X) + (8 * X^2 * (X - 1) * (48050 - TT * 41.4) +      4 * X^2 * (1 - 2 * X) * (32075 - TT * 30.3) + 12 * Y^2 * (1 - Y) * (136900 - TT * 111.965) -      6 * Y^2 * (1 - 2 * Y) * (33020 - TT * 2.5008)) / 8.314 / TT) / 2.303;
    warning(" THERMOOXOMETER, X(i)- after LINDSLEY ,SPENCER (1982)");
    return(data.frame(Name=paste0(Ilm$Name,"-", Mag$Name),LGFO2_calc=LGFO2));
}

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

#sensors.f.mag_ilm_anderson_lindsley_1985_c.test <- function(){
#
#}