Cpx-Opx: Sen, Jones, 1989

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

source("utilities_tpf.r");

sensors.t.cpx_opx_sen_jones_1989 <- function(Cpx, Opx, usePFU=FALSE, ...){

    row_names = paste0(Cpx$Name,"-", Opx$Name);

    if(!usePFU) {
        Opx <- tpf.recalc_minaral("OPX",Opx);
        Cpx <- tpf.recalc_minaral("CPX",Cpx);
    }

    A1 <- (1 - (Cpx$Al_c4_pfu + Cpx$Al_c6_pfu) / 2 - Cpx$Ti_pfu - Cpx$Cr_pfu - Cpx$Fe_p3_pfu) * Cpx$Mg_pfu;
    A2 <- (1 - (Cpx$Al_c4_pfu + Cpx$Al_c6_pfu) / 2 - Cpx$Ti_pfu - Cpx$Cr_pfu - Cpx$Fe_p3_pfu + 1 - Cpx$Mn_pfu - Cpx$Ca_pfu - Cpx$Na_pfu);
    A3 <- (1 - Cpx$Mn_pfu - Cpx$Ca_pfu - Cpx$Na_pfu) * Cpx$Mg_pfu;
    A4 <- (1 - (Cpx$Al_c4_pfu + Cpx$Al_c6_pfu) / 2 - Cpx$Ti_pfu - Cpx$Cr_pfu - Cpx$Fe_p3_pfu + 1 - Cpx$Mn_pfu - Cpx$Ca_pfu - Cpx$Na_pfu);
    B1 <- (1 - (Opx$Al_c4_pfu + Opx$Al_c6_pfu) / 2 - Opx$Ti_pfu - Opx$Cr_pfu - Opx$Fe_p3_pfu) * Opx$Mg_pfu;
    B2 <- (1 - (Opx$Al_c4_pfu + Opx$Al_c6_pfu) / 2 - Opx$Ti_pfu - Opx$Cr_pfu - Opx$Fe_p3_pfu + 1 - Opx$Mn_pfu - Opx$Ca_pfu - Opx$Na_pfu);
    B3 <- (1 - Opx$Mn_pfu - Opx$Ca_pfu - Opx$Na_pfu) * Opx$Mg_pfu;
    B4 <- (1 - (Opx$Al_c4_pfu + Opx$Al_c6_pfu) / 2 - Opx$Ti_pfu - Opx$Cr_pfu - Opx$Fe_p3_pfu + 1 - Opx$Mn_pfu - Opx$Ca_pfu - Opx$Na_pfu);
    A <- A1 / A2 * A3 / A4;
    B <- B1 / B2 * B3 / B4;
    T_eq3 <- 4900 / (1.807 - log(A / B)) - 273;


    A <- (1 - (Cpx$Al_c6_pfu) - Cpx$Ti_pfu - Cpx$Cr_pfu - Cpx$Fe_p3_pfu) * Cpx$Mg_pfu / (Cpx$Mg_pfu + Cpx$Fe_p2_pfu) * Cpx$Ca_pfu;
    B <- (1 - (Opx$Al_c6_pfu) - Opx$Ti_pfu - Opx$Cr_pfu - Opx$Fe_p3_pfu) * Opx$Mg_pfu / (Opx$Mg_pfu + Opx$Fe_p2_pfu) * Opx$Ca_pfu;
    T_eq4 <- 7045 / (2.47 - log(B / A)) - 273;

    return(data.frame(Name=row_names,T_calc_eq3=T_eq3,T_calc_eq4=T_eq4));
}

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

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