OrbitalError: Get the orbital error components

In EarthSystemDiagnostics/psem: Proxy Spectral Error Model

Get the orbital error components

Description

Get the orbital error components

Usage

OrbitalError(
  nu = NULL,
  tau_s,
  tau_b,
  tau_p,
  delta_t,
  T,
  sig.sq_c,
  phi_c,
  delta_phi_c,
  nu_c,
  sig.sq_a,
  phi_a,
  nu_a,
  N
)

Arguments

nu	frequency
tau_s	sediment slice thickness in years (layer.width / sedimentation rate)
tau_b	timescale of bioturbation (bioturbation depth / sedimentation rate) (L/sr)
tau_p	length of proxy carrier "growth season" (< 12 months)
delta_t	sampling frequency of the sediment core / climate timeseries
T	length of a finite time series, e.g. 1e04
sig.sq_c	variance of the seasonal cycle
phi_c	phase of carrier growth season relative to the maximum of the seasonal cycle
delta_phi_c	Uncertainty in the phase of the signal carrier production. A value between 0 and 2Pi.
nu_c	1/1 (yr) = frequency of the seasonal cycle
sig.sq_a	variance of precessionary amplitude modulation
phi_a	phase of precessionary amplitude modulation relative to centre of finite timeseries of length T
nu_a	1/tau_a = frequency of the orbital variation, e.g. for precession 1/23e03 yrs
N	number of signal carriers (e.g. individual foraminifera)

Value

A list with the frequency domain and time domain error components

Examples

spec.pars <- GetSpecPars("Mg_Ca", phi_c = pi, delta_phi_c = pi, sig.sq_a = 0.1, tau_p = 1/12)
do.call(OrbitalError, spec.pars[names(spec.pars) %in% names(formals(OrbitalError))])

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