Различия
Здесь показаны различия между двумя версиями данной страницы.
Следующая версия | Предыдущая версия | ||
sensors:p:cpx_pl_qtz_raith_1983 [2018/11/12 15:08] 127.0.0.1 внешнее изменение |
sensors:p:cpx_pl_qtz_raith_1983 [2021/09/23 15:21] f0ma |
||
---|---|---|---|
Строка 31: | Строка 31: | ||
source("utilities_tpf.r"); | source("utilities_tpf.r"); | ||
- | sensors.p.cpx_pl_qtz_raith_1983 <- function(Cpx, Pl, Qtz, Temp, ...){ | + | sensors.p.cpx_pl_qtz_raith_1983 <- function(Cpx, Pl, Temp, ...){ |
- | T <- Temp; | + | Name <-paste0(Cpx$Name,"-", Pl$Name); |
+ | T <- Temp; | ||
Pl <- tpf.recalc_minaral("PL",Pl); | Pl <- tpf.recalc_minaral("PL",Pl); | ||
Cpx <- tpf.recalc_minaral("CPX",Cpx); | Cpx <- tpf.recalc_minaral("CPX",Cpx); | ||
Строка 59: | Строка 60: | ||
L <- (1 - Y)^2 * (2075 + 9318 * Y); | L <- (1 - Y)^2 * (2075 + 9318 * Y); | ||
M <- 26885 * (Z * (1 - Z))^2; | M <- 26885 * (Z * (1 - Z))^2; | ||
+ | P <- NA; | ||
if (Z>0.01999) { | if (Z>0.01999) { | ||
P <- (5360 + K + 12864 * Z * (1 - Z) - M - L) / 0.349 / 1000; | P <- (5360 + K + 12864 * Z * (1 - Z) - M - L) / 0.349 / 1000; | ||
Строка 66: | Строка 68: | ||
warning(" CANNOT BE USED FOR (AL-CR-2*TI-FE3+NA)/2 IN CPX <0.02"); | warning(" CANNOT BE USED FOR (AL-CR-2*TI-FE3+NA)/2 IN CPX <0.02"); | ||
} | } | ||
- | return(data.frame(Name=paste0(Pl$Name,"-", Cpx$Name),Temp=Temp,P_calc=P)); | + | return(data.frame(Name=Name,Temp=Temp,P_calc=P,Z_calc=Z)); |
} | } | ||
sensors.p.cpx_pl_qtz_raith_1983.args <- function(){ | sensors.p.cpx_pl_qtz_raith_1983.args <- function(){ | ||
ars <- list(); | ars <- list(); | ||
- | ars$minerals <- c("Cpx", "Pl", "Qtz"); | + | ars$minerals <- c("Cpx", "Pl"); |
ars$ranges <- c("Temp"); | ars$ranges <- c("Temp"); | ||
return(ars); | return(ars); |