Различия
Здесь показаны различия между двумя версиями данной страницы.
| Предыдущая версия справа и слева Предыдущая версия Следующая версия | Предыдущая версия | ||
|
sensors:t:hbl_pl_spear_1981_b [2016/08/12 07:45] 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; | ||
| - | K <- Hbl$Mg_pfu+Hbl$Fe_p2_pfu+Hbl$Fe_p3_pfu+Hbl$Ti_pfu+Hbl$Cr_pfu+Hbl$Al_c4_pfu+Hbl$Al_c6_pfu+Hbl$Si_pfu; | + | K <- Hbl$Mg_pfu + Hbl$Fe_p2_pfu + Hbl$Fe_p3_pfu + Hbl$Ti_pfu + Hbl$Cr_pfu + Hbl$Al_c4_pfu + Hbl$Al_c6_pfu + Hbl$Si_pfu; |
| - | N <- 13/K; | + | N <- 13 / K; |
| - | 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 Pl Hbl }} | ||