Различия

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

--- sensors:t:electrum_sp_shikazono [2016/02/09 14:14]
admin created
+++ sensors:t:electrum_sp_shikazono [2018/11/12 15:08] (текущий)
@@ Строка 5: / Строка 5: @@
 === Details ===
-^ Article | Shikazono N. A comparison of temperatures estimated from the electrum – sphalerite – pyrite – argentite assemblage and filling temperatures of fluid implications from epithermal Au – Ag  vein –type deposits in Japan  Econ. Geol. 1985. V. 80. No 5. P. 1415 – 1424.   |
+^ Article | <ref:shikazono1985> |
-^ Reference | {{:sources:shikazono1985.pdf}} |
 ^ Collibration | In this study, sphalerite and electrum from Sanru, Chitose, Sado, Takatama, Seigoshi, Yugashima, Toi, Nebazawa, Takeno, and Kushikino were analyzed by an electron microprobe analyzer  |
+^ Ranges | N_Ag > 0.3 |
+^ Input | atomic wtp |
 ^ Accuracy | T ±20 °C (for T > 300 °C)  |
-^[[:ppy:mod_sensors#t_elect_spf_shikazono|Calculate]] ^ ^
+<use:sensors.t.electrum_sp_shikazono>
+<test:sensors.t.electrum_sp_shikazono>
+<WRAP center round todo 60%>
+Для улучшения этой статьи следует:
+  * Добавить подробное описание
+</WRAP>
@@ Строка 15: / Строка 22: @@
 <Code linenums lang-r>
-sensors.t.elect_spf_shikazono <- function(Elect, Spf)
+source('recalc_wtp_at_to_atp.r')
+sensors.t.electrum_sp_shikazono <- function(Elect, Spf)
 {
-    dElect <- selectNames(Elect,c('Au_wtp','Ag_wtp')) / 100;
+    aElect <- recalc.wtp_at_to_atp(Elect);
-    dSpf <- selectNames(Spf,c('Fe_wtp')) / 100;
+    aSpf <- recalc.wtp_at_to_atp(Spf);
-    dElect <- rowApply('/',dElect,c(196.96657,107.8682));
-    names(dElect) <- c('Au_moles','Ag_moles');
-    dSpf <- dSpf / 55.845;
-    N_Ag = dElect$Ag_moles/(dElect$Ag_moles+dElect$Au_moles);
+    dElect <- selectNames(aElect,c('Au_atp','Ag_atp'));
-    X_FeS = dSpf;
+    dSpf <- selectNames(aSpf,c('Fe_atp','Zn_atp'));
+    N_Ag = dElect$Ag_atp/(dElect$Ag_atp+dElect$Au_atp);
+    X_FeS = dSpf$Fe_atp/(dSpf$Fe_atp+dSpf$Zn_atp);
     N = 28765 + 22600*(1 - N_Ag)^2 - 6400*(1 - N_Ag)^3;
@@ Строка 41: / Строка 49: @@
               (4*(1 - N_Ag)^2*(5650 - 1600*(1 - N_Ag) - 1.375*Temp)));
-    return(data.frame(Name=paste0(Elect$Name,'-',Spf$Name), N_Ag=N_Ag,X_FeS =X_FeS, T_calc_ElectSpf_Shikazono_eq6 = Tcalc6, T_calc_ElectSpf_Shikazono_eq7 = Tcalc7,log_f_S=log_f_S));
+    return(data.frame(Name=paste0(Elect$Name,'-',Spf$Name), N_Ag=N_Ag, X_FeS =X_FeS,
+           T_calc_ElectSp_Shikazono_eq6 = Tcalc6,
+           T_calc_ElectSp_Shikazono_eq7 = Tcalc7,
+           log_f_S=log_f_S));
 }
-sensors.t.elect_spf_shikazono.args <- function() {
+sensors.t.electrum_sp_shikazono.args <- function() {
     ars <- list();
-    ars$minerals = c('Elect','Spf');
+    ars$minerals = c('Elect','Sp');
     ars$ranges = c();
     return(ars);
 }
-sensors.t.elect_spf_shikazono.test <- function () {
+sensors.t.electrum_sp_shikazono.test <- function () {
   Elect <- data.frame(Name='Elect',Au_wtp=2, Ag_wtp=8);
-  Spf <- data.frame(Name='Spf',Fe_wtp=21.23);
+  Spf <- data.frame(Name='Sp',Fe_wtp=21.23,Zn_wtp=21.23);
-  result <- sensors.t.elect_spf_shikazono(Elect, Spf);
+  result <- sensors.t.electrum_sp_shikazono(Elect, Spf);
-  expect_equal(result$T_calc_ElectSpf_Shikazono_eq6, 312, tolerance = 3);
+  expect_equal(result$T_calc_ElectSp_Shikazono_eq6, 312, tolerance = 3);
-  expect_equal(result$T_calc_ElectSpf_Shikazono_eq7, 297, tolerance = 3);
+  expect_equal(result$T_calc_ElectSp_Shikazono_eq7, 297, tolerance = 3);
 }
 </Code>