====== Opx: Opx components, Thompson, 1982 ======

Oxide end members, Thompson’s (1982) : Fm6O6, Ca2Si2O6, Na4Si2O6, Ti3O6, Si3O6, Al4O6 and Cr4O6.

=== Details ===

^ Article |  |
^ Recalc | [[recalc:px_putirka]] |
^ Results | En, Di, NaAlSi2O6, FmTiAlSiO6, CrAl2SiO6, FmAl2SiO6, CaFmSi2O6, Fm2Si2O6 |

<use:interpmin.opx_thompson>
<test:interpmin.opx_thompson>

==== Код ====
<Code linenums lang-r>
source('recalc_px_putirka.r');

interpmin.opx_thompson <- function (wtp_data, ...)
{
p <- recalc.px_putirka(wtp_data);
comps <- data.frame(Name=p$Name);

comps$X_NaAlSi2O6 <- p$Na_pfu;
comps$X_FmTiAlSiO6 <- p$Ti_pfu;
comps$X_CrAl2SiO6 <- p$Cr_pfu;
comps$X_FmAl2SiO6 <- p$Al_c6_pfu - comps$X_NaAlSi2O6 - comps$X_CrAl2SiO6;
comps$X_CaFmSi2O6 <- p$Ca_pfu;
comps$X_Fm2Si2O6 <- ((p$Fe_pfu+p$Mg_pfu+p$Mn_pfu)-comps$X_FmTiAlSiO6-comps$X_FmAl2SiO6-comps$X_CaFmSi2O6)/2;

comps$X_En <- comps$X_Fm2Si2O6*(p$Mg_pfu/(p$Mg_pfu+p$Mn_pfu+p$Fe_pfu));
comps$X_Di <- comps$X_CaFmSi2O6*(p$Mg_pfu/(p$Mg_pfu+p$Mn_pfu+p$Fe_pfu));

return(comps);
}

interpmin.opx_thompson.test <- function (){
warning("No test routine!");
}
</Code>