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

Thomas Laepple

Examples

RelativeAmplitudeModulation(34)
#> $mean.amp
#> [1] 303.4541
#> 
#> $max.amp
#> [1] 348.3554
#> 
#> $min.amp
#> [1] 257.7186
#> 
#> $relOrbitalAmplitude
#> [1] 0.2986839
#> 
#> $relOrbitalMean
#> [1] 0.2570214
#> 
#> $sig_a
#> [1] 0.1046287
#> 
#> $sig.sq_a
#> [1] 0.01094716
#> 
if (FALSE) {
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)))
}