LorenzUV

The 2-scale/layer/speed coupled version of Lorenz-96.

See wilks2005 - U: large amp, low frequency vars: convective events - V: small amp, high frequency vars: large-scale synoptic events

Typically, the DA system will only use the truncated system (containing only the U variables), where the V's are parameterized as model noise, while the truth is simulated by the full system.

Stochastic parmateterization: Wilks: benefit of including stochastic noise negligible unless its temporal auto-corr is taken into account (as AR(1)) (but spatial auto-corr can be neglected). But AR(1) noise is technically difficult because DAPPER is built around the Markov assumption.

Modules:

Name	Description
L96	A 1D emulator of chaotic atmospheric behaviour.
LP	On-line (live) plots of the DA process for various models and methods.
demo	Demonstrate the Lorenz two-speed/scale/layer model.
illust_LorenzUV	Prettier, static illustration of Lorenz two-speed/scale/layer model.
illust_parameterizations	Illusrate parameterizations.
lorenz96	As in lorenz1996predictability.
wilks05	Uses nU, J, F as in mods.LorenzUV ie. from wilks2005.

model_instance

Use OOP to facilitate having multiple parameter settings simultaneously.

Default parameters from wilks2005.

dxdt(x)

Compute full (coupled) dxdt.

dxdt_parameterized(x, t)

Compute truncated dxdt with parameterization of fast variables (V).

dxdt_trunc(x)

Compute truncated dxdt: slow variables (U) only.

reversible(fun)

Reverse input/output (instead of manipulating indices).