Различия

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

--- sensors:t:amp_pl_holland_blundy [2015/10/20 10:41]
admin
+++ sensors:t:amp_pl_holland_blundy [2018/11/12 15:07] (текущий)
@@ Строка 5: / Строка 5: @@
 === Details ===
-^ Article | Holland, T. , Blundy, J. Non-ideal interactions in calcic amphiboles and their bearing on amphibole-plagioclase thermometry  |
+^ Article | <ref:holland1994> |
-^ Reference | {{:sources:holland1994.pdf}} |
 ^ Collibration |  |
-^ Presition | ± 50 °C |
+^ Precision | ± 50 °C |
 ^ Ranges | 1 < P < 16 kbar,  400 °C < T < 1100 °C |
 ^ Renew sensors | Blundy and Holland (1990),  Kohn and Spear (1989, 1990).  |
-^[[:ppy:mod_sensors#t_amp_pl_holland_blundy|Calculate]] ^ ^
+<use:sensors.t.amp_pl_holland_blundy>
+<test:sensors.t.amp_pl_holland_blundy>
+<WRAP center round todo 60%>
+  * <del>заголовок</del>
+  * <del>краткое описание</del>
+  * <del>статья</del>
+  * <del>ссылка</del>
+  * карточка
+  * подробное описание
+  * <del>код</del>
+  * тест
+</WRAP>
@@ Строка 18: / Строка 29: @@
 <Code linenums lang-r>
 source('recalc_amp_holland_blundy.r');
-source('recalc_pl_8ox.r');
+source('recalc_pl_nm.r');
 source('interpmin_amp_holland_blundy.r');
@@ Строка 26: / Строка 37: @@
     Amp2 <- interpmin.amp_holland_blundy(Amp);
     Amp <- cbind(Amp1,Amp2);
-    Pl <- recalc.pl_8ox(Pl);
+    Pl <- recalc.pl_nm(Pl);
-    X_An <- Pl$Ca_pfu;
+    Pl$X_An <- Pl$Ca_pfu;
-    X_Alb <- 1 - X_An;
+    Pl$X_Ab <- 1 - Pl$X_An;
+    K_ed_tr <- (27/256)*(Amp$X_vac_A*Amp$X_Si_T1*Pl$X_Ab)/(Amp$X_Na_A*Amp$X_Al_T1);
-=((-76,95+B103*0,79+39,4*B$96+22,4*B$99+(41,5-2,89*B103)*B$94)/(-0,065-0,0083144*LN(B$100)))-273,15
+    T_A <-  (
-=((78,44+3-33,6*B$97-(66,8-2,92*B103)*B$94+78,5*B$93+9,4*B$96)/(0,0721-0,0083144*LN((27*B$97*B$92*B$39)/(64*B$98*B$93*B$38))))-273,15
+            (-76.95+P*0.79+39.4*Amp$X_Na_A+22.4*Amp$X_K_A+(41.5-2.89*P)*Amp$X_Al_M2) /
+            (-0.065-0.0083144*log(K_ed_tr))
+            )-273.15;
+    T_B <- (
+           (78.44+3-33.6*Amp$X_Na_M4-(66.8-2.92*P)*Amp$X_Al_M2+78.5*Amp$X_Al_T1+9.4*Amp$X_Na_A) /
+           (0.0721-0.0083144*
+               log(
+                   (27*Amp$X_Na_M4*Amp$X_Si_T1*Pl$X_An)/
+                   (64*Amp$X_Ca_M4*Amp$X_Al_T1*Pl$X_Ab)
+                   )
+           )
+           )-273.15;
+    return (data.frame(Name=paste0(Amp$Name,'-',Pl$Name), P_seq=P, T_calc_Holland_Blundy_A=T_A, T_calc_Holland_Blundy_B=T_B));
 }