View source: R/additional-functions.R
RelativeAmplitudeModulation | R Documentation |
Return the fraction of the orbital variations in the seasonal amplitude relative to the mean seasonal amplitude
RelativeAmplitudeModulation(
latitude,
maxTimeKYear = 100,
minTimeKYear = 1,
bPlot = FALSE
)
latitude |
Latitude in degN |
maxTimeKYear |
Maximum time in kyr BP which should be analysed |
minTimeKYear |
Minimum time in kyr BP which should be analysed |
bPlot |
logical, plot the max, min and mean daily insolation |
fraction of orbital vs. seasonal amplitude
Thomas Laepple
RelativeAmplitudeModulation(34)
## Not run:
library(dplyr)
library(tidyr)
df <- crossing(lat = seq(-90, 90, by = 10), maxT = c(23, 100, 1000)) %>%
group_by(lat, maxT) %>%
mutate(sig_a = RelativeAmplitudeModulation(lat, maxTimeKYear = maxT, minTimeKYear = 0)$sig_a) %>%
ungroup() %>%
mutate(A_a = sig_a * sqrt(2))
df %>%
ggplot(aes(x = lat, y = A_a)) +
geom_line(aes(colour = factor(maxT)))
df %>%
ggplot(aes(x = lat, y = sig_a^2)) +
geom_line(aes(colour = factor(maxT)))
df %>%
ggplot(aes(x = lat, y = sig_a)) +
geom_line(aes(colour = factor(maxT)))
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.