Opx-L: Nielsen, Drake, 1979

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

source("utilities_tpf.r");

sensors.t.opx_l_nielsen_drake_1979 <- function(Opx, L, ...){
    Opx <- tpf.recalc_minaral("OPX",Opx);
    L <- tpf.recalc_minaral("L",L);

    C <- Opx$Ca_pfu / (Opx$Ca_pfu + Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + 2 * Opx$Na_pfu + Opx$Cr_pfu + Opx$Ti_pfu + 0.5 * (Opx$Al_c4_pfu + Opx$Al_c6_pfu - Opx$Na_pfu));
    M <- Opx$Mg_pfu / (Opx$Ca_pfu + Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + 2 * Opx$Na_pfu + Opx$Cr_pfu + Opx$Ti_pfu + 0.5 * (Opx$Al_c4_pfu + Opx$Al_c6_pfu - Opx$Na_pfu));
    F1 <- (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu);
    F2 <- (Opx$Ca_pfu + Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + 2 * Opx$Na_pfu + Opx$Cr_pfu + Opx$Ti_pfu + 0.5 * (Opx$Al_c4_pfu + Opx$Al_c6_pfu - Opx$Na_pfu));
    X <- L$Ca_pfu / (L$Ca_pfu + L$Mg_pfu + L$Fe_p2_pfu + L$Fe_p3_pfu + L$Al_c4_pfu + L$Al_c6_pfu - L$Na_pfu - L$K_pfu);
    Y <- L$Mg_pfu / (L$Ca_pfu + L$Mg_pfu + L$Fe_p2_pfu + L$Fe_p3_pfu + L$Al_c4_pfu + L$Al_c6_pfu - L$Na_pfu - L$K_pfu);
    Z <- (L$Fe_p2_pfu + L$Fe_p3_pfu) / (L$Ca_pfu + L$Mg_pfu + L$Fe_p2_pfu + L$Fe_p3_pfu + L$Al_c4_pfu + L$Al_c6_pfu - L$Na_pfu - L$K_pfu);
    S <- L$Si_pfu / (L$Na_pfu + L$K_pfu + L$Si_pfu);
    F <- F1 / F2;
    K <- log(C * M * F / X / Y / Z / S);
    T <- 37800 / (K + 26.25) - 273;
    return(data.frame(Name=paste0(L$Name,"-", Opx$Name),T_calc=T));
}

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

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