Различия

Здесь показаны различия между двумя версиями данной страницы.

--- sensors:p:grt_cpx_beyer [2017/06/02 10:14]
admin
+++ sensors:p:grt_cpx_beyer [2018/11/12 15:08] (текущий)
@@ Строка 34: / Строка 34: @@
 sensors.p.grt_cpx_beyer <- function(Grt, Cpx, Temp, usePFU = FALSE, ...){
-    row_names <- paste0(Grt$Name,"-",Cpx$Name);
+    presults <- data.frame(Name = paste0(Grt$Name,"-",Cpx$Name));
     if(!usePFU){
         Grt <- recalc.grt_nm(Grt);
-        Cpx <- recalc.px_nm(Grt);
+        Cpx <- recalc.px_nm(Cpx);
     }
@@ Строка 98: / Строка 98: @@
     Grt$Pyr_RTLn_X3Y2 <- 8.314*T*log(Grt$X_Mg^3);
-    Grt$Grs_RTlna_MgAl <-  3*Grt$Pyr_RTg + Grt$Pyr_RTLn_X3Y2;
+    Grt$Pyr_RTlna_MgAl <-  3*Grt$Pyr_RTg + Grt$Pyr_RTLn_X3Y2;
     Cpx$W_MgFe <- 0;
@@ Строка 145: / Строка 145: @@
                              Cpx$CaTs_RTLn_tet + Cpx$CaTs_Recip;
-    return(Cpx);
+    presults$RTLnK <- Cpx$CaTs_RTLna_CaAlAl +Cpx$Dio_RTLna_CaMgSi -1/3*Grt$Pyr_RTlna_MgAl -2/3*Grt$Grs_RTlna_CaAl;
+    presults$K <- exp(presults$RTLnK/(8.314*T));
+    presults$dGt <- 16500-23.5*T;
+    presults$P <- -(presults$dGt+8.314*T*log(presults$K))/(0.719*1000);
+    # if NA -> "no Al4 in cpx!"
+    return(data.frame(Name = presults$Name, T=Temp, P_calc = presults$P));
 }
@@ Строка 165: / Строка 174: @@
 </Code>
+{{tag> MA_igroup Grt Cpx fixed}}