| Version 11 (modified by , 3 years ago) (diff) |
|---|
Available options for the different filter algorithms
Contents of this page
There are different operations for each of the filter algorithms that need to be specified in the call to pdaf_init. To display the available options in a compiled assimilation program, one can use run with the specification subtype=-1. In this case pdaf_init wil display the available options for the selected filter algorithm and exit with an error status of -1.
Below we list the options as they are displayed using subtype=-1.
SEEK (filtertype=0)
PDAF Available options for SEEK: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: Evolve unit modes with finite difference approx. of TLM PDAF 1: like 0 with modes scaled by eigenvalues, unit U PDAF 2: Fixed basis vectors; variable U matrix PDAF 3: Fixed covariance matrix (V and U kept constant) PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): Interval for re-diagonalization of P (>0); optional: default 1 PDAF param_int(4): 1 for incremental updating, 0 else; optional: default 0 PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF param_real(2): epsilon for finite-difference approx. of TLM, required PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the SEEK filter ++++++++++
SEIK (filtertype=1)
PDAF Available options for SEIK: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: full ensemble integration; left-sided application of T PDAF 1: full ensemble integration; right-sided application of T PDAF 2: Fixed error space basis PDAF 3: Fixed state covariance matrix PDAF 4: Implementation with explicit ensemble transformation PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): not used PDAF param_int(4): 1 for incremental updating, 0 else; optional, default: 0 PDAF param_int(5): Type of forgetting factor; optional, default: 0 PDAF 0: fixed forgetting factor PDAF 1: adaptive forgetting factor (experimental) PDAF param_int(6): Type of ensemble transformation matrix; optional, default: 0 PDAF 0: deterministic omega PDAF 1: random orthonormal omega orthogonal to (1,...,1)^T PDAF 2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T PDAF (experimental; for random transformations, 1 is recommended) PDAF param_int(7): Type of transformation matrix square root; optional, default: 0 PDAF (Only relevant for subtype/=3) PDAF 0: symmetric square root PDAF 1: Cholesky decomposition PDAF param_int(8): Application of observation operator H, default: 0 PDAF 0: Apply H to ensemble mean to compute residual PDAF 1: Apply H to all ensemble states and then compute residual from mean of these PDAF param_int(8)=1 is the recomended choice for nonlinear H PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF --- Internal parameter (defined inside PDAF) --- PDAF Nm1vsN: Normalization of covariance matrix; default: 1 PDAF 0: normalization with 1/(Ensemble size) PDAF (original SEIK, mainly for compatibility with older studies) PDAF 1: normalization with 1/(Ensemble size - 1) PDAF (sample covariance matrix consistent with other EnKFs) PDAF +++++++++ End of option overview for the SEIK filter ++++++++++
EnKF (filtertype=2)
PDAF Available options for EnKF: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: Full ensemble integration; analysis for 2*dim_obs>dim_ens PDAF 1: Full ensemble integration; analysis for 2*dim_obs<=dim_ens PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): maximum rank for inversion of HPH^T, optional, default=0 PDAF (for =0, HPH is inverted by solving the representer equation) PDAF (if set to >=ensemble size, it is reset to ensemble size - 1) PDAF param_int(4): not used PDAF param_int(5): Size of lag for smoothing PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF (>=1; not larger than total number of processors) PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the EnKF ++++++++++
LSEIK (filtertype=3)
PDAF Available options for LSEIK: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: full ensemble integration; left-sided application of T PDAF 2: Fixed error space basis PDAF 3: Fixed state covariance matrix PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): not used PDAF param_int(4): 1 for incremental updating, 0 else; optional, default: 0 PDAF param_int(5): Type of forgetting factor; optional, default: 0 PDAF 0: fixed forgetting factor PDAF 1: adaptive forgetting factor for full domain (experimental) PDAF 2: locally adaptive forgetting factor (experimental) PDAF param_int(6): Type of ensemble transformation matrix; optional, default: 0 PDAF 0: deterministic omega PDAF 1: random orthonormal omega orthogonal to (1,...,1)^T PDAF 2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T PDAF (experimental; for random transformations, 1 is recommended) PDAF param_int(7): Type of transformation matrix square root; optional, default: 0 PDAF (Only relevant for subtype/=3) PDAF 0: symmetric square root PDAF 1: Cholesky decomposition PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF --- Internal parameter (defined inside PDAF) --- PDAF Nm1vsN: Normalization of covariance matrix; default: 1 PDAF 0: normalization with 1/(Ensemble size) PDAF (original SEIK, mainly for compatibility with older studies) PDAF 1: normalization with 1/(Ensemble size - 1) PDAF (sample covariance matrix consistent with other EnKFs) PDAF +++++++++ End of option overview for the LSEIK filter ++++++++++
ETKF (filtertype=4)
PDAF Available options for ETKF: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: full ensemble integration; apply T-matrix analogously to SEIK PDAF 1: full ensemble integration; formulation without T matrix PDAF 2: Fixed error space basis; analysis with T-matrix PDAF 3: Fixed state covariance matrix; analysis with T-matrix PDAF 5: Offline mode; analysis with T-matrix PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): Size of lag for smoothing PDAF param_int(4): not used PDAF param_int(5): Type of forgetting factor; optional, default: 0 PDAF 0: fixed forgetting factor PDAF 1: adaptive forgetting factor (experimental) PDAF param_int(6): Type of ensemble transformation matrix; optional, default: 0 PDAF 0: deterministic transformation PDAF 2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T PDAF param_int(7): not used PDAF param_int(8): Application of observation operator H, default: 0 PDAF 0: Apply H to ensemble mean to compute residual PDAF 1: Apply H to all ensemble states and then compute residual from mean of these PDAF param_int(8)=1 is the recomended choice for nonlinear H PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the ETKF ++++++++++
LETKF (filtertype=5)
PDAF Available options for LETKF: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: full ensemble integration; apply T-matrix analogously to SEIK PDAF 2: Fixed error space basis; analysis with T-matrix PDAF 3: Fixed state covariance matrix; analysis with T-matrix PDAF 5: Offline mode; analysis with T-matrix PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): Size of lag for smoothing PDAF param_int(4): not used PDAF param_int(5): Type of forgetting factor; optional, default: 0 PDAF 0: fixed forgetting factor PDAF 1: adaptive forgetting factor for full domain (experimental) PDAF 2: locally adaptive forgetting factor (experimental) PDAF param_int(6): Type of ensemble transformation matrix; optional, default: 0 PDAF 0: deterministic transformation PDAF 2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the LETKF ++++++++++
ESTKF (filtertype=6)
PDAF Available options for ESTKF: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: Standard implementation with ensemble integration PDAF 2: Fixed error space basis PDAF 3: Fixed state covariance matrix PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): Size of lag for smoothing PDAF param_int(4): not used PDAF param_int(5): Type of forgetting factor; optional, default: 0 PDAF 0: fixed forgetting factor PDAF 1: adaptive forgetting factor (experimental) PDAF param_int(6): Type of ensemble transformation matrix; optional, default: 0 PDAF 0: deterministic omega PDAF 1: random orthonormal omega orthogonal to (1,...,1)^T PDAF 2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T PDAF (experimental; for random transformations, 0 or 1 are recommended) PDAF param_int(7): Type of transformation matrix square root; optional, default: 0 PDAF 0: symmetric square root PDAF 1: Cholesky decomposition PDAF param_int(8): Application of observation operator H, default: 0 PDAF 0: Apply H to ensemble mean to compute residual PDAF 1: Apply H to all ensemble states and then compute residual from mean of these PDAF param_int(8)=1 is the recomended choice for nonlinear H PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the ESTKF ++++++++++
LESTKF (filtertype=7)
PDAF Available options for LESTKF: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: Standard implementation with ensemble integration PDAF 2: Fixed error space basis PDAF 3: Fixed state covariance matrix PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): Size of lag for smoothing PDAF param_int(4): not used PDAF param_int(5): Type of forgetting factor; optional, default: 0 PDAF 0: fixed forgetting factor PDAF 1: adaptive forgetting factor for full domain (experimental) PDAF 2: locally adaptive forgetting factor (experimental) PDAF param_int(6): Type of ensemble transformation matrix; optional, default: 0 PDAF 0: deterministic omega PDAF 1: random orthonormal omega orthogonal to (1,...,1)^T PDAF 2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T PDAF (experimental; for random transformations, 1 is recommended) PDAF param_int(7): Type of transformation matrix square root; optional, default: 0 PDAF 0: symmetric square root PDAF 1: Cholesky decomposition PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the LESTKF ++++++++++
LEnKF (filtertype=8)
PDAF Available options for LEnKF: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: Full ensemble integration; analysis with covariance localization PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): maximum rank for inversion of HPH^T, optional, default=0 PDAF (for =0, HPH is inverted by solving the representer equation) PDAF (if set to >=ensemble size, it is reset to ensemble size - 1) PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF (>=1; not larger than total number of processors) PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the LEnKF ++++++++++
NETF (filtertype=9)
PDAF Available options for NETF: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: Standard implementation with ensemble integration PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): Size of lag for smoothing PDAF param_int(4): not used PDAF param_int(5): Type of forgetting factor; optional, default: 0 PDAF 0: forgetting factor on forecast ensemble PDAF 2: forgetting factor on analysis ensemble PDAF param_int(6): Type of ensemble transformation matrix; optional, default: 0 PDAF 0: random orthonormal matrix orthogonal to (1,...,1)^T PDAF 1: deterministic transformation PDAF param_int(7): Type of weights inflation; optional, default: 0 PDAF 0: no weights inflation PDAF 1: inflate so that N_eff/N > param_real(2) PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF param_real(2): Limit for weigts inflation N_eff/N > param_real(2), optional, default=0.0 PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the NETF ++++++++++
LNETF (filtertype=10)
PDAF Available options for LNETF: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: Standard implementation with ensemble integration PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): Size of lag for smoothing PDAF param_int(4): not used PDAF param_int(5): Type of forgetting factor; optional, default: 0 PDAF 0: forgetting factor on forecast ensemble PDAF 1: forgetting factor on forecast ensemble only observed domains PDAF 2: forgetting factor on analysis ensemble PDAF 3: forgetting factor on analysis ensemble only observed domains PDAF param_int(6): Type of ensemble transformation matrix; optional, default: 0 PDAF 0: random orthonormal matrix orthogonal to (1,...,1)^T PDAF 1: deterministic transformation PDAF param_int(7): Type of weights inflation; optional, default: 0 PDAF 0: no weights inflation PDAF 1: inflate so that N_eff/N > param_real(2) PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required PDAF param_real(2): Limit for weigts inflation N_eff/N > param_real(2), optional, default=0.0 PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the LNETF ++++++++++
GENOBS (filtertype=11)
PDAF Available options for GENOBS: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: Standard implementation with ensemble integration PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Forgetting factor (usually >0 and <=1), required, but not used PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF =1 for GENOBS; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for GENOBS ++++++++++
PF (filtertype=12)
PDAF Available options for PF: PDAF --- Sub-types (Parameter subtype) --- PDAF 0: Standard implementation with ensemble integration PDAF 5: Offline mode PDAF --- Integer parameters (Array param_int) --- PDAF param_int(1): Dimension of state vector (>0), required PDAF param_int(2): Ensemble size (>0), required PDAF param_int(3): Resampling type, optional PDAF 1: probabilistic resamping (default) PDAF 2: stochastic universal resampling PDAF 3: residual resampling PDAF param_int(4): Type of ensemble perturbations, optional PDAF 0: no perturbations (default) PDAF 1: constant standard deviation PDAF 2: relative to ensemble standard deviation PDAF --- Floating point parameters (Array param_real) --- PDAF param_real(1): Ensemble pert. level (>0), required, only used if param_int(4)>0 PDAF --- Further parameters --- PDAF n_modeltasks: Number of parallel model integration tasks PDAF >=1 for subtypes 0 and 1; not larger than total number of processors PDAF =1 required for subtypes 2 and 3 PDAF screen: Control verbosity of PDAF PDAF 0: no outputs PDAF 1: basic output (default) PDAF 2: 1 plus timing output PDAF 3: 2 plus debug output PDAF +++++++++ End of option overview for the PF ++++++++++
