Cpx-Opx: Kretz, 1982 <A>

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

source("utilities_tpf.r");

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

    X <- Cpx$Fe_p2_pfu / (Cpx$Fe_p2_pfu + Cpx$Mg_pfu);
    Y <- Cpx$Ca_pfu / (Cpx$Ca_pfu + Cpx$Mg_pfu + Cpx$Fe_p2_pfu);
    TH <- (1000 / (0.468 + 0.246 * X - 0.123 * log(1 - 2 * Y))) - 273.15;
    if (TH>1080.0) {
    T <- TH;
    warning(" ONLY CPX COMPOSITION IS INVOLVED IN CALCULATIONS, ACCURACY: +/-60 C");
    }
    TL <- (1000 / (0.054 + 0.608 * X - 0.304 * log(1 - 2 * Y))) - 273.15;
    if (TL<1080.0) {
    T <- TL;
    warning(" ONLY CPX COMPOSITION IS INVOLVED IN CALCULATIONS, ACCURACY: +/-60 C");
    }
    if (TH<1080.0) {
    if (TL>1080.0) {
    warning(" IS NOT TO BE USED FOR THESE COMPOSITIONS");
    }
    }
    return(data.frame(Name=paste0(Cpx$Name),T_calc=T));
}

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

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