Различия

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

--- sensors:t:opx_cpx_fonarev_graphchikov [2015/10/01 14:00]
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.
@@ Строка 9: / Строка 9: @@
 ^ Calibration |  |
 ^ Accuracy |  |
-^ Ranges | 900-1250 °C |
+^ 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>
+Multiple sensors with different constrains:
+Number 1:
+  * IS NOT TO BE USED FOR (AL2O3+TIO2+MNO+NA2O) IN CPX >5 W
+  * IS NOT TO BE USED FOR FE/(FE+MG+CA) IN CPX < 0.1
+  * IS NOT TO BE USED FOR T<900 C
+  * IS NOT TO BE USED FOR T>1250 C
+Number 2:
+  * IS NOT TO BE USED FOR T>900 C
+  * IS NOT TO BE USED FOR (AL2O3+MNO+TIO2+NA2O) IN CPX >5 W
+  * IS NOT TO BE USED FOR FE/(FE+MG+CA) IN CPX < 0.1
+Number 3:
+  * FOR MNO IN OPX  >4 W.P.
+  * IS NOT TO BE USED FOR T>900 C
+  * IS NOT TO BE USED FOR MNO IN OPX <4 W.
+  * IS NOT TO BE USED FOR FE/(FE+MG+CA) IN CPX < 0.1
+Number 4:
+  * ONLY CPX COMPOSITION IS INVOLVED IN CALCULATIONS
+  * IS NOT TO BE USED FOR FE/(FE+MG+CA) IN CPX < 0.1
+  * IS NOT TO BE USED FOR (AL2O3+MNO+TIO2+NA2O) IN CPX >5 W
+  * IS NOT TO BE USED FOR T>1250 C
+  * IS NOT TO BE USED FOR T<890 C
 === Code ===
@@ Строка 43: / Строка 79: @@
     M = (2*G*1000/(-B+sqrt(B**2-4*A*G)))-273.15;
-    T_1 = ifelse( Z>5.0 | X<0.1 | M<900.001 | M>1250.001, -1, M);
+    T_1 = ifelse( Z>5.0 | X<0.1 | M<900.001 | M>1250.001, NA, M);
     #Method2 T in 600-900
@@ Строка 54: / Строка 90: @@
     N = ifelse(N > 850 & M < 900.001,(M+N)/2,N);
-    T_2 = ifelse( Z>5.0 | X<0.1 | N<600 | N>900, -1, N);
+    T_2 = ifelse( Z>5.0 | X<0.1 | N<600 | N>900, NA, N);
+    #Method3 T in 600-900
+    X = CpxP$Fe_pfu/(CpxP$Mg_pfu+CpxP$Fe_pfu+CpxP$Ca_pfu);
+    Y = OpxP$Fe_pfu/(OpxP$Mg_pfu+OpxP$Fe_pfu+OpxP$Ca_pfu);
+    Z = Opx$MnO_wtp;
+    C = CpxP$Fe_pfu/CpxP$Mn_pfu;
+    O = OpxP$Fe_pfu/OpxP$Mn_pfu;
+    LK = log(C/O);
+    K = 1.0221-0.2816*LK;
+    A = 5.6066*X+13.6767*X^2-23.6366*X^3-60.7866*X^4-Y*(1-(0.35-X)*0.01*P);
+    B = (-9.9203*X-25.1949*X^2+30.5070*X^3+148.9344*X^4)*K;
+    G = (6.9830*X+10.0505*X^2-9.1097*X^3-85.1719*X^4)*K^2;
+    N = (2*G*1000/(-B+sqrt(B^2-4*A*G)))-273.15;
+    A = -15.9135*X+112.8452*X^2-232.7006*X^3+156.3120*X^4-Y*(1-(0.5-X)*0.01*P);
+    B = (44.2986*X-313.6564*X^2+660.4657*X^3-446.5655*X^4)*K;
+    G = (-27.1202*X+212.7500*X^2-462.5763*X^3+316.4015*X^4)*K^2;
+    M = (2*G*1000/(-B+sqrt(B**2-4*A*G)))-273.15;
+    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 = M,T_calc_OpxCpx_FonarevGraphchikov2 = 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_FonarevGraphchikov4 = T_4));
 }