Opx-Mag-Qtz: Graphchikov, 1992

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

source("utilities_tpf.r");

sensors.f.opx_mag_qtz_graphchikov_1992 <- function(Opx, Mag, Qtz, Temp, P, ...){
T <- Temp;
    Opx <- tpf.recalc_minaral("OPX",Opx);

    X <- (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu) / (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mg_pfu + Opx$Ca_pfu + Opx$Mn_pfu + Opx$Al_c6_pfu);
    Y <- Opx$Mg_pfu / (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mg_pfu + Opx$Ca_pfu + Opx$Mn_pfu + Opx$Al_c6_pfu);
    Z <- Opx$Mn_pfu / (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mg_pfu + Opx$Ca_pfu + Opx$Mn_pfu + Opx$Al_c6_pfu);
    M <- Opx$Mn_pfu / (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mg_pfu + Opx$Mn_pfu);
    V <-  - 32.544 + 0.038462 * (T + 273.15);
    if (P<V) {
    TT <- T + 273.15;
    FO2 <- 2.606 * ( - log(X) - (( - 4350 + 5051 * 1000 / TT) * Y^2 + 2400 * Z^2 + Y * Z * ( - 1950          + 5051 * 1000 / TT)) /          1.9872 / TT - ( - 3.2637 + 9.7534 * 1000 / TT + 0.05675 * (1000 / TT) * (1000 / TT)));
    LgFO2 <- FO2 + 2.606 * (0.1281 * P * 1000 / 1.9872 / TT);
    }
    if (P>V) {
    TT <- T + 273.15;
    FO2 <-  - log(X) - (Y * Y * ( - 4350 + 5051 * 1000 / TT) + 2400 * (Z^2) + Y * Z * ( - 1950     + (5051 * 1000) / TT)) / (1.9872 * TT)     - ( - 3.5115 + 9.963 * 1000 / TT + 0.05675 * (1000 / TT) * (1000 / TT));
    LgFO2 <- 2.606 * (FO2 + 0.1153 * P * 1000 / (1.9872 * TT));
    }
    return(data.frame(Name=paste0(Opx$Name),Temp=Temp,P=P,LGFO2_calc=LGFO2));
}

sensors.f.opx_mag_qtz_graphchikov_1992.args <- function(){
    ars <- list();
    ars$minerals <- c("Opx", "Mag", "Qtz");
    ars$ranges <- c("Temp", "P");
    return(ars);
}

#sensors.f.opx_mag_qtz_graphchikov_1992.test <- function(){
#
#}