sensors:t:grt_cpx_krogh_ravna [pin.ifz.ru]

Эта страница только для чтения. Вы можете посмотреть её исходный текст, но не можете его изменить. Сообщите администратору, если считаете, что это неправильно.
====== Grt-Cpx: Ravna, 2000 ======

Калиброван на 360 образцах с широким покрытием давлений температур и составов из 27 экспериментальных наборов данных и 49 природных высоко-Mn гранулитах из Ruby Range, Montana, USA, и Karnataka, India с коэффициентом корреляции: (r2=0.98).

Равновесие Fe2+–Mg между гранатом и клинопироксеном не подвержено влиянию Na в клинопироксене в диапазоне X_Cpx_Na = 0 - 0.51.

=== Details ===


^ Article | <ref:ravna2000>  |
^ Calibration | 360 образцов с широким покрытием давлений температур и составов из 27 экспериментальных наборов данных и 49 природных высоко-Mn гранулитах из Ruby Range, Montana, USA, и Karnataka, India |
^ Accuracy | ±83 ºС |
^ Ranges | 600-1740 ºС, 1-16 kbar |
^ Renew | Berman et al. (1995), Ai (1994), Krogh (1988), Powell (1985), Ellis & Green (1979) |

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

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

=== Код метода ===

<Code linenums lang-r>
source('recalc_grt_12ox_norm.r');
source('recalc_cpx_4cat_norm.r');
source('recalc_common_cation.r');

sensors.t.grt_cpx_krogh_ravna <- function(Grt, Cpx, P, usePFU = FALSE, ...){

    Name <- paste0(Grt$Name,"-",Cpx$Name);

    if (!usePFU) Grt <- recalc.grt_12ox_norm(Grt);
    if (!usePFU) Cpx <- recalc.cpx_4cat_norm(Cpx);

    X_Grt_Ca <- Grt$Ca_pfu / (Grt$Ca_pfu + Grt$Mn_pfu + Grt$Fe_p2_pfu + Grt$Mg_pfu);
    X_Grt_Mn <- Grt$Mn_pfu / (Grt$Ca_pfu + Grt$Mn_pfu + Grt$Fe_p2_pfu + Grt$Mg_pfu);
    X_Grt_Mg <- Grt$Mg_pfu / (Grt$Mg_pfu + Grt$Fe_p2_pfu);

    Kd <- (Grt$Fe_p2_pfu / Grt$Mg_pfu) / (Cpx$Fe_p2_pfu / Cpx$Mg_pfu);
    
    Tcalc <- ( 1939.9 +
               3270 * X_Grt_Ca - 1396 * X_Grt_Ca^2 +
               3319 * X_Grt_Mn - 3535 * X_Grt_Mn^2 +
               1105 * X_Grt_Mg - 3561 * X_Grt_Mg^2 + 2324 * X_Grt_Mg^3 +
               169.4 * (P/10) ) / (log(Kd)+1.223) - 273 ;

    return (data.frame(Name=Name,P = P, T_calc_grt_cpx_krogh_ravna=Tcalc));
}

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

sensors.t.grt_cpx_krogh_ravna.test <- function () {
  
  Grt <- data.frame(Name='Grt',SiO2_wtp=37.50,TiO2_wtp=0.03, Al2O3_wtp=20.21, FeO_wtp=26.3, MnO_wtp=2.94, MgO_wtp=3.25, CaO_wtp=9.45);
  Cpx <- data.frame(Name='Cpx',SiO2_wtp=50.95,TiO2_wtp=0.23, Al2O3_wtp=2.42, FeO_wtp=12.35, MnO_wtp=0.33, MgO_wtp=10.93, CaO_wtp=21.39, Na2O=0.62);

  result <- sensors.t.grt_cpx_krogh_ravna (Grt, Cpx, 10);
  
  expect_equal(result$T_calc_grt_cpx_krogh_ravna, 711, tolerance = 1);
}
</Code>