Skip to content

QG⚓︎

Quasi-geostraphic 2D flow. Described in detail by sakov2008b.

Adapted from Pavel Sakov's enkf-matlab package.

More info:

  • governing_eqn.png
  • demo.py
  • ψ (psi) is the stream function (i.e. surface elevation)
  • Doubling time "between 25 and 50"
  • Note Sakov's trick of increasing RKH2 from 2.0e-12 to 2.0e-11 to stabilize the ensemble integration, which may be necessary for EnKF's with small N. See example in counillon2009.

Modules:

Name Description
LP

On-line (live) plots of the DA process for various models and methods.

counillon2009

Reproduce experiments from counillon2009a.

demo

Demonstrate the QG (quasi-geostrophic) model.

f90

This dir contains the Fortran-90 code of the model.

illust_obs

Stream function and observation time series for QG (quasi-geostrophic) model.

modelling

Contains models included with DAPPER.

sakov2008

Reproduce results from sakov2008b.

model_config ⚓︎

Define model.

Helps ensure consistency between prms file (that Fortran module reads) and Python calls to step(), for example for dt.

__init__(name, prms, mp=True) ⚓︎

Use prms={} to get the default configuration.

step(E, t, dt) ⚓︎

Vector and 2D-array (ens) input, with multiproc for ens case.

step_1(x0, t, dt) ⚓︎

Step a single state vector.