Cpx-L: Nielsen, Drake, 1979

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

source("utilities_tpf.r");

sensors.t.cpx_l_nielsen_drake_1979 <- function(Cpx, L, ...){
    L <- tpf.recalc_minaral("L",L);
    Cpx <- tpf.recalc_minaral("CPX",Cpx);

    C <- Cpx$Ca_pfu / (Cpx$Ca_pfu + Cpx$Mg_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + 2 * Cpx$Na_pfu + Cpx$Cr_pfu + Cpx$Ti_pfu + 0.5 * (Cpx$Al_c4_pfu + Cpx$Al_c6_pfu - Cpx$Na_pfu));
    M <- Cpx$Mg_pfu / (Cpx$Ca_pfu + Cpx$Mg_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + 2 * Cpx$Na_pfu + Cpx$Cr_pfu + Cpx$Ti_pfu + 0.5 * (Cpx$Al_c4_pfu + Cpx$Al_c6_pfu - Cpx$Na_pfu));
    F1 <- (Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu);
    F2 <- (Cpx$Ca_pfu + Cpx$Mg_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + 2 * Cpx$Na_pfu + Cpx$Cr_pfu + Cpx$Ti_pfu + 0.5 * (Cpx$Al_c4_pfu + Cpx$Al_c6_pfu - Cpx$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 <- 14600 / (K + 9.52) - 273;
    return(data.frame(Name=paste0(L$Name,"-", Cpx$Name),T_calc=T));
}

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

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