Opx-Cpx: Brey and Kohler, 1990

source('recalc_px_putirka.r');

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

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


Tcalc <- -273.15 + ((
                   23664 + (
                      24.9+126.3*Cpx$Fe_pfu/(Cpx$Fe_pfu+Cpx$Mg_pfu))*P
                      )/(
                      13.38 + (
                          log((1-Cpx$Ca_pfu/(1-Cpx$Na_pfu))/(1-Opx$Ca_pfu/(1-Opx$_Na_pfu)))
                      )^2 + 
                      11.59*Opx$Fe_pfu/(Opx$Fe_pfu+Opx$Mg_pfu)
                     ));

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

Ed_check <- Ed_check > 0.95 && Ed_check < 1.23;

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

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