Opx-Cpx: Putirka, 2009 RiMG

source('recalc_px_putirka.r');
source('interpmin_opx_thompson.r');
source('interpmin_cpx.r');

sensors.t.opx_cpx_putirka <- function(Opx, Cpx, P=10, usePFU = FALSE, ...){

if (!usePFU) Opx <- recalc.px_putirka(Opx);
if (!usePFU) Cpx <- recalc.px_putirka(Cpx);

interpmin.opx_thompson

Tcalc <- 10^4/(11.2-1.96*log(CJ17/EO17)-3.3*BX17-25.8*CH17+33.2*EA17-23.6*ED17-2.08*EP17-8.33*EQ17-0.05*P);

Eq_check <- (Cpx$Fe_pfu/Cpx$Mg_pfu)/(Opx$Fe_pfu/Opx$Mg_pfu);

Eq_check <- Eq_check > 0.95 && Eq_check < 1.23;

return (data.frame(Name=paste0(Opx$Name,'-',Cpx$Name),P_seq=P,T_calc_Opx_Cpx_Putirka=Tcalc,Eq_check=Eq_check));
}

sensors.t.opx_cpx_putirka.args <- function() {
    ars <- list();
    ars$minerals = c('Opx','Cpx');
    ars$ranges = c('P');
    return(ars);
}