Grt-Opx: Aranovich, Kosyakova, 1987 <B>

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

source("utilities_tpf.r");

sensors.p.grt_opx_aranovich_kosyakova_1987_b <- function(Grt, Opx, Temp, ...){
T <- Temp;
    Opx <- tpf.recalc_minaral("OPX",Opx);
    Grt <- tpf.recalc_minaral("GRT",Grt);

    A <- Grt$Al_pfu / (Grt$Al_pfu + Grt$Cr_pfu + Grt$Ti_pfu + Grt$Fe_p3_pfu);
    B <- Grt$Cr_pfu / (Grt$Al_pfu + Grt$Cr_pfu + Grt$Ti_pfu + Grt$Fe_p3_pfu);
    C <- Grt$Mg_pfu / (Grt$Mg_pfu + Grt$Fe_p2_pfu + Grt$Mn_pfu + Grt$Ca_pfu + Grt$Na_pfu + Grt$K_pfu);
    D <- Grt$Fe_p2_pfu / (Grt$Mg_pfu + Grt$Fe_p2_pfu + Grt$Mn_pfu + Grt$Ca_pfu + Grt$Na_pfu + Grt$K_pfu);
    F <- Grt$Ca_pfu / (Grt$Mg_pfu + Grt$Fe_p2_pfu + Grt$Mn_pfu + Grt$Ca_pfu + Grt$Na_pfu + Grt$K_pfu);
    X1 <- ((Opx$Al_c4_pfu + Opx$Al_c6_pfu - abs(Opx$Na_pfu - Opx$Cr_pfu) - 2 * Opx$Ti_pfu) / 2);
    X2 <- ((Opx$Al_c4_pfu + Opx$Al_c6_pfu - abs(Opx$Na_pfu - Opx$Cr_pfu) - 2 * Opx$Ti_pfu) / 2 + Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Ca_pfu);
    Y <- Opx$Mg_pfu / (Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Ca_pfu + (Opx$Al_c4_pfu + Opx$Al_c6_pfu - abs(Opx$Na_pfu - Opx$Cr_pfu) - 2 * Opx$Ti_pfu) / 2);
    Z <- Opx$Fe_p2_pfu / (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mg_pfu + Opx$Ca_pfu + (Opx$Al_c4_pfu + Opx$Al_c6_pfu - abs(Opx$Na_pfu - Opx$Cr_pfu) - 2 * Opx$Ti_pfu) / 2);
    X <- X1 / X2;
    TT <- T + 273.15;
    G <- log(Y^3 * X / C^3 / A^2);
    U <- log(D * Y / C / Z);
    W <- 1.987 * U + 0.49 * B - 1.688 * (Y - Z) - 1.34 * X + 2.654 + F * (3 * C + 3.436 * D + 1.466 *      F + 0.658 * (D - C));
    J <- 1237 * (3 * X + Y - 4 * X * Y) + 6441 * Z * (1 - 4 * X);
    I <- (0.3285 - 0.14245 * (3 * X^2 + (Z + Y)^2)) * 1000;
    S <- 1218 * B - 2372 * (Y - Z) - 5204 * X + F * (6600 * C + 7322 * D + 3226 * F + 1506.5 * (D - C));
    R <- 3 * F * F * (1 - 2 * C);
    RR <- 6 * C * F * (1 - C);
    P <- (R * (3300 - 1.5 * TT) + RR * (3300 - 1.5 * TT) -      6 * F * D * (F * (74 - 0.034 * TT) + D * ( - 3661 + 1.718 * TT)) + 6 * D * F * (1 - 2 * C) *      (1506.5 - 0.658 * TT) + (1 - C) * B * (29229 - 11.88 * TT) - D * B * (25575 - 10.4 * TT) -      Z * (3 - 4 * Y) * ( - 2372 + 1.688 * TT) + J - 1.987 * TT * G - 2785 - 1.869 * TT - 5.88 * (TT - 970) +      5.88 * TT * log(TT / 970)) / I;
    warning(" ACCURACY +-1-8 KBAR");
    return(data.frame(Name=paste0(Grt$Name,"-", Opx$Name),Temp=Temp,P_calc=P));
}

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

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