l96_integrate: Integrate state forward with lorenz96

In robertsy/lorenzr: Lorenz Models

Description

Time integration with lorenz96: implemented with a euler scheme in fortran90.

Usage

l96_integrate(state, nsteps = 1, deltat = 1/1000, ndim = 40,
  Forcing = 8, ...)

Arguments

state	initial state (can be an ensemble matrix of ndim x K)
nsteps	how many steps ahead
deltat	time interval for model integration
ndim	dimension of system
Forcing	external forcing on dynamics
...	additional arguments passed to lorenz_observe (e.g. nobs in case obs_type=partial)

Value

state at the end of integration (either a vector or a matrix depending on initial state)

Examples

current_state <- rnorm(40)
par(mfrow=c(4,2), mai=c(0.2,0.4,.5,.2))
for (i in 1:8){
 plot.ts(current_state, main=paste('Time=',i), xlab='', ylab='')
 current_state <- l96_integrate( current_state, 500, deltat=1/1000)
}

robertsy/lorenzr documentation built on May 27, 2019, 10:33 a.m.