Различия

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

--- sensors:t:amp_pl_holland_blundy [2015/10/20 10:24]
admin created
+++ sensors:t:amp_pl_holland_blundy [2018/11/12 15:07] (текущий)
@@ Строка 1: / Строка 1: @@
 ====== Amp-Pl: Holland, Blundy, 1994 ======
+For edenite + 4 quartz = tremolite + albite (A) and edenite + albite = richterite + anorthite (B) equilibrium.
 === 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 | ±  °C |
+^ Precision | ± 50 °C |
-^ Ranges | < P <   kbar,   °C < T <  °C |
+^ Ranges | 1 < P < 16 kbar,  400 °C < T < 1100 °C |
-^ Renew sensors |  |
+^ 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>
@@ Строка 17: / Строка 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');
+sensors.t.amp_pl_holland_blundy <- function(Amp, Pl, P, usePFU = FALSE, ...){
+    Amp1 <- recalc.amp_holland_blundy(Amp);
+    Amp2 <- interpmin.amp_holland_blundy(Amp);
+    Amp <- cbind(Amp1,Amp2);
+    Pl <- recalc.pl_nm(Pl);
+    Pl$X_An <- Pl$Ca_pfu;
+    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);
+    T_A <-  (
+            (-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));
+}
+sensors.t.amp_pl_holland_blundy.args <- function() {
+    ars <- list();
+    ars$minerals = c('Amp','Pl');
+    ars$ranges = c('P');
+    return(ars);
+}
 </Code>