Различия
Здесь показаны различия между двумя версиями данной страницы.
sensors:p:grt_opx_aranovich_kosyakova_1987_b [2016/08/12 08:04] 127.0.0.1 внешнее изменение |
sensors:p:grt_opx_aranovich_kosyakova_1987_b [2018/11/12 15:08] |
||
---|---|---|---|
Строка 1: | Строка 1: | ||
- | ====== Grt-Opx: Aranovich, Kosyakova, 1987 <B> ====== | ||
- | |||
- | === Details === | ||
- | <WRAP center round info 60%> | ||
- | This sensor from [[:about_tpf|TPF]] collection | ||
- | </WRAP> | ||
- | |||
- | ^ Article | | | ||
- | ^ Collibration | | | ||
- | ^ Presition | | | ||
- | ^ Ranges | | | ||
- | ^ Renew sensors | | | ||
- | |||
- | <use:sensors.p.grt_opx_aranovich_kosyakova_1987_b> | ||
- | <test:sensors.p.grt_opx_aranovich_kosyakova_1987_b> | ||
- | |||
- | <WRAP center round todo 60%> | ||
- | * <del>заголовок</del> | ||
- | * краткое описание | ||
- | * статья | ||
- | * ссылка | ||
- | * карточка | ||
- | * подробное описание | ||
- | * <del>код</del> | ||
- | * тест | ||
- | </WRAP> | ||
- | |||
- | === Source === | ||
- | <Code linenums lang-r> | ||
- | source("utilities_tpf.r"); | ||
- | |||
- | sensors.p.grt_opx_aranovich_kosyakova_1987_b <- function(Grt, Opx, Temp, ...){ | ||
- | T <- Temp; | ||
- | Opx <- tpf.recalc_minaral("OPX",Opx); | ||
- | Grt <- tpf.recalc_minaral("GRT",Grt); | ||
- | |||
- | A <- Grt$Al_pfu / (Grt$Al_pfu + Grt$Cr_pfu + Grt$Ti_pfu + Grt$Fe_p3_pfu); | ||
- | B <- Grt$Cr_pfu / (Grt$Al_pfu + Grt$Cr_pfu + Grt$Ti_pfu + Grt$Fe_p3_pfu); | ||
- | C <- Grt$Mg_pfu / (Grt$Mg_pfu + Grt$Fe_p2_pfu + Grt$Mn_pfu + Grt$Ca_pfu + Grt$Na_pfu + Grt$K_pfu); | ||
- | D <- Grt$Fe_p2_pfu / (Grt$Mg_pfu + Grt$Fe_p2_pfu + Grt$Mn_pfu + Grt$Ca_pfu + Grt$Na_pfu + Grt$K_pfu); | ||
- | F <- Grt$Ca_pfu / (Grt$Mg_pfu + Grt$Fe_p2_pfu + Grt$Mn_pfu + Grt$Ca_pfu + Grt$Na_pfu + Grt$K_pfu); | ||
- | X1 <- ((Opx$Al_c4_pfu + Opx$Al_c6_pfu - abs(Opx$Na_pfu - Opx$Cr_pfu) - 2 * Opx$Ti_pfu) / 2); | ||
- | X2 <- ((Opx$Al_c4_pfu + Opx$Al_c6_pfu - abs(Opx$Na_pfu - Opx$Cr_pfu) - 2 * Opx$Ti_pfu) / 2 + Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Ca_pfu); | ||
- | Y <- Opx$Mg_pfu / (Opx$Mg_pfu + Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Ca_pfu + (Opx$Al_c4_pfu + Opx$Al_c6_pfu - abs(Opx$Na_pfu - Opx$Cr_pfu) - 2 * Opx$Ti_pfu) / 2); | ||
- | Z <- Opx$Fe_p2_pfu / (Opx$Fe_p2_pfu + Opx$Fe_p3_pfu + Opx$Mg_pfu + Opx$Ca_pfu + (Opx$Al_c4_pfu + Opx$Al_c6_pfu - abs(Opx$Na_pfu - Opx$Cr_pfu) - 2 * Opx$Ti_pfu) / 2); | ||
- | X <- X1 / X2; | ||
- | TT <- T + 273.15; | ||
- | G <- log(Y^3 * X / C^3 / A^2); | ||
- | U <- log(D * Y / C / Z); | ||
- | W <- 1.987 * U + 0.49 * B - 1.688 * (Y - Z) - 1.34 * X + 2.654 + F * (3 * C + 3.436 * D + 1.466 * F + 0.658 * (D - C)); | ||
- | J <- 1237 * (3 * X + Y - 4 * X * Y) + 6441 * Z * (1 - 4 * X); | ||
- | I <- (0.3285 - 0.14245 * (3 * X^2 + (Z + Y)^2)) * 1000; | ||
- | S <- 1218 * B - 2372 * (Y - Z) - 5204 * X + F * (6600 * C + 7322 * D + 3226 * F + 1506.5 * (D - C)); | ||
- | R <- 3 * F * F * (1 - 2 * C); | ||
- | RR <- 6 * C * F * (1 - C); | ||
- | P <- (R * (3300 - 1.5 * TT) + RR * (3300 - 1.5 * TT) - 6 * F * D * (F * (74 - 0.034 * TT) + D * ( - 3661 + 1.718 * TT)) + 6 * D * F * (1 - 2 * C) * (1506.5 - 0.658 * TT) + (1 - C) * B * (29229 - 11.88 * TT) - D * B * (25575 - 10.4 * TT) - Z * (3 - 4 * Y) * ( - 2372 + 1.688 * TT) + J - 1.987 * TT * G - 2785 - 1.869 * TT - 5.88 * (TT - 970) + 5.88 * TT * log(TT / 970)) / I; | ||
- | warning(" ACCURACY +-1-8 KBAR"); | ||
- | return(data.frame(Name=paste0(Grt$Name,"-", Opx$Name),Temp=Temp,P_calc=P)); | ||
- | } | ||
- | |||
- | sensors.p.grt_opx_aranovich_kosyakova_1987_b.args <- function(){ | ||
- | ars <- list(); | ||
- | ars$minerals <- c("Grt", "Opx"); | ||
- | ars$ranges <- c("Temp"); | ||
- | return(ars); | ||
- | } | ||
- | |||
- | #sensors.p.grt_opx_aranovich_kosyakova_1987_b.test <- function(){ | ||
- | # | ||
- | #} | ||
- | </Code> | ||
- | |||
- | {{tag> tpf Grt Opx }} |