Различия

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

--- spdgrams:ree_primitive_mantle [2016/03/11 09:54]
admin создано
+++ spdgrams:ree_primitive_mantle [2018/11/12 15:07] (текущий)
@@ Строка 1: / Строка 1: @@
-====== Primitive Mantle ======
+====== Primitive Mantle, Sun-McDonough, 1989 ======
 === Details ===
-<data plot>
+| Reference | <ref:sun1989> |
-Reference [wiki]:
+| Data | bulk |
-Data: bulk
+| Input | ppm |
-Input: ppm
-</data>
+<use:spdgrams.ree_primitive_mantle>
+<test:spdgrams.ree_primitive_mantle>
 <WRAP center round todo 60%>
-  * заголовок
+  * <del>заголовок</del>
   * краткое описание
-  * статья
+  * <del>статья</del>
-  * ссылка
+  * <del>ссылка</del>
   * <del>карточка</del>
   * подробное описание
-  * код
+  * <del>код</del>
   * тест
 </WRAP>
@@ Строка 25: / Строка 27: @@
 <Code linenums lang-r>
-spdgrams.ree_primitive_mantle <- function(ppm_data)
+require('ggplot2');
+cmp <- read.csv("constants_oceanic_basalts_compositions.csv", sep=";", stringsAsFactors=F);
+spdgrams.ree_primitive_mantle <- function(ppm_data, el_order = NULL)
 {
-primitive_mantle = data.frame(Cs = c(0.032), Rb = c(0.635), Ba = c(6.989), Th = c(0.085),
+ns <- el_order;
-                              U = c(0.021), Nb = c(0.713), Ta = c(0.041), La = c(0.687),
-                              Ce = c(1.775), Pb = c(0.185), Pr = c(0.276), Sr = c(21.1),
+if (is.null(el_order))
-                              Nd = c(1.354), Zr = c(11.2), Hf = c(0.309), Sm = c(0.444),
+{
-                              Eu = c(0.168), Gd = c(0.596), Tb = c(0.108), Dy = c(0.737),
+ns <- c('Cs', 'Rb', 'Ba', 'Th', 'U', 'Nb', 'Ta', 'La', 'Ce', 'Pb', 'Pr',
-                              Y = c(4.55), Ho = c(0.164), Er = c(0.48), Tm = c(0.074),
+        'Sr', 'Nd', 'Zr', 'Hf', 'Sm', 'Eu', 'Gd', 'Tb', 'Dy', 'Y', 'Ho',
-                              Yb = c(0.493), Lu = c(0.074));
+        'Er', 'Tm', 'Yb', 'Lu');
+}
+pmsc <- c();
-ns <- names(primitive_mantle);
+for (n in ns)
+  pmsc <- c(pmsc, subset(cmp, element == n)$primitive_mantle);
 ds <- selectNames(ppm_data,c(paste0(ns,'_ppm')));
+names(ds) <- ns;
 ret <- c();
-ret$values <- rowApply(`/`,ds,primitive_mantle[1,]);
+values <- cbind(ppm_data$Name, rowApply(`/`,ds,pmsc));
-ret$names <- ppm_data$Name;
+names(values)[1] <- 'Name';
+ret$values <- values;
 ret$grades <- ns;
+ret$labels <- c('REE', 'Sample/Primitive mantle');
 return(ret);
 }
 </Code>
+{{tag>spdgrams mantle}}