Различия

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

--- sensors:p:hbl_thomas_ernst_1990 [2017/01/17 13:44]
127.0.0.1 внешнее изменение
+++ sensors:p:hbl_thomas_ernst_1990 [2018/11/12 15:08] (текущий)
@@ Строка 1: / Строка 1: @@
 ====== Hbl: Thomas, Ernst, 1990 ======
+ONLY FOR HBL+BT+EP+SPL+MELT+KFS+PL+QTZ+FE-TI OXIDE
 === Details ===
@@ Строка 7: / Строка 9: @@
 </WRAP>
-^ Article |  |
+^ Article | <ref:thomas1990> |
-^ Collibration |  |
+^ Callibration | 10 experements |
-^ Precision |  |
+^ Precision | ±1.0 kbar |
-^ Ranges |  |
+^ Ranges | 6-12 kbar |
-^ Renew sensors |  |
+^ Renew sensors | HAMMARSTROM and ZEN ( 1987), HOLLISTER and OTHERS (1988);  |
 <use:sensors.p.hbl_thomas_ernst_1990>
@@ Строка 18: / Строка 20: @@
 <WRAP center round todo 60%>
   * <del>заголовок</del>
-  * краткое описание
+  * <del>краткое описание</del>
-  * статья
+  * <del>статья</del>
-  * ссылка
+  * <del>ссылка</del>
-  * карточка
+  * <del>карточка</del>
   * подробное описание
   * <del>код</del>
-  * тест
+  * <del>тест</del>
 </WRAP>
@@ Строка 31: / Строка 33: @@
 source("utilities_tpf.r");
-sensors.p.hbl_thomas_ernst_1990 <- function(Hbl, ...){
+sensors.p.hbl_thomas_ernst_1990 <- function(Hbl, usePFU=FALSE, ...){
-    Hbl <- tpf.recalc_minaral("HBL",Hbl);
+    row_names = Hbl$Name;
-    XAL <- Hbl$Al_c4_pfu + Hbl$Al_c6_pfu;
+    if (!usePFU) {
-    PP <-  - 6.23 + 5.34 * XAL;
+        Hbl <- tpf.recalc_minaral("HBL",Hbl);
-    if (PP<12.00) {
-    P <- PP;
-    warning(" ONLY FOR HBL+BT+EP+SPL+MELT+KFS+PL+QTZ+FE-TI OXIDE");
     }
-    if (PP>12.00) {
-    warning(" CANNOT BE USED FOR THESE COMPOSITIONS");
+    XAL <- Hbl$Al_c4_pfu + Hbl$Al_c6_pfu;
-    }
+    PP <-  -6.23 + 5.34 * XAL;
-    return(data.frame(Name=paste0(Hbl$Name),P_calc=P));
+    P <- ifelse(PP<12,PP,NA);
+    return(data.frame(Name=row_names,P_calc=P));
 }
@@ Строка 53: / Строка 56: @@
 }
-#sensors.p.hbl_thomas_ernst_1990.test <- function(){
+sensors.p.hbl_thomas_ernst_1990.test <- function(){
-#
+    hbl <- data.frame(Name=c('HBL'), SiO2_wtp=c(39.44), TiO2_wtp=c(1.16), Al2O3_wtp=c(15.02), Cr2O3_wtp=c(0.00),
-#}
+                        FeO_wtp=c(18.11), MnO_wtp=c(0.36), MgO_wtp=c(7.90), CaO_wtp=c(11.45), Na2O_wtp=c(1.92), K2O_wtp=c(1.58));
+    res <- sensors.p.hbl_thomas_ernst_1990(hbl);
+    expect_equal(object = res[1,]$P_calc, 8, tolerance = 1);
+}
 </Code>
-{{tag> tpf Hbl }}
+{{tag> tpf Hbl AP_igroup fixed}}