propagate_d: Propagate determinstic model forward over a series ot...

In mcrucifix/iceages: Ice age models and analysis

Propagate determinstic model forward over a series ot timesteps

Description

Essentially a wrapper to propagate_1step_D. Propagate determinstic model forward over a series ot timesteps. The numerical schemeis supplied by the model (see original codes)

Usage

propagate_d(model, times, init, par, Astro, ...)

Arguments

model

One of the deterministic models out of the packaged model list (data(models)) or user-supplied.

times	Vector of times. times[1] at which the model trajectory is being evalueted. Is used as the initial condition
init	State vector of initial conditions

par	Parameter vector
Astro	Astronomical forcing list, as supplied by routine read_astro
...	Other parameters supplied to propagate_1step_D

Value

Matrix, with columns os follows:

• Column 1 .. n:States at specified times (n is the system dimension)

• Column n+1:precession according to the astronomical solution

• Column n+2:obliquity according to the astronomical solution

Author(s)

M. Crucifix (2012)

See Also

propagate_s (stochastic version)

Examples

## check convergence of time series
data(models)

times = seq(-80,0,0.2)
deltat = c(0.01, 0.02, 0.1)

# uses 12 components of precession and obliquity
Astro <- read_astro(12,12)

sol = list()
for (i in seq(along=deltat))
 {
  sol[[i]] <- propagate_d (models$vdp_d, models$vdp_d$spar, init=c(0.1, 0.1), 
                          times=times, deltat=deltat[i], Astro=Astro)
 }

plot(times, sol[[1]][,1], type='l')
for (i in seq(along=sol)[-1]) 
 { 
  lines(times, sol[[i]][,1], type='l', col=i)
 }

mcrucifix/iceages documentation built on Jan. 11, 2023, 9:17 p.m.