====== Opx-Cpx: Bertrand, Mercier, 1985 ======

This geothermometer is directly applicable to the natural system as illustrated by the remarkable agreement between the experimental conditions and the calculated temperatures for reequilibration experiments on natural systems. 

=== Details ===

^ Source | <ref:bertrand1985> |
^ Callibration | 36 points from different sources |
^ Accuracy | 2% |
^ Ranges | 810-1500°C and 0.0-10.0 GPa |
^ Renew sensors |  |

<use:sensors.t.opx_cpx_bertrand_mercier>
<test:sensors.t.opx_cpx_bertrand_mercier>

<WRAP center round todo 60%>
  * <del>заголовок</del>
  * <del>краткое описание</del>
  * <del>статья</del>
  * <del>ссылка</del>
  * <del>карточка</del>
  * подробное описание
  * <del>код</del>
  * тест
</WRAP>
=== Code ===

<Code linenums lang-r>
source("utilities_tpf.r");
 
sensors.t.opx_cpx_bertrand_mercier <- function(Opx, Cpx, P, usePFU=FALSE, ...){

    row_names <- paste0(Opx$Name,"-",Cpx$Name);

    if(!usePFU) {
        Opx <- tpf.recalc_minaral("OPX",Opx);
        Cpx <- tpf.recalc_minaral("CPX",Cpx);
    }

    X = Opx$Na_pfu;
    Y = Opx$Ca_pfu;
    Z = Cpx$Na_pfu;
    F = Cpx$Ca_pfu;

    K = Cpx$Fe_pfu /(Cpx$Fe_pfu + Cpx$Mg_pfu);
    
    A = 36273+399*P/10;
    C = 1-Y/(1-X);
    B = F/(1-Z)+(-0.77+0.001*T)*K;
    L = log((1-B)/C);
    T = A/(19.31-8.314*L-12.15*B**2);
    T = T-273;
    
    return(data.frame(Name=row_names,P=P, T_calc=T));
}

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

</Code>

{{tag> AP_igroup fixed}}