Grt-Cpx-Pl-Qtz: Fonarev, Unpubl

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

source("utilities_tpf.r");

sensors.p.grt_cpx_pl_qtz_fonarev_unpubl <- function(Grt, Cpx, Pl, Temp, ...){
    Name <- paste0(Pl$Name,"-", Cpx$Name,"-", Grt$Name);
    T <- Temp;
    Grt <- tpf.recalc_minaral("GRT2",Grt);
    Pl <- tpf.recalc_minaral("PL",Pl);
    Cpx <- tpf.recalc_minaral("CPX",Cpx);

    B <- Cpx$Mg_pfu / (Cpx$Mg_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + Cpx$Mn_pfu + Cpx$Ca_pfu + Cpx$Al_c6_pfu + Cpx$Na_pfu);
    A <- (Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu) / (Cpx$Mg_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + Cpx$Mn_pfu + Cpx$Ca_pfu + Cpx$Al_c6_pfu + Cpx$Na_pfu);
    C <- Cpx$Ca_pfu / (Cpx$Mg_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + Cpx$Mn_pfu + Cpx$Ca_pfu + Cpx$Al_c6_pfu + Cpx$Na_pfu);
    L <- Pl$Ca_pfu / (Pl$Ca_pfu + Pl$Na_pfu + Pl$K_pfu);
    M <- Pl$Na_pfu / (Pl$Ca_pfu + Pl$Na_pfu + Pl$K_pfu);
    X <- Grt$Fe_p2_pfu / (Grt$Fe_p2_pfu + Grt$Mg_pfu + Grt$Ca_pfu + Grt$Mn_pfu);
    Y <- Grt$Ca_pfu / (Grt$Fe_p2_pfu + Grt$Mg_pfu + Grt$Ca_pfu + Grt$Mn_pfu);
    Z <- Grt$Mg_pfu / (Grt$Fe_p2_pfu + Grt$Mg_pfu + Grt$Ca_pfu + Grt$Mn_pfu);
    E <- Cpx$TiO2_wtp + Cpx$MnO_wtp + Cpx$Al2O3_wtp + Cpx$Na2O_wtp + Cpx$Cpx$K_pfu2O_wtp + Cpx$Cr2O3_wtp;
    TT <- T + 273;
    AA <- 3 * log(X) + 3 * ((Z^2 * (0.2 + X * 4.6) + Y^2 * (3.5 - 0.0014 * TT - 7.1 * X) + Y * Z * (Y * 4.3       - Z * 1.8 - 0.7)) / 0.001987 / TT);
    AG <- 3 * log(Y) + 3 * ((Z^2 * (4 - 0.0015 * TT - Y * 6.1) + X^2 * ( - 0.1 - 0.0014 * TT + Y * 7.1) +       X * Z * (3.4 - 0.0029 * TT - 2.8 * X + 1.8 * Z)) / 0.001987 / TT);
    ACF <- log(A) + (C^2 * ( - 3.3 + 0.0062 * TT + (0.8 - 0.0007 * TT) * (4 * A - 1)) + B^2 * (11.9 - 0.0122 * TT +        ( - 8.4 + 0.0085 * TT) * (4 * A - 1)) + B * C * (11.9 - 0.0122 * TT - 3.3 + 0.0062 * TT -        16.4 + 0.0209 * TT + 2 - 0.0018 * TT) + 2 * A * B * C * (0.8 - 0.0007 * TT - 8.4 + 0.0085 * TT +        16.4 - 0.0209 * TT)) / 0.0019872 / TT;
    AP <- log(L) + (M^2 * ((1980 - 1.526 * TT) + 2 * L * (4880 - 2.348 * TT)) / 1.9872 / TT);
    K <- 3 * AP + 6 * ACF - AG - 2 * AA;
    P <- 12.7 + 0.0068 * TT - 11.48 * 1000 / TT - 1.09 * TT * K / 1000;
    if (E>5.0) {
    warning(" CANNOT BE USED FOR (TIO2+MNO+NA2O+K2O+CR2O3+AL2O3) IN CPX >5 WT.");
    }
    return(data.frame(Name=Name,Temp=Temp,P_calc=P));
}

sensors.p.grt_cpx_pl_qtz_fonarev_unpubl.args <- function(){
    ars <- list();
    ars$minerals <- c("Grt", "Cpx", "Pl");
    ars$ranges <- c("Temp");
    return(ars);
}

#sensors.p.grt_cpx_pl_qtz_fonarev_unpubl.test <- function(){
#
#}