Различия

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

--- classifydiag:la_y_nb_cabanis [2016/03/04 14:49]
admin создано
+++ classifydiag:la_y_nb_cabanis [2018/11/12 15:07] (текущий)
@@ Строка 1: / Строка 1: @@
-====== La-Y-Nb: Continental crust magmatic mixtures, Cabanis ======
+====== La-Y-Nb: Continental crust magmatic mixtures, Cabanis, 1989 ======
-<ref:cabanis1989>
+Continental crust magmatic mixtures. NMORB - Normal mid-ocean-rige basalt, EMORB - Enriched mid-ocean-rige basalt, VAC - volcanic arc tholeiites, * - back arc basin, Cont. - continental.
+==== Details ====
+| Reference | <ref:cabanis1989> |
+| Data | bulk |
+| Input | La_ppm, Y_ppm, Nb_ppm |
+<use:classifydiag.la_y_nb_cabanis>
+<test:classifydiag.la_y_nb_cabanis>
 <WRAP center round todo 60%>
-  * заголовок
+  * <del>заголовок</del>
   * краткое описание
   * статья
-  * ссылка
+  * <del>ссылка</del>
-  * карточка
+  * <del>карточка</del>
   * подробное описание
   * код
   * тест
 </WRAP>
+=== Source ===
+<Code linenums lang-r>
+require('ggtern');
+classifydiag.la_y_nb_cabanis <- function (wtp_data, ...)
+{
+ds <- selectNames(wtp_data,c('Name','La_ppm', 'Y_ppm', 'Nb_ppm'))
+pts <- data.frame(Name=ds$Name,
+                  X = ds$La_ppm/10,
+                  Y = ds$Y_ppm/15,
+                  Z = ds$Nb_ppm/8);
+labels <- c("La/10", "Y/15", "Nb/8");
+lines <- data.frame(X1 = c(0.00, 0.32, 0.43, 0.27, 0.12, 0.16, 0.22, 0.29),
+                    Y1 = c(1.00, 0.68, 0.57, 0.56, 0.80, 0.68, 0.50, 0.28),
+                    Z1 = c(0.00, 0.00, 0.00, 0.17, 0.08, 0.17, 0.28, 0.43),
+                    X2 = c(0.62, 0.24, 0.34, 0.19, 0.38, 0.00, 0.00, 0.00),
+                    Y2 = c(0.00, 0.62, 0.47, 0.60, 0.00, 0.76, 0.60, 0.42),
+                    Z2 = c(0.38, 0.15, 0.20, 0.21, 0.62, 0.24, 0.40, 0.58));
+marks  <- data.frame(Label=c("VAT", "NMORB",  "*", "EMORB", "Calc-alkali", "Cont.", "Alkaline", "Intercont. rifts"),
+                          X=c(0.24,    0.07, 0.20,    0.10,          0.70,    0.45,       0.17,              0.18),
+                          Y=c(0.70,    0.80, 0.63,    0.45,          0.10,    0.10,       0.20,              0.10),
+                          Z=c(0.07,    0.17, 0.17,    0.43,          0.20,    0.45,       0.63,              0.72));
+ret <- c();
+ret$type <- "ternogram";
+ret$points <- pts;
+ret$labels  <- labels;
+ret$lines <- lines;
+ret$marks <- marks;
+return(ret);
+}
+</Code>
+{{tag>classifydiag wholerock}}