Bt-Opx: Aranovich, 1988

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

source("utilities_tpf.r");

sensors.t.bt_opx_aranovich_1988_input <- function(Bt, Opx, P, ...){
    Opx <- tpf.recalc_minaral("OPX",Opx);
    Bt <- tpf.recalc_minaral("BT",Bt);

    X <- (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu) / (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mg_pfu + Opx$Al_pfu / 2);
    Y <- Opx$Mg_pfu / (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mg_pfu + Opx$Al_pfu / 2);
    Z <- (Opx$Al_pfu / 2) / (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mg_pfu + Opx$Al_pfu / 2);
    M <- (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu) / (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu + Bt$Mg_pfu);
    N <- Bt$Mg_pfu / (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu + Bt$Mg_pfu);
    A <- (Bt$Al_c4_pfu + Bt$Al_c6_pfu - 1) / 2;
    W <- Bt$Ti_pfu;
    FF <- F;
    L <- FF / 2;
    T <- (3055 + 0.015 * P * 1000 - 1903 * (A + W) + 2372 * (Y - X) +      5204 * Z) / (1.987 * log(X * N / Y / M) + 2.531 - 0.287 * (A + W) + 1.688 * (Y - X) - 1.466 * L) - 273;
    warning(" THE QUANTITY OF 'F' IN BT IS USED. ");
    return(data.frame(Name=paste0(Bt$Name,"-", Opx$Name),P=P,T_calc=T));
}

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

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