ModelSpectrum: Create a model SST spectrum based on the observed spectra and...

In EarthSystemDiagnostics/psem: Proxy Spectral Error Model

Create a model SST spectrum based on the observed spectra and a powerlaw spectrum

Description

Create a model SST spectrum based on the observed spectra and a powerlaw spectrum

Usage

ModelSpectrum(
  freq = NULL,
  latitude,
  spectable = psem:::spectable,
  beta = 1,
  variable = c("temperature", "d18O", "T_deg_Mg_Ca", "T_deg_Uk37"),
  freq.match.lower = 0.03,
  freq.match.upper = 0.1,
  bPlot = FALSE
)

Arguments

freq	desired frequency vector
latitude	latitude in degN
spectable	lookup table with spectra by latitude list(spec[iLat,iFreq],lat[iLat],freq[iFreq])
beta	desired scaling
variable	climate variable for which to return the power spectrum, currently temperature or d18O
freq.match.lower	lower frequency boundary at which the variance is matched and the spectrum is switched from instrumental to powerlaw (default 0.03 as lowest reliable estimate from the instrumental record given the detrending and the multitaper bias)
freq.match.upper	higher frequency boundary at which the variance is matched (default 0.1 = decadal)
bPlot	TRUE = provide diagnostic plot

Value

model spectrum list(spec,freq)

Author(s)

Thomas Laepple

Examples

ModelSpectrum(latitude = 20,beta=1, bPlot=TRUE)
ModelSpectrum(freq = GetNu(T = 10000, delta_t = 100), latitude = 20, beta=1, bPlot=FALSE)
ModelSpectrum(freq = c(0.1, GetNu(T = 10000, delta_t = 100), 0.1), latitude = 20, beta=1, bPlot=FALSE)
ModelSpectrum(freq = c(0.6, 1), latitude = 20, beta=1, bPlot=FALSE)

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