Различия

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

--- sensors:t:opx_cpx_fonarev_graphchikov [2015/10/01 15:05]
admin
+++ sensors:t:opx_cpx_fonarev_graphchikov [2018/11/12 15:08] (текущий)
@@ Строка 1: / Строка 1: @@
-====== Opx-Cpx termometer Fonarev-Graphchikov, 1991 ======
+====== Opx-Cpx: Fonarev, Graphchikov, 1991 ======
 Two-pyroxene thermometry: a critical evaluation.
@@ Строка 11: / Строка 11: @@
 ^ Ranges | 600-1250 °C |
 ^ Renew sensors |  |
-^[[:ppy:mod_sensors#t_opx_cpx_fonarev_graphchikov|Рассчитать]] ^ ^
+<use:sensors.t.opx_cpx_fonarev_graphchikov>
+<test:sensors.t.opx_cpx_fonarev_graphchikov>
+<WRAP center round todo 60%>
+  * <del>заголовок</del>
+  * краткое описание
+  * статья
+  * ссылка
+  * карточка
+  * подробное описание
+  * <del>код</del>
+  * тест
+</WRAP>
@@ Строка 86: / Строка 98: @@
     Y = OpxP$Fe_pfu/(OpxP$Mg_pfu+OpxP$Fe_pfu+OpxP$Ca_pfu);
     Z = Opx$MnO_wtp;
-    CpxP$Mn_pfu = CpxP$Mn_pfu + 0.0001;
-    OpxP$Mn_pfu = OpxP$Mn_pfu + 0.0001;
     C = CpxP$Fe_pfu/CpxP$Mn_pfu;
     O = OpxP$Fe_pfu/OpxP$Mn_pfu;
@@ Строка 107: / Строка 117: @@
     T_3 = ifelse( Z<4.0 | X<0.10 | N<600.001 | M>900.001, NA, N);
     T_3 = ifelse(T_3<850, N, (M+N)/2);
+    #Method 4
+    X = CpxP$Fe_pfu/(CpxP$Mg_pfu+CpxP$Fe_pfu+CpxP$Ca_pfu);
+    Y = CpxP$Ca_pfu/(CpxP$Mg_pfu+CpxP$Fe_pfu+CpxP$Ca_pfu);
+    Z = Cpx$Al2O3_wtp+Cpx$TiO2_wtp+Cpx$MnO_wtp+Cpx$Na2O_wtp;
+    A = 0.0435-2.0422*X+14.2661*X^2-38.4323*X^3+34.4206*X^4-Y*(1-(0.5-X)*0.012*P);
+    B = 0.6859+1.2269*X-37.7330*X^2+119.7173*X^3-109.3990*X^4;
+    G = -0.2832+0.9887*X+27.5921*X^2-97.3991*X^3+88.8668*X^4;
+    Q = (2*G*1000/(-B+sqrt(B^2-4*G*A)))-273.15;
+    T_4 = ifelse( Z>5.0 | X<0.10 | Q<890 | Q>1250, NA, N);
     return (data.frame(Name=paste0(Opx$Name,'-',Cpx$Name),P = P,
                        T_calc_OpxCpx_FonarevGraphchikov1 = T_1,
                        T_calc_OpxCpx_FonarevGraphchikov2 = T_2,
-                       T_calc_OpxCpx_FonarevGraphchikov3 = Y));
+                       T_calc_OpxCpx_FonarevGraphchikov4 = T_4));
 }