Гранат-клинопироксеновый термометр Д. Накамура, 2009

Это старая версия документа!

Термометр с учетом неидеальности Fe-Mg обмена в гранате.

Подробнее

Информация о публикации	D. Nakamura, A new formulation of garnet–clinopyroxene geothermometer based on accumulation and statistical analysis of a large experimental data set, J. metamorphic Geol., 2009, 27, 495–508
Ссылка	nakamura2009.pdf
Калибровка	Экспериментальные данные. 333 пары в графитовых капсулах
Точность	±74 ºС
Диапазон значений	800-1820 ºС
Обновляет сенсоры	grt_cpx_krogh_ravna

Код метода

source('recalc_grt_8cat_norm.r');
source('recalc_cpx_4cat_norm_jd.r');

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

    if (!usePFU) Grt <- recalc.grt_8cat_norm(Grt);
    if (!usePFU) Cpx <- recalc.cpx_4cat_norm_jd(Cpx);
    

    X_Grt_Ca <- Grt$Ca_pfu / (Grt$Mg_pfu + Grt$Fe_pfu + Grt$Ca_pfu);
    X_Grt_Fe <- Grt$Fe_pfu / (Grt$Mg_pfu + Grt$Fe_pfu + Grt$Ca_pfu);
    X_Grt_Mg <- Grt$Mg_pfu / (Grt$Mg_pfu + Grt$Fe_pfu + Grt$Ca_pfu);

    X_Cpx_Ca <- Cpx$Ca_pfu / (Cpx$Mg_pfu + Cpx$Fe_p2_pfu + Cpx$Ca_pfu);
    X_Cpx_Fe <- Cpx$Fe_pfu / (Cpx$Mg_pfu + Cpx$Fe_p2_pfu + Cpx$Ca_pfu);
    X_Cpx_Mg <- Cpx$Mg_pfu / (Cpx$Mg_pfu + Cpx$Fe_p2_pfu + Cpx$Ca_pfu);
    
    
    return (data.frame(Name=paste0(Grt$Name,'-',Cpx$Name),P = P, X_Grt_Ca = X_Grt_Ca, X_Grt_Mg = X_Grt_Mg, X_Grt_Fe = X_Grt_Fe, X_Cpx_Ca = X_Cpx_Ca, X_Cpx_Fe = X_Cpx_Fe, X_Cpx_Mg = X_Cpx_Mg));
}

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

sensors.t.grt_cpx_nakamura.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, 711, tolerance = 1);
}