====== Cpx-Opx: Kretz, 1982 <A> ======

Cpx slope of the solvus surface based on Natural data for associated Ca pyroxene (Cpx) and orthopyroxene (Opx) or pigeonite are combined with experimental data for Fe-free pyroxenes.

=== Details ===

<WRAP center round info 60%>
This sensor from [[:about_tpf|TPF]] collection
</WRAP>

^ Article | <ref:kretz1982> |
^ Callibration | Experimetal and data from Skaergaard and Bushveld |
^ Precision | ±60°C |
^ Ranges | 750-1400 °C |
^ See also | [[sensors:t:cpx_opx_kretz_1982_b|Kretz, 1982 <B>]] |

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

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

=== Source ===
<Code linenums lang-r>
source("utilities_tpf.r");

sensors.t.cpx_opx_kretz_1982_a <- function(Cpx, usePFU = FALSE, ...){
    
    row_names = Cpx$Name;
    
    if(!usePFU) 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;
    TL <- (1000 / (0.054 + 0.608 * X - 0.304 * log(1 - 2 * Y))) - 273.15;
    
    Ts <- data.frame(TH = TH, TL = TL, T = rep(NA, nrow(Cpx)));
    Ts$T <- ifelse(Ts$TH>1080.0,Ts$TH,ifelse(Ts$TL<1080.0,TL,NA));

    return(data.frame(Name=row_names,T_calc=Ts$T));
}

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

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

{{tag> tpf Cpx AP_igroup fixed}}