====== Opx-Cpx: Putirka, 2009 ======

Global regression, based on the partitioning of enstatite + ferrosilite (Fm2Si2O6 = EnFs; FmO = FeO + MgO + MnO) between clinopyroxene and orthopyroxene. Cpx with Mg#>0.75.

=== Details ===

^ Reference | <ref:putirka2008> |
^ Collibration | 483 points from different sources |
^ Precision | ± 60 °C  |
^ Ranges | P < 70 kbar,  700 °C < T < 2000 °C |
^ Renew sensors | [[:sensors:t:opx_cpx_brey_kohler|Brey and Köhler (1990) model T(BKN)]] |
<use:sensors.t.opx_cpx_putirka>
<test:sensors.t.opx_cpx_putirka>

<WRAP center round todo 60%>
Для улучшения этой статьи следует:
  * Добавить подробное описание
  * Добавить тест
</WRAP>


=== Source ===

<Code linenums lang-r>
source('recalc_px_putirka.r');
source('interpmin_opx_thompson.r');
source('interpmin_cpx_putirka.r');

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

stopifnot(usePFU == F);

CpxCmp <-interpmin.cpx_putirka(Cpx);
Cpx <- recalc.px_putirka(Cpx);

OpxCmp <- interpmin.opx_thompson(Opx);
Opx <- recalc.px_putirka(Opx);

Tcalc <- 10^4/(11.2-1.96*log(CpxCmp$X_EnFs/OpxCmp$X_Fm2Si2O6)-3.3*Cpx$Ca_pfu-25.8*CpxCmp$X_CrCaTs+33.2*Opx$Mn_pfu-23.6*Opx$Na_pfu-2.08*OpxCmp$X_En-8.33*OpxCmp$X_Di-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);
}
</Code>