Различия
Здесь показаны различия между двумя версиями данной страницы.
Следующая версия | Предыдущая версия | ||
sensors:t:hbl_pl_spear_1981_b [2016/08/12 07:41] 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 | | | ||
Строка 29: | Строка 29: | ||
=== Source === | === Source === | ||
<Code linenums lang-r> | <Code linenums lang-r> | ||
- | 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(Pl$Name,'-', Hbl$Name),P=P,T_calc=T)); | + | return(data.frame(Name=paste0(Pl$Name,"-", Hbl$Name),P=P,T_calc=T)); |
} | } | ||
sensors.t.hbl_pl_spear_1981_b.args <- function(){ | sensors.t.hbl_pl_spear_1981_b.args <- function(){ | ||
ars <- list(); | ars <- list(); | ||
- | ars$minerals = c('Hbl', 'Pl'); | + | ars$minerals <- c("Hbl", "Pl"); |
- | ars$ranges = c('P'); | + | ars$ranges <- c("P"); |
return(ars); | return(ars); | ||
} | } | ||
- | #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 }} |