Различия
Здесь показаны различия между двумя версиями данной страницы.
Предыдущая версия справа и слева Предыдущая версия Следующая версия | Предыдущая версия | ||
sensors:t:hbl_pl_spear_1981_b [2016/08/12 07:49] admin |
sensors:t:hbl_pl_spear_1981_b [2018/11/12 15:07] (текущий) |
||
---|---|---|---|
Строка 9: | Строка 9: | ||
^ Article | | | ^ Article | | | ||
^ Collibration | | | ^ Collibration | | | ||
- | ^ Presition | | | + | ^ Precision | | |
^ Ranges | | | ^ Ranges | | | ||
^ Renew sensors | | | ^ Renew sensors | | | ||
Строка 31: | Строка 31: | ||
source("utilities_tpf.r"); | source("utilities_tpf.r"); | ||
- | #main code | ||
sensors.t.hbl_pl_spear_1981_b <- function(Hbl, Pl, P, ...){ | sensors.t.hbl_pl_spear_1981_b <- function(Hbl, Pl, P, ...){ | ||
Pl <- tpf.recalc_minaral("PL",Pl); | Pl <- tpf.recalc_minaral("PL",Pl); | ||
Hbl <- tpf.recalc_minaral("HBL",Hbl); | Hbl <- tpf.recalc_minaral("HBL",Hbl); | ||
- | X <- Pl$Ca_pfu / Pl$Na_pfu; | + | X <- Pl$Na_pfu / (Pl$Na_pfu + Pl$K_pfu + Pl$Ca_pfu); |
Y <- Hbl$Ca_pfu; | Y <- Hbl$Ca_pfu; | ||
Z <- Hbl$Na_pfu; | Z <- Hbl$Na_pfu; | ||
Строка 43: | Строка 42: | ||
Y <- Y * N; | Y <- Y * N; | ||
Z <- Z * N; | Z <- Z * N; | ||
- | if (Y<2.0) { | + | if (Y<2) { |
if (Y + Z>2) { | if (Y + Z>2) { | ||
- | N <- 2 - Y; | + | N <- Z - 2 + Y; |
- | T <- (0.24 * P * 1000 - 47000) / (2 * 1.9872 * log(X * N / Y) - 44.27) - 273; | + | T <- ( - 0.257 * P * 1000 + 7075) / ( - 1.9872 * log(N / (1 - N) / X) + 6.86) - 273; |
} | } | ||
if (Y + Z<2) { | if (Y + Z<2) { | ||
- | T <- (0.24 * P * 1000 - 47000) / (2 * 1.9872 * log(X * Z / Y) - 44.27) - 273; | + | warning(" CANNOT BE USED FOR (CA+NA) IN HBL < 2"); |
} | } | ||
} | } | ||
- | if (Y>2.0) { | + | if (Y>2) { |
warning(" CANNOT BE USED FOR CA IN HBL > 2"); | warning(" CANNOT BE USED FOR CA IN HBL > 2"); | ||
} | } | ||
- | return(data.frame(Name=paste0(Hbl$Name,"-", Pl$Name),P=P,T_calc=T)); | + | return(data.frame(Name=paste0(Pl$Name,"-", Hbl$Name),P=P,T_calc=T)); |
} | } | ||
Строка 65: | Строка 64: | ||
} | } | ||
- | #test stub | ||
#sensors.t.hbl_pl_spear_1981_b.test <- function(){ | #sensors.t.hbl_pl_spear_1981_b.test <- function(){ | ||
+ | # | ||
#} | #} | ||
</Code> | </Code> | ||
- | {{tag> tpf Hbl Pl }} | + | {{tag> tpf Pl Hbl }} |