Grt-Opx: Perchuk, 1977

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

source("utilities_tpf.r");
 
sensors.p.grt_opx_perchuk_1977 <- function(Grt, Opx, Temp, ...){
T <- Temp;
    Opx <- tpf.recalc_minaral("OPX",Opx);
 
    N <- Opx$Al2O3_wtp;
    TT <- T + 273;
    if (N<4.000009) {
    P <- 15 + (26.5 - 0.045 * T) * (log(N) - 3.778 + 3150 / TT);
    }
    if (N>4.000009) {
    P <- 20 + (27.53 - 0.028 * T) * (log(N) - 4.175 + 3077 / TT);
    }
    warning(" FOR LHERZOLITES");
    return(data.frame(Name=paste0(Opx$Name),Temp=Temp,P_calc=P));
}
 
sensors.p.grt_opx_perchuk_1977.args <- function(){
    ars <- list();
    ars$minerals <- c("Grt", "Opx");
    ars$ranges <- c("Temp");
    return(ars);
}
 
#sensors.p.grt_opx_perchuk_1977.test <- function(){
#
#}