Различия

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

--- sensors:p:opx_cpx_putirka [2015/10/08 10:41]
admin создано
+++ sensors:p:opx_cpx_putirka [2018/11/12 15:08] (текущий)
@@ Строка 1: / Строка 1: @@
-====== Opx-Cpx: Putirka, 2009 RiMG ======
+====== Opx-Cpx: Putirka, 2009 ======
-Global regression, based on the partitioning of enstatite + ferrosilite (Fm2Si2O6 = EnFs; FmO = FeO + MgO + MnO) between clinopyroxene and orthopyroxene. Cpx with Mg#>0.75.
+Global regression, based on the partitioning of enstatite + ferrosilite (Fm2Si2O6 = EnFs; FmO = FeO + MgO + MnO) between clinopyroxene and orthopyroxene.
 === Details ===
-^ Article | Keith D. Putirka Thermometers and Barometers for Volcanic Systems, Mineralogy & Geochemistry Vol. 69 pp. 61-120, 2008 |
+^ Article | <ref:putirka2008> |
-^ Reference | {{:sources:putirka2008.pdf}} |
+^ Collibration | 483 points from different sources |
-^ Collibration | 487 points with Mg# in Cpx >0.75 from large dataset |
+^ Precision | ± 3.2 kbar  |
-^ Presition | ± 46 °C  |
 ^ Ranges | P < 70 kbar,  700 °C < T < 2000 °C |
-^ Renew sensors | [[:sensors:p:opx_cpx_brey_kohler|Brey and Köhler (1990) model T(BKN)]] |
+^ Renew sensors |  |
-^[[:ppy:mod_sensors#p_opx_cpx_putirka|Calculate]] ^ ^
+<use:sensors.p.opx_cpx_putirka>
+<test:sensors.p.opx_cpx_putirka>
+<WRAP center round todo 60%>
+  * <del>заголовок</del>
+  * <del>краткое описание</del>
+  * <del>статья</del>
+  * <del>ссылка</del>
+  * <del>карточка</del>
+  * подробное описание
+  * <del>код</del>
+  * тест
+</WRAP>
 === Source ===
@@ Строка 20: / Строка 32: @@
 source('interpmin_cpx_putirka.r');
-sensors.p.opx_cpx_putirka <- function(Opx, Cpx, T=1000, usePFU = FALSE, ...){
+sensors.p.opx_cpx_putirka <- function(Opx, Cpx, Temp=1000, usePFU = FALSE, ...){
 stopifnot(usePFU == F);
@@ Строка 30: / Строка 42: @@
 Opx <- recalc.px_putirka(Opx);
-Pcalc <- 10^4/(11.2-1.96*log(CpxCmp$X_EnFs/OpxCmp$X_Fm2Si2O6)-3.3*Cpx$Ca_pfu-25.8*CpxCmp$X_CrCaTs+33.2*Opx$Mn_pfu-23.6*Opx$Na_pfu-2.08*OpxCmp$X_En-8.33*OpxCmp$X_Di-0.05*P);
+Kf <- Opx$Ca_pfu/(1-Cpx$Ca_pfu);
+Opx_En_Act <- (0.5*Opx$Mg_pfu/(Opx$Ca_pfu+0.5*Opx$Mg_pfu+0.5*(Opx$Fe_pfu-Opx$Fe_p3_pfu)+Opx$Mn_pfu+Opx$Na_pfu))*
+              (0.5*Opx$Mg_pfu/(Opx$Fe_p3_pfu+0.5*Opx$Fe_p2_pfu+Opx$Al_c6_pfu+Opx$Ti_pfu+Opx$Cr_pfu+0.5*Opx$Mg_pfu));
+Cpx_Di_Act <- Cpx$Ca_pfu/(Cpx$Ca_pfu+0.5*Cpx$Mg_pfu+0.5*Cpx$Fe_p2_pfu+Cpx$Mn_pfu+Cpx$Na_pfu);
+Pcalc <- -94.25 +0.045*Temp +187.7*Opx$Al_c6_pfu +246.8*OpxCmp$X_Fm2Si2O6 -212.5*OpxCmp$X_En +
+.5*Opx_En_Act -69.4*CpxCmp$X_EnFs -133.9*Cpx_Di_Act -1.66/Kf;
 Eq_check <- (Cpx$Fe_pfu/Cpx$Mg_pfu)/(Opx$Fe_pfu/Opx$Mg_pfu);
@@ Строка 36: / Строка 56: @@
 Eq_check <- Eq_check > 0.95 && Eq_check < 1.23;
-return (data.frame(Name=paste0(Opx$Name,'-',Cpx$Name),T_seq=T,P_calc_Opx_Cpx_Putirka=Pcalc,Eq_check=Eq_check));
+return (data.frame(Name=paste0(Opx$Name,'-',Cpx$Name),T_seq=Temp,P_calc_Opx_Cpx_Putirka=Pcalc,Eq_check=Eq_check));
 }
-sensors.t.opx_cpx_putirka.args <- function() {
+sensors.p.opx_cpx_putirka.args <- function() {
     ars <- list();
     ars$minerals = c('Opx','Cpx');
-    ars$ranges = c('T');
+    ars$ranges = c('Temp');
     return(ars);
 }
 </Code>
+{{tag> AP_igroup fixed}}