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

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

source("utilities_tpf.r");

sensors.f.mag_ilm_anderson_lindsley_1985_d <- 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;
    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 STORMER(1983) ");
    return(data.frame(Name=paste0(Ilm$Name,"-", Mag$Name),LGFO2_calc=LGFO2));
}

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

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