====== Pl-L: Kudo, Weill, 1970 ====== === Details === This sensor from [[:about_tpf|TPF]] collection ^ Article | | ^ Collibration | | ^ Precision | | ^ Ranges | | ^ Renew sensors | | * заголовок * краткое описание * статья * ссылка * карточка * подробное описание * код * тест === Source === source("utilities_tpf.r"); sensors.t.pl_l_kudo_weill_1970 <- function(Pl, L, P, ...){ Pl <- tpf.recalc_minaral("PL",Pl); L <- tpf.recalc_minaral("L",L); X <- Pl$Na_pfu / Pl$Ca_pfu; K <- L$Ca_pfu / (L$Ca_pfu + L$Al_c4_pfu + L$Al_c6_pfu + L$Na_pfu + L$Si_pfu); L <- (L$Al_c4_pfu + L$Al_c6_pfu) / (L$Ca_pfu + L$Al_c4_pfu + L$Al_c6_pfu + L$Na_pfu + L$Si_pfu); M <- L$Na_pfu / (L$Ca_pfu + L$Al_c4_pfu + L$Al_c6_pfu + L$Na_pfu + L$Si_pfu); N <- L$Si_pfu / (L$Ca_pfu + L$Al_c4_pfu + L$Al_c6_pfu + L$Na_pfu + L$Si_pfu); F <- K + L - M - N; R <- M * N / K / L; if (P<0.005) { Q <- 17.24 + log(R / X); T <- (Q + sqrt(Q^2 + 4 * 0.01034 * 12900 * F)) / 2 / 0.01034 - 273; } if (P>0.4) { if (P<0.6) { Q <- 17.86 + log(R / X); T <- (Q + sqrt(Q^2 + 4 * 0.01105 * 12900 * F)) / 2 / 0.01105 - 273; } } if (P>0.9) { if (P<1.1) { Q <- 17.67 + log(R / X); T <- (Q + sqrt(Q^2 + 4 * 0.01114 * 12900 * F)) / 2 / 0.01114 - 273; } } if (P>4.9) { if (P<5.1) { Q <- 16.63 + log(R / X); T <- (Q + sqrt(Q^2 + 4 * 0.01218 * 12900 * F)) / 2 / 0.01218 - 273; } } if (P>5.1) { warning(" IS NOT TO BE USED FOR P>5.1 KBAR"); } return(data.frame(Name=paste0(L$Name,"-", Pl$Name),P=P,T_calc=T)); } sensors.t.pl_l_kudo_weill_1970.args <- function(){ ars <- list(); ars$minerals <- c("Pl", "L"); ars$ranges <- c("P"); return(ars); } #sensors.t.pl_l_kudo_weill_1970.test <- function(){ # #} {{tag> tpf L Pl }}