RelativeAmplitudeModulation: Return the fraction of the orbital variations in the seasonal...

In EarthSystemDiagnostics/psem: Proxy Spectral Error Model

Return the fraction of the orbital variations in the seasonal amplitude relative to the mean seasonal amplitude

Description

Return the fraction of the orbital variations in the seasonal amplitude relative to the mean seasonal amplitude

Usage

RelativeAmplitudeModulation(
  latitude,
  maxTimeKYear = 100,
  minTimeKYear = 1,
  bPlot = FALSE
)

Arguments

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

Value

fraction of orbital vs. seasonal amplitude

Author(s)

Thomas Laepple

Examples

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)

EarthSystemDiagnostics/psem documentation built on March 3, 2024, 10:24 p.m.