Opx-Grt-Cpx: Raith, 1983

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

source("utilities_tpf.r");

sensors.t.opx_grt_cpx_raith_1983 <- function(Opx, Grt, Cpx, ...){
    Opx <- tpf.recalc_minaral("OPX",Opx);
    Grt <- tpf.recalc_minaral("GRT",Grt);

    X <- Opx$Fe_p2_pfu + Opx$Fe_p3_pfu;
    SUM1 <- Opx$Si_pfu * 2 + Opx$Al_c4_pfu * 3 / 2 + Opx$Al_c6_pfu * 3 / 2 + Opx$Cr_pfu * 3 / 2;
    SUM2 <- Opx$Ti_pfu * 2 + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu * 3 / 2 + Opx$Mn_pfu + Opx$Mg_pfu + Opx$Ca_pfu + Opx$Na_pfu / 2 + Opx$K_pfu / 2;
    S <- Opx$Si_pfu + Opx$Al_pfu + Opx$Cr_pfu + Opx$Ti_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mn_pfu + Opx$Mg_pfu + Opx$Ca_pfu + Opx$Na_pfu + Opx$K_pfu;
    Z <- Opx$Mg_pfu;
    Y <- (Grt$Fe_p2_pfu + Grt$Fe_p3_pfu) / Grt$Mg_pfu;
    SUM <- SUM1 + SUM2;
    SK <- 4 / S;
    SUM <- SK * SUM;
    SO <- 12 - SUM * 2;
    if (SO>0.00) {
    CFE3 <- SO;
    }
    if (SO<0.00) {
    CFE3 <- 0;
    }
    F <- X - CFE3;
    K <- log(Y / (F / Z));
    T <- 1648 / (K + 0.334) - 273;
    warning(" CAN BE USED ONLY FOR BASIC GRANULITE AND ECLOGITE, ACCURACY +/- 70 C");
    return(data.frame(Name=paste0(Grt$Name,"-", Opx$Name),T_calc=T));
}

sensors.t.opx_grt_cpx_raith_1983.args <- function(){
    ars <- list();
    ars$minerals <- c("Opx", "Grt", "Cpx");
    ars$ranges <- c();
    return(ars);
}

#sensors.t.opx_grt_cpx_raith_1983.test <- function(){
#
#}