SimulateYearFromInsolation: Simulate a year from a parametric relationship to insolation

In EarthSystemDiagnostics/orbitalforcing: Orbital Forcing from Insolation

Simulate a year from a parametric relationship to insolation

Description

Simulate a year from a parametric relationship to insolation

Usage

SimulateYearFromInsolation(
  kyear = 0,
  transfer,
  latitude,
  bPolynomial = FALSE,
  anomalyRef = NULL
)

Arguments

kyear	point in time (kyr BP)
transfer	transfer function
latitude	latitude of the insolation calculation
bPolynomial	(TRUE) use polynomial, FALSE = linear
anomalyRef	NULL for absolute mode or vector of 365 anomaly correction values

Value

365 predicted values for each day of the year

Author(s)

Thomas Laepple

Examples

library(orbitalforcing)
clim <- orbitalforcing:::SelSpace3D(sat.ncep.clim, lat1 = 40, lon1 = 200)
transfer <- EstimateTransferfunctionInsolation(clim, 40, bPlot = FALSE, b3plot = FALSE)
plot(1:365, SimulateYearFromInsolation(0, transfer, 40), type = "l",
     main = "predicted at 0 and 10kyr BP", xlab = "day of the year")
lines(1:365, SimulateYearFromInsolation(10, transfer, 40), col = "red")

EarthSystemDiagnostics/orbitalforcing documentation built on March 24, 2022, 11:25 a.m.