Править страницу Ссылки сюда Это старая версия документа! Bt-Opx: Sengupta, 1990 Details This sensor from TPF collection Article Collibration Precision Ranges Renew sensors Использовать Проверить заголовок краткое описание статья ссылка карточка подробное описание код тест Source source("utilities_tpf.r"); sensors.t.bt_opx_sengupta_1990 <- function(Bt, Opx, P, ...){ Name <- paste0(Bt$Name,"-", Opx$Name); Opx <- tpf.recalc_minaral("OPX",Opx); Bt <- tpf.recalc_minaral("BT",Bt); XM <- Opx$Mg_pfu / (Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu); XF <- (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu) / (Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu); KO <- (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu) / Opx$Mg_pfu; XT <- Bt$Ti_pfu / (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu + Bt$Mg_pfu + Bt$Ti_pfu + Bt$Al_c6_pfu); XA <- Bt$Al_c6_pfu / (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu + Bt$Mg_pfu + Bt$Ti_pfu + Bt$Al_c6_pfu); KB <- (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu) / Bt$Mg_pfu; K <- KO / KB; T1 <- (4130 + 603 * (XM - XF) - 4423 * XT - 3595 * XA + 0.017 * P * 1000) / (1.9872 * log(K) + 3.27) - 273; T2 <- (4130 - 4423 * XT - 3595 * XA + 0.017 * P * 1000) / (1.9872 * log(K) + 3.27) - 273; if (T2>899) { T <- T2; } if (T1<901) { T <- T1; } return(data.frame(Name=Name,P=P,T_calc=T)); } sensors.t.bt_opx_sengupta_1990.args <- function(){ ars <- list(); ars$minerals <- c("Bt", "Opx"); ars$ranges <- c("P"); return(ars); } #sensors.t.bt_opx_sengupta_1990.test <- function(){ # #} tpf, Bt, Opx sensors/t/bt_opx_sengupta_1990.1582108042.txt.gz Последние изменения: 2020/02/19 10:27 — f0ma Войти