Различия
Здесь показаны различия между двумя версиями данной страницы.
| Следующая версия | Предыдущая версия | ||
|
spdgrams:ree_nmorb [2016/03/11 10:30] admin создано |
spdgrams:ree_nmorb [2018/11/12 15:07] (текущий) |
||
|---|---|---|---|
| Строка 1: | Строка 1: | ||
| - | ====== N-MORB ====== | + | ====== N-MORB, Sun-McDonough, 1989 ====== |
| === Details === | === Details === | ||
| - | <data plot> | + | | Reference | <ref:sun1989> | |
| - | Reference [wiki]: | + | | Data | bulk | |
| - | Data: bulk | + | | Input | ppm | |
| - | Input: ppm | + | |
| - | </data> | + | <use:spdgrams.ree_nmorb> |
| + | <test:spdgrams.ree_nmorb> | ||
| <WRAP center round todo 60%> | <WRAP center round todo 60%> | ||
| - | * заголовок | + | * <del>заголовок</del> |
| * краткое описание | * краткое описание | ||
| - | * статья | + | * <del>статья</del> |
| - | * ссылка | + | * <del>ссылка</del> |
| * <del>карточка</del> | * <del>карточка</del> | ||
| * подробное описание | * подробное описание | ||
| - | * код | + | * <del>код</del> |
| * тест | * тест | ||
| </WRAP> | </WRAP> | ||
| Строка 25: | Строка 26: | ||
| <Code linenums lang-r> | <Code linenums lang-r> | ||
| - | require('ggplot'); | + | require('ggplot2'); |
| - | spdgrams.ree_primitive_mantle <- function(ppm_data) | + | cmp <- read.csv("constants_oceanic_basalts_compositions.csv", sep=";", stringsAsFactors=F); |
| + | |||
| + | spdgrams.ree_nmorb <- function(ppm_data, el_order = NULL) | ||
| { | { | ||
| - | primitive_mantle = data.frame(Cs = c(0.007), Rb = c(0.56), Ba = c(6.3), Th = c(0.12), | + | ns <- el_order; |
| - | U = c(0.047), Nb = c(2.33), Ta = c(0.132), La = c(2.5), | + | |
| - | Ce = c(7.5), Pb = c(0.3), Pr = c(1.32), Sr = c(90), | + | if (is.null(el_order)) |
| - | Nd = c(7.3), Zr = c(74), Hf = c(2.05), Sm = c(2.63), | + | { |
| - | Eu = c(1.02), Gd = c(3.68), Tb = c(0.67), Dy = c(4.55), | + | ns <- c('Cs', 'Rb', 'Ba', 'Th', 'U', 'Nb', 'Ta', 'La', 'Ce', 'Pb', 'Pr', |
| - | Y = c(28), Ho = c(1.01), Er = c(2.97), Tm = c(0.456), | + | 'Sr', 'Nd', 'Zr', 'Hf', 'Sm', 'Eu', 'Gd', 'Tb', 'Dy', 'Y', 'Ho', |
| - | Yb = c(3.05), Lu = c(0.455)); | + | 'Er', 'Tm', 'Yb', 'Lu'); |
| + | } | ||
| - | ns <- names(primitive_mantle); | + | pmsc <- c(); |
| + | |||
| + | for (n in ns) | ||
| + | pmsc <- c(pmsc, subset(cmp, element == n)$N_MORB); | ||
| ds <- selectNames(ppm_data,c(paste0(ns,'_ppm'))); | ds <- selectNames(ppm_data,c(paste0(ns,'_ppm'))); | ||
| Строка 45: | Строка 52: | ||
| ret <- c(); | ret <- c(); | ||
| - | values <- cbind(ppm_data$Name, rowApply(`/`,ds,primitive_mantle[1,])); | + | values <- cbind(ppm_data$Name, rowApply(`/`,ds,pmsc)); |
| names(values)[1] <- 'Name'; | names(values)[1] <- 'Name'; | ||
| ret$values <- values; | ret$values <- values; | ||
| ret$grades <- ns; | ret$grades <- ns; | ||
| - | ret$labels <- c('REE', 'Sample/Primitive mantle'); | + | ret$labels <- c('REE', 'Sample/N-MORB'); |
| return(ret); | return(ret); | ||
| } | } | ||
| + | |||
| </Code> | </Code> | ||
| + | |||
| + | {{tag>spdgrams oceanic}} | ||