Различия

Здесь показаны различия между двумя версиями данной страницы.

--- sensors:t:grt_cpx_nakamura [2015/09/22 12:16]
admin
+++ sensors:t:grt_cpx_nakamura [2018/11/12 15:07] (текущий)
@@ Строка 1: / Строка 1: @@
-====== Гранат-клинопироксеновый термометр Д. Накамура, 2009 ======
+====== Grt-Cpx: Nakamura, 2009 ======
 Термометр с учетом неидеальности Fe-Mg обмена в гранате.
@@ Строка 5: / Строка 5: @@
 === Подробнее ===
-^ Информация о публикации | 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 |
+^ Reference: | <ref:nakamura2009> |
-^ Ссылка | {{:sources:nakamura2009.pdf}} |
+^ Collibration: | Экспериментальные данные. 333 пары в графитовых капсулах |
-^ Калибровка | Экспериментальные данные. 333 пары в графитовых капсулах |
+^ Accuracy: | ±74 ºС |
-^ Точность | ±74 ºС |
+^ Ranges: | 800-1820 ºС 15-75 Kbar |
-^ Диапазон значений | 800-1820 ºС |
+^ Renew: | [[sensors:t:grt_cpx_krogh_ravna|Krogh Ravna, 2000]] |
-^ Обновляет сенсоры | [[sensors:t:grt_cpx_krogh_ravna]] |
+<use:sensors.t.grt_cpx_nakamura>
+<test:sensors.t.grt_cpx_nakamura>
+<WRAP center round todo 60%>
+Для улучшения этой статьи следует:
+  * Добавить подробное описание
+</WRAP>
 === Код метода ===
-<code>
+<Code linenums lang-r>
 source('recalc_grt_8cat_norm.r');
 source('recalc_cpx_4cat_norm_jd.r');
@@ Строка 28: / Строка 33: @@
     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);
+    X_Cpx_Fe <- Cpx$Fe_p2_pfu / (Cpx$Al_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + Cpx$Mg_pfu);
+    X_Cpx_Mg <- Cpx$Mg_pfu / (Cpx$Al_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + Cpx$Mg_pfu);
+    Kd <- (Grt$Fe_pfu / Grt$Mg_pfu) / (Cpx$Fe_p2_pfu / Cpx$Mg_pfu);
+    X_Prp <- Grt$Mg_pfu / (Grt$Fe_pfu + Grt$Mn_pfu + Grt$Mg_pfu + Grt$Ca_pfu);
+    X_Alm <- Grt$Fe_pfu / (Grt$Fe_pfu + Grt$Mn_pfu + Grt$Mg_pfu + Grt$Ca_pfu);
+    X_Sps <- Grt$Mn_pfu / (Grt$Fe_pfu + Grt$Mn_pfu + Grt$Mg_pfu + Grt$Ca_pfu);
+    X_Grs <- Grt$Ca_pfu / (Grt$Fe_pfu + Grt$Mn_pfu + Grt$Mg_pfu + Grt$Ca_pfu);
+    A <- 0.5*X_Grs*(X_Prp - X_Alm - X_Sps);
+    B <- 0.5*X_Grs*(X_Prp - X_Alm + X_Sps);
+    C <- 0.5*(X_Grs + X_Sps)*(X_Prp - X_Alm);
+    Tcalc <- (  2784 + 14.52 * P +
+               (2610 + 1.44*P)*(2*X_Grs*X_Prp - A) +
+               (1183 + 6.98*P)*(X_Grs^2 - A) -
+*(2*X_Grs*X_Alm + B) +
+               (814.6 + 3.61*P)*(X_Grs^2 + B) -
+               (254.6 + 8.42*P)*(2*X_Prp*X_Alm - X_Alm^2 + C) -
+.6*(X_Prp^2 - 2*X_Prp*X_Alm + C) +
+*X_Sps -
+*(X_Cpx_Mg - X_Cpx_Fe)  )  /
+             ( log(Kd) + 1.431 + 0.695*(2*X_Grs*X_Prp + X_Grs^2 - 2*A) +
+.203 * (X_Grs^2 - 2*X_Grs*X_Alm) +
+.922 * X_Sps ) -
+;
-    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));
+    return (data.frame(Name=paste0(Grt$Name,'-',Cpx$Name),P = P, T_calc_GrtCpx_Nakamura = Tcalc));
 }
@@ Строка 46: / Строка 75: @@
 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);
+  Grt <- data.frame(Name='Grt',SiO2_wtp=40.23,TiO2_wtp=0.55, Al2O3_wtp=21.87, FeO_wtp=19.03, MnO_wtp=0.40, MgO_wtp=10.07, CaO_wtp=7.27);
-  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);
+  Cpx <- data.frame(Name='Cpx',SiO2_wtp=52.86,TiO2_wtp=1.11, Al2O3_wtp=10.24, FeO_wtp=6.31, MgO_wtp=9.92, CaO_wtp=14.49, Na2O=4.65);
-  result <- sensors.t.grt_cpx_krogh_ravna (Grt, Cpx, 10);
+  result <- sensors.t.grt_cpx_nakamura (Grt, Cpx, 25);
-  expect_equal(result, 711, tolerance = 1);
+  expect_equal(result$T_calc_GrtCpx_Nakamura, 1036, tolerance = 1);
 }
-</code>
+</Code>