View source: R/main_interface.R
propagate_d | R Documentation |
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)
propagate_d(model, times, init, par, Astro, ...)
model |
One of the deterministic models out of the packaged model list ( |
times |
Vector of times. |
init |
State vector of initial conditions |
par |
Parameter vector |
Astro |
Astronomical forcing list, as supplied by routine |
... |
Other parameters supplied to |
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
M. Crucifix (2012)
propagate_s
(stochastic version)
## 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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.