====== Grt-Bt: Indares, Martignole, 1985 ====== === Details === This sensor from [[:about_tpf|TPF]] collection ^ Article | | ^ Collibration | | ^ Precision | | ^ Ranges | | ^ Renew sensors | | * заголовок * краткое описание * статья * ссылка * карточка * подробное описание * код * тест === Source === source("utilities_tpf.r"); sensors.t.grt_bt_indares_martignole_1985_a <- function(Grt, Bt, P, ...){ Bt <- tpf.recalc_minaral("BT1",Bt); Grt <- tpf.recalc_minaral("GRT",Grt); Y <- (Grt$Fe_p2_pfu + Grt$Fe_p3_pfu) / Grt$Mg_pfu; Z <- Grt$Ca_pfu / (Grt$Mg_pfu + Grt$Ca_pfu + Grt$Mn_pfu + Grt$Fe_p2_pfu + Grt$Fe_p3_pfu); N <- (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu) / Bt$Mg_pfu; X <- Bt$Al_c6_pfu / (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu + Bt$Mg_pfu + Bt$Mn_pfu + Bt$Ti_pfu + Bt$Al_c6_pfu); M <- Bt$Ti_pfu / (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu + Bt$Mg_pfu + Bt$Mn_pfu + Bt$Ti_pfu + Bt$Al_c6_pfu); K <- log(N / Y); T <- (12454 + 0.057 * P * 1000 + 3 * ( - 454 * X - 6767 * M) + 9900 * Z) / (4.662 - 5.9616 * K + 4.5 * Z) - 273; return(data.frame(Name=paste0(Grt$Name,"-", Bt$Name),P=P,T_calc=T)); } sensors.t.grt_bt_indares_martignole_1985_a.args <- function(){ ars <- list(); ars$minerals <- c("Grt", "Bt"); ars$ranges <- c("P"); return(ars); } #sensors.t.grt_bt_indares_martignole_1985_a.test <- function(){ # #} {{tag> tpf Grt Bt }}