Changes between Version 1 and Version 2 of AvailableOptionsforInitPDAFinPDAF3


Ignore:
Timestamp:
Mar 24, 2025, 11:11:36 AM (8 months ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • AvailableOptionsforInitPDAFinPDAF3

    v1 v2  
    9191
    9292
    93 == Option outputs using subtype=-1 ==
    94 
    95 === SEIK (filtertype=PDAF_DA_SEIK=1) ===
    96 
    97 {{{
    98 PDAF     Available options for SEIK:
    99 PDAF     --- Sub-types (Parameter subtype) ---
    100 PDAF       0: full ensemble integration; left-sided application of T
    101 PDAF       1: full ensemble integration; right-sided application of T
    102 PDAF       2: full ensemble integration; explicit ensemble transformation
     93== Option outputs using subtype=-1 - Local Ensemble Kalman filters ==
     94
     95
     96
     97
     98=== LESTKF (filtertype=PDAF_DA_LESTKF=7) ===
     99
     100{{{
     101PDAF     Available options for LESTKF:
     102PDAF     --- Sub-types (Parameter subtype) ---
     103PDAF       0: Standard implementation with ensemble integration
    103104PDAF       10: Fixed error space basis
    104105PDAF       11: Fixed state covariance matrix
    … …  
    106107PDAF       param_int(1): Dimension of state vector (>0), required
    107108PDAF       param_int(2): Ensemble size (>0), required
    108 PDAF       param_int(3): not used
     109PDAF       param_int(3): dim_lag
     110PDAF           Size of smoothing lag (>=0), optional
     111PDAF            0: no smoothing (default)
     112PDAF            >0: apply smoother up to specified lag
    109113PDAF       param_int(4): not used
    110114PDAF       param_int(5) type_forget
    111115PDAF           Type of forgetting factor; optional
    112116PDAF            0: fixed forgetting factor (default)
    113 PDAF            1: adaptive forgetting factor (experimental)
     117PDAF            1: adaptive forgetting factor for full domain
     118PDAF            2: locally adaptive forgetting factor
    114119PDAF       param_int(6) type_trans
    115120PDAF           Type of ensemble transformation matrix; optional
    … …  
    120125PDAF       param_int(7) type_sqrt
    121126PDAF           Type of transformation matrix square root; optional
    122 PDAF            (Only relevant for subtype/=11)
    123127PDAF            0: symmetric square root (default)
    124128PDAF            1: Cholesky decomposition
    … …  
    137141PDAF     --- Further parameters ---
    138142PDAF       n_modeltasks: Number of parallel model integration tasks
    139 PDAF           >=1 for subtypes 0, 1 and 2; not larger than total number of processors
     143PDAF           >=1 for subtype 0; not larger than total number of processors
    140144PDAF           =1 required for subtypes 10 and 11
    141145PDAF       screen: Control verbosity of PDAF
    … …  
    144148PDAF           2: 1 plus timing output
    145149PDAF           3: 2 plus debug output
    146 PDAF     --- Internal parameter (defined inside PDAF) ---
    147 PDAF       Nm1vsN: Normalization of covariance matrix; default: 1
    148 PDAF           0: normalization with 1/(Ensemble size)
    149 PDAF              (original SEIK, mainly for compatibility with older studies)
    150 PDAF           1: normalization with 1/(Ensemble size - 1)
    151 PDAF           (sample covariance matrix consistent with other EnKFs)
    152 PDAF     +++++++++ End of option overview for the SEIK filter ++++++++++
    153 }}}
    154 
    155 
    156 === EnKF (filtertype=PDAF_DA_ENKF=2) ===
    157 
    158 {{{
    159 PDAF     Available options for EnKF:
    160 PDAF     --- Sub-types (Parameter subtype) ---
    161 PDAF       0: Full ensemble integration; analysis for 2*dim_obs>dim_ens
    162 PDAF       1: Full ensemble integration; analysis for 2*dim_obs<=dim_ens
     150PDAF     +++++++++ End of option overview for the LESTKF ++++++++++
     151}}}
     152
     153
     154=== LETKF (filtertype=PDAF_DA_LETKF=5) ===
     155
     156{{{
     157PDAF     Available options for LETKF:
     158PDAF     --- Sub-types (Parameter subtype) ---
     159PDAF       0: full ensemble integration;  apply T-matrix analogously to SEIK
     160PDAF       1: full ensemble integration; formulation cf. Hunt et al. (2007) without T matrix
     161PDAF       10: Fixed error space basis; analysis with T-matrix
     162PDAF       11: Fixed state covariance matrix; analysis with T-matrix
    163163PDAF     --- Integer parameters (Array param_int) ---
    164164PDAF       param_int(1): Dimension of state vector (>0), required
    … …  
    166166PDAF       param_int(3): dim_lag
    167167PDAF           Size of smoothing lag (>=0), optional
    168 PDAF           0: no smoothing (default)
    169 PDAF           >0: apply smoother up to specified lag
     168PDAF            0: no smoothing (default)
     169PDAF            >0: apply smoother up to specified lag
     170PDAF       param_int(4): not used
     171PDAF       param_int(5) type_forget
     172PDAF           Type of forgetting factor; optional
     173PDAF            0: fixed forgetting factor (default)
     174PDAF            1: adaptive forgetting factor for full domain
     175PDAF            2: locally adaptive forgetting factor
     176PDAF       param_int(6) type_trans
     177PDAF           Type of ensemble transformation matrix; optional
     178PDAF            0: deterministic transformation (default)
     179PDAF            2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T
     180PDAF       param_int(8): observe_ens
     181PDAF           Application of observation operator H, optional
     182PDAF            0: Apply H to ensemble mean to compute innovation
     183PDAF            1: Apply H to ensemble states; then compute innovation from their mean (default)
     184PDAF               param_int(8)=1 is the recomended choice for nonlinear H
     185PDAF       param_int(9): type_obs_init
     186PDAF           Initialize observations before or after call to prepoststep_pdaf
     187PDAF            0: Initialize observations before call to prepoststep_pdaf
     188PDAF            1: Initialize observations after call to prepoststep_pdaf (default)
     189PDAF     --- Floating point parameters (Array param_real) ---
     190PDAF       param_real(1): forget
     191PDAF           Forgetting factor (usually >0 and <=1), required
     192PDAF     --- Further parameters ---
     193PDAF       n_modeltasks: Number of parallel model integration tasks
     194PDAF           >=1 for subtypes 0 and 1; not larger than total number of processors
     195PDAF           =1 required for subtypes 10 and 11
     196PDAF       screen: Control verbosity of PDAF
     197PDAF           0: no outputs
     198PDAF           1: basic output (default)
     199PDAF           2: 1 plus timing output
     200PDAF           3: 2 plus debug output
     201PDAF     +++++++++ End of option overview for the LETKF ++++++++++
     202}}}
     203
     204
     205=== LEnKF (filtertype=PDAF_DA_LENKF=8) ===
     206
     207{{{
     208PDAF     Available options for LEnKF:
     209PDAF     --- Sub-types (Parameter subtype) ---
     210PDAF       0: Standard EnKF analysis with covariance localization
     211PDAF     --- Integer parameters (Array param_int) ---
     212PDAF       param_int(1): Dimension of state vector (>0), required
     213PDAF       param_int(2): Ensemble size (>0), required
     214PDAF       param_int(3): not used
    170215PDAF       param_int(4): rank_ana_enkf
    171216PDAF           maximum rank for inversion of HPH^T, optional, default=0
    … …  
    195240PDAF           2: 1 plus timing output
    196241PDAF           3: 2 plus debug output
    197 PDAF     +++++++++ End of option overview for the EnKF ++++++++++
     242PDAF     +++++++++ End of option overview for the LEnKF ++++++++++
     243}}}
     244
     245
     246
     247=== ENSRF/EAKF (filtertype=PDAF_DA_ENSRF=13) ===
     248
     249{{{
     250PDAF     Available options for ENSRF:
     251PDAF     --- Sub-types (Parameter subtype) ---
     252PDAF       0: ENSRF with serial observation processing (cf. Houtekamer/Mitchell, 2002)
     253PDAF       1: EAKF/2-step local least squares filter (cf. Anderson, 2003)
     254PDAF     --- Integer parameters (Array param_int) ---
     255PDAF       param_int(1): Dimension of state vector (>0), required
     256PDAF       param_int(2): Ensemble size (>0), required
     257PDAF       param_int(3): not used
     258PDAF       param_int(4): not used
     259PDAF       param_int(5): not used
     260PDAF       param_int(6): not used
     261PDAF       param_int(7): not used
     262PDAF       param_int(8): observe_ens
     263PDAF           Application of observation operator H, optional
     264PDAF            0: Apply H to ensemble mean to compute innovation
     265PDAF            1: Apply H to ensemble states; then compute innovation from their mean (default)
     266PDAF               param_int(8)=1 is the recomended choice for nonlinear H
     267PDAF       param_int(9): type_obs_init
     268PDAF           Initialize observations before or after call to prepoststep_pdaf
     269PDAF            0: Initialize observations before call to prepoststep_pdaf
     270PDAF            1: Initialize observations after call to prepoststep_pdaf (default)
     271PDAF     --- Floating point parameters (Array param_real) ---
     272PDAF       param_real(1): forget
     273PDAF           Forgetting factor (usually >0 and <=1), required
     274PDAF     --- Further parameters ---
     275PDAF       n_modeltasks: Number of parallel model integration tasks
     276PDAF           (>=1; not larger than total number of processors)
     277PDAF       screen: Control verbosity of PDAF
     278PDAF           0: no outputs
     279PDAF           1: basic output (default)
     280PDAF           2: 1 plus timing output
     281PDAF           3: 2 plus debug output
     282PDAF     +++++++++ End of option overview for the ENSRF ++++++++++
    198283}}}
    199284
    … …  
    260345
    261346
    262 === ETKF (filtertype=PDAF_DA_ETKF=4) ===
    263 
    264 {{{
    265 PDAF     Available options for ETKF:
    266 PDAF     --- Sub-types (Parameter subtype) ---
    267 PDAF       0: full ensemble integration; apply T-matrix analogously to SEIK
    268 PDAF       1: full ensemble integration; formulation cf. Hunt et al. (2007) without T matrix
    269 PDAF       10: Fixed error space basis; analysis with T-matrix
    270 PDAF       11: Fixed state covariance matrix; analysis with T-matrix
    271 PDAF     --- Integer parameters (Array param_int) ---
    272 PDAF       param_int(1): Dimension of state vector (>0), required
    273 PDAF       param_int(2): Ensemble size (>0), required
    274 PDAF       param_int(3): dim_lag
    275 PDAF           Size of smoothing lag (>=0), optional
    276 PDAF            0: no smoothing (default)
    277 PDAF            >0: apply smoother up to specified lag
    278 PDAF       param_int(4): not used
    279 PDAF       param_int(5) type_forget
    280 PDAF           Type of forgetting factor; optional
    281 PDAF            0: fixed forgetting factor (default)
    282 PDAF            1: adaptive forgetting factor (experimental)
    283 PDAF       param_int(6) type_trans
    284 PDAF           Type of ensemble transformation matrix; optional
    285 PDAF            0: deterministic transformation (default)
    286 PDAF            2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T
    287 PDAF       param_int(7): not used
    288 PDAF       param_int(8): observe_ens
    289 PDAF           Application of observation operator H, optional
    290 PDAF            0: Apply H to ensemble mean to compute innovation
    291 PDAF            1: Apply H to ensemble states; then compute innovation from their mean (default)
    292 PDAF               param_int(8)=1 is the recomended choice for nonlinear H
    293 PDAF       param_int(9): type_obs_init
    294 PDAF           Initialize observations before or after call to prepoststep_pdaf
    295 PDAF           0: Initialize observations before call to prepoststep_pdaf
    296 PDAF           1: Initialize observations after call to prepoststep_pdaf (default)
    297 PDAF     --- Floating point parameters (Array param_real) ---
    298 PDAF       param_real(1): forget
    299 PDAF           Forgetting factor (usually >0 and <=1), required
    300 PDAF     --- Further parameters ---
    301 PDAF       n_modeltasks: Number of parallel model integration tasks
    302 PDAF           >=1 for subtypes 0 and 1; not larger than total number of processors
    303 PDAF           =1 required for subtypes 10 and 11
    304 PDAF       screen: Control verbosity of PDAF
    305 PDAF           0: no outputs
    306 PDAF           1: basic output (default)
    307 PDAF           2: 1 plus timing output
    308 PDAF           3: 2 plus debug output
    309 PDAF     +++++++++ End of option overview for the ETKF ++++++++++
    310 }}}
    311 
    312 
    313 
    314 
    315 === LETKF (filtertype=PDAF_DA_LETKF=5) ===
    316 
    317 {{{
    318 PDAF     Available options for LETKF:
    319 PDAF     --- Sub-types (Parameter subtype) ---
    320 PDAF       0: full ensemble integration;  apply T-matrix analogously to SEIK
    321 PDAF       1: full ensemble integration; formulation cf. Hunt et al. (2007) without T matrix
    322 PDAF       10: Fixed error space basis; analysis with T-matrix
    323 PDAF       11: Fixed state covariance matrix; analysis with T-matrix
    324 PDAF     --- Integer parameters (Array param_int) ---
    325 PDAF       param_int(1): Dimension of state vector (>0), required
    326 PDAF       param_int(2): Ensemble size (>0), required
    327 PDAF       param_int(3): dim_lag
    328 PDAF           Size of smoothing lag (>=0), optional
    329 PDAF            0: no smoothing (default)
    330 PDAF            >0: apply smoother up to specified lag
    331 PDAF       param_int(4): not used
    332 PDAF       param_int(5) type_forget
    333 PDAF           Type of forgetting factor; optional
    334 PDAF            0: fixed forgetting factor (default)
    335 PDAF            1: adaptive forgetting factor for full domain
    336 PDAF            2: locally adaptive forgetting factor
    337 PDAF       param_int(6) type_trans
    338 PDAF           Type of ensemble transformation matrix; optional
    339 PDAF            0: deterministic transformation (default)
    340 PDAF            2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T
    341 PDAF       param_int(8): observe_ens
    342 PDAF           Application of observation operator H, optional
    343 PDAF            0: Apply H to ensemble mean to compute innovation
    344 PDAF            1: Apply H to ensemble states; then compute innovation from their mean (default)
    345 PDAF               param_int(8)=1 is the recomended choice for nonlinear H
    346 PDAF       param_int(9): type_obs_init
    347 PDAF           Initialize observations before or after call to prepoststep_pdaf
    348 PDAF            0: Initialize observations before call to prepoststep_pdaf
    349 PDAF            1: Initialize observations after call to prepoststep_pdaf (default)
    350 PDAF     --- Floating point parameters (Array param_real) ---
    351 PDAF       param_real(1): forget
    352 PDAF           Forgetting factor (usually >0 and <=1), required
    353 PDAF     --- Further parameters ---
    354 PDAF       n_modeltasks: Number of parallel model integration tasks
    355 PDAF           >=1 for subtypes 0 and 1; not larger than total number of processors
    356 PDAF           =1 required for subtypes 10 and 11
    357 PDAF       screen: Control verbosity of PDAF
    358 PDAF           0: no outputs
    359 PDAF           1: basic output (default)
    360 PDAF           2: 1 plus timing output
    361 PDAF           3: 2 plus debug output
    362 PDAF     +++++++++ End of option overview for the LETKF ++++++++++
    363 }}}
    364 
    365 
    366 
     347== Option outputs using subtype=-1 - Global Ensemble Kalman filters ==
    367348
    368349=== ESTKF (filtertype=PDAF_DA_ESTKF=6) ===
    … …  
    422403
    423404
    424 
    425 === LESTKF (filtertype=PDAF_DA_LESTKF=7) ===
    426 
    427 {{{
    428 PDAF     Available options for LESTKF:
    429 PDAF     --- Sub-types (Parameter subtype) ---
    430 PDAF       0: Standard implementation with ensemble integration
    431 PDAF       10: Fixed error space basis
    432 PDAF       11: Fixed state covariance matrix
     405=== ETKF (filtertype=PDAF_DA_ETKF=4) ===
     406
     407{{{
     408PDAF     Available options for ETKF:
     409PDAF     --- Sub-types (Parameter subtype) ---
     410PDAF       0: full ensemble integration; apply T-matrix analogously to SEIK
     411PDAF       1: full ensemble integration; formulation cf. Hunt et al. (2007) without T matrix
     412PDAF       10: Fixed error space basis; analysis with T-matrix
     413PDAF       11: Fixed state covariance matrix; analysis with T-matrix
    433414PDAF     --- Integer parameters (Array param_int) ---
    434415PDAF       param_int(1): Dimension of state vector (>0), required
    … …  
    442423PDAF           Type of forgetting factor; optional
    443424PDAF            0: fixed forgetting factor (default)
    444 PDAF            1: adaptive forgetting factor for full domain
    445 PDAF            2: locally adaptive forgetting factor
     425PDAF            1: adaptive forgetting factor (experimental)
     426PDAF       param_int(6) type_trans
     427PDAF           Type of ensemble transformation matrix; optional
     428PDAF            0: deterministic transformation (default)
     429PDAF            2: use product of 0 with random orthonomal matrix with eigenvector (1,...,1)^T
     430PDAF       param_int(7): not used
     431PDAF       param_int(8): observe_ens
     432PDAF           Application of observation operator H, optional
     433PDAF            0: Apply H to ensemble mean to compute innovation
     434PDAF            1: Apply H to ensemble states; then compute innovation from their mean (default)
     435PDAF               param_int(8)=1 is the recomended choice for nonlinear H
     436PDAF       param_int(9): type_obs_init
     437PDAF           Initialize observations before or after call to prepoststep_pdaf
     438PDAF           0: Initialize observations before call to prepoststep_pdaf
     439PDAF           1: Initialize observations after call to prepoststep_pdaf (default)
     440PDAF     --- Floating point parameters (Array param_real) ---
     441PDAF       param_real(1): forget
     442PDAF           Forgetting factor (usually >0 and <=1), required
     443PDAF     --- Further parameters ---
     444PDAF       n_modeltasks: Number of parallel model integration tasks
     445PDAF           >=1 for subtypes 0 and 1; not larger than total number of processors
     446PDAF           =1 required for subtypes 10 and 11
     447PDAF       screen: Control verbosity of PDAF
     448PDAF           0: no outputs
     449PDAF           1: basic output (default)
     450PDAF           2: 1 plus timing output
     451PDAF           3: 2 plus debug output
     452PDAF     +++++++++ End of option overview for the ETKF ++++++++++
     453}}}
     454
     455
     456
     457=== EnKF (filtertype=PDAF_DA_ENKF=2) ===
     458
     459{{{
     460PDAF     Available options for EnKF:
     461PDAF     --- Sub-types (Parameter subtype) ---
     462PDAF       0: Full ensemble integration; analysis for 2*dim_obs>dim_ens
     463PDAF       1: Full ensemble integration; analysis for 2*dim_obs<=dim_ens
     464PDAF     --- Integer parameters (Array param_int) ---
     465PDAF       param_int(1): Dimension of state vector (>0), required
     466PDAF       param_int(2): Ensemble size (>0), required
     467PDAF       param_int(3): dim_lag
     468PDAF           Size of smoothing lag (>=0), optional
     469PDAF           0: no smoothing (default)
     470PDAF           >0: apply smoother up to specified lag
     471PDAF       param_int(4): rank_ana_enkf
     472PDAF           maximum rank for inversion of HPH^T, optional, default=0
     473PDAF            for =0, HPH is inverted by solving the representer equation
     474PDAF            allowed range is 0 to ensemble size - 1
     475PDAF       param_int(5): not used
     476PDAF       param_int(6): not used
     477PDAF       param_int(7): not used
     478PDAF       param_int(8): observe_ens
     479PDAF           Application of observation operator H, optional
     480PDAF            0: Apply H to ensemble mean to compute innovation
     481PDAF            1: Apply H to ensemble states; then compute innovation from their mean (default)
     482PDAF               param_int(8)=1 is the recomended choice for nonlinear H
     483PDAF       param_int(9): type_obs_init
     484PDAF           Initialize observations before or after call to prepoststep_pdaf
     485PDAF            0: Initialize observations before call to prepoststep_pdaf
     486PDAF            1: Initialize observations after call to prepoststep_pdaf (default)
     487PDAF     --- Floating point parameters (Array param_real) ---
     488PDAF       param_real(1): forget
     489PDAF           Forgetting factor (usually >0 and <=1), required
     490PDAF     --- Further parameters ---
     491PDAF       n_modeltasks: Number of parallel model integration tasks
     492PDAF           (>=1; not larger than total number of processors)
     493PDAF       screen: Control verbosity of PDAF
     494PDAF           0: no outputs
     495PDAF           1: basic output (default)
     496PDAF           2: 1 plus timing output
     497PDAF           3: 2 plus debug output
     498PDAF     +++++++++ End of option overview for the EnKF ++++++++++
     499}}}
     500
     501
     502=== SEIK (filtertype=PDAF_DA_SEIK=1) ===
     503
     504{{{
     505PDAF     Available options for SEIK:
     506PDAF     --- Sub-types (Parameter subtype) ---
     507PDAF       0: full ensemble integration; left-sided application of T
     508PDAF       1: full ensemble integration; right-sided application of T
     509PDAF       2: full ensemble integration; explicit ensemble transformation
     510PDAF       10: Fixed error space basis
     511PDAF       11: Fixed state covariance matrix
     512PDAF     --- Integer parameters (Array param_int) ---
     513PDAF       param_int(1): Dimension of state vector (>0), required
     514PDAF       param_int(2): Ensemble size (>0), required
     515PDAF       param_int(3): not used
     516PDAF       param_int(4): not used
     517PDAF       param_int(5) type_forget
     518PDAF           Type of forgetting factor; optional
     519PDAF            0: fixed forgetting factor (default)
     520PDAF            1: adaptive forgetting factor (experimental)
    446521PDAF       param_int(6) type_trans
    447522PDAF           Type of ensemble transformation matrix; optional
    … …  
    452527PDAF       param_int(7) type_sqrt
    453528PDAF           Type of transformation matrix square root; optional
     529PDAF            (Only relevant for subtype/=11)
    454530PDAF            0: symmetric square root (default)
    455531PDAF            1: Cholesky decomposition
    … …  
    468544PDAF     --- Further parameters ---
    469545PDAF       n_modeltasks: Number of parallel model integration tasks
    470 PDAF           >=1 for subtype 0; not larger than total number of processors
     546PDAF           >=1 for subtypes 0, 1 and 2; not larger than total number of processors
    471547PDAF           =1 required for subtypes 10 and 11
    472548PDAF       screen: Control verbosity of PDAF
    … …  
    475551PDAF           2: 1 plus timing output
    476552PDAF           3: 2 plus debug output
    477 PDAF     +++++++++ End of option overview for the LESTKF ++++++++++
    478 }}}
    479 
    480 
    481 
    482 
    483 === LEnKF (filtertype=PDAF_DA_LENKF=8) ===
    484 
    485 {{{
    486 PDAF     Available options for LEnKF:
    487 PDAF     --- Sub-types (Parameter subtype) ---
    488 PDAF       0: Standard EnKF analysis with covariance localization
    489 PDAF     --- Integer parameters (Array param_int) ---
    490 PDAF       param_int(1): Dimension of state vector (>0), required
    491 PDAF       param_int(2): Ensemble size (>0), required
    492 PDAF       param_int(3): not used
    493 PDAF       param_int(4): rank_ana_enkf
    494 PDAF           maximum rank for inversion of HPH^T, optional, default=0
    495 PDAF            for =0, HPH is inverted by solving the representer equation
    496 PDAF            allowed range is 0 to ensemble size - 1
    497 PDAF       param_int(5): not used
    498 PDAF       param_int(6): not used
    499 PDAF       param_int(7): not used
    500 PDAF       param_int(8): observe_ens
    501 PDAF           Application of observation operator H, optional
    502 PDAF            0: Apply H to ensemble mean to compute innovation
    503 PDAF            1: Apply H to ensemble states; then compute innovation from their mean (default)
    504 PDAF               param_int(8)=1 is the recomended choice for nonlinear H
    505 PDAF       param_int(9): type_obs_init
    506 PDAF           Initialize observations before or after call to prepoststep_pdaf
    507 PDAF            0: Initialize observations before call to prepoststep_pdaf
    508 PDAF            1: Initialize observations after call to prepoststep_pdaf (default)
    509 PDAF     --- Floating point parameters (Array param_real) ---
    510 PDAF       param_real(1): forget
    511 PDAF           Forgetting factor (usually >0 and <=1), required
    512 PDAF     --- Further parameters ---
    513 PDAF       n_modeltasks: Number of parallel model integration tasks
    514 PDAF           (>=1; not larger than total number of processors)
    515 PDAF       screen: Control verbosity of PDAF
    516 PDAF           0: no outputs
    517 PDAF           1: basic output (default)
    518 PDAF           2: 1 plus timing output
    519 PDAF           3: 2 plus debug output
    520 PDAF     +++++++++ End of option overview for the LEnKF ++++++++++
    521 }}}
    522 
    523 
    524 
     553PDAF     --- Internal parameter (defined inside PDAF) ---
     554PDAF       Nm1vsN: Normalization of covariance matrix; default: 1
     555PDAF           0: normalization with 1/(Ensemble size)
     556PDAF              (original SEIK, mainly for compatibility with older studies)
     557PDAF           1: normalization with 1/(Ensemble size - 1)
     558PDAF           (sample covariance matrix consistent with other EnKFs)
     559PDAF     +++++++++ End of option overview for the SEIK filter ++++++++++
     560}}}
     561
     562
     563
     564
     565
     566
     567
     568== Option outputs using subtype=-1 - Nonlinear DA Methods ==
    525569
    526570=== NETF (filtertype=PDAF_DA_NETF=9) ===
    … …  
    760804
    761805
    762 === ENSRF/EAKF (filtertype=PDAF_DA_ENSRF=13) ===
    763 
    764 {{{
    765 PDAF     Available options for ENSRF:
    766 PDAF     --- Sub-types (Parameter subtype) ---
    767 PDAF       0: ENSRF with serial observation processing (cf. Houtekamer/Mitchell, 2002)
    768 PDAF       1: EAKF/2-step local least squares filter (cf. Anderson, 2003)
    769 PDAF     --- Integer parameters (Array param_int) ---
    770 PDAF       param_int(1): Dimension of state vector (>0), required
    771 PDAF       param_int(2): Ensemble size (>0), required
    772 PDAF       param_int(3): not used
    773 PDAF       param_int(4): not used
    774 PDAF       param_int(5): not used
    775 PDAF       param_int(6): not used
    776 PDAF       param_int(7): not used
    777 PDAF       param_int(8): observe_ens
    778 PDAF           Application of observation operator H, optional
    779 PDAF            0: Apply H to ensemble mean to compute innovation
    780 PDAF            1: Apply H to ensemble states; then compute innovation from their mean (default)
    781 PDAF               param_int(8)=1 is the recomended choice for nonlinear H
    782 PDAF       param_int(9): type_obs_init
    783 PDAF           Initialize observations before or after call to prepoststep_pdaf
    784 PDAF            0: Initialize observations before call to prepoststep_pdaf
    785 PDAF            1: Initialize observations after call to prepoststep_pdaf (default)
    786 PDAF     --- Floating point parameters (Array param_real) ---
    787 PDAF       param_real(1): forget
    788 PDAF           Forgetting factor (usually >0 and <=1), required
    789 PDAF     --- Further parameters ---
    790 PDAF       n_modeltasks: Number of parallel model integration tasks
    791 PDAF           (>=1; not larger than total number of processors)
    792 PDAF       screen: Control verbosity of PDAF
    793 PDAF           0: no outputs
    794 PDAF           1: basic output (default)
    795 PDAF           2: 1 plus timing output
    796 PDAF           3: 2 plus debug output
    797 PDAF     +++++++++ End of option overview for the ENSRF ++++++++++
    798 }}}
    799 
    800 
    801 
     806
     807== Option outputs using subtype=-1 - Observation Geneeration ==
    802808
    803809=== GENOBS (filtertype=PDAF_DA_GENOBS=100) ===
    … …  
    830836
    831837
     838== Option outputs using subtype=-1 - 3D-Var Methods ==
    832839
    833840=== 3D-Var (filtertype=PDAF_DA_3DVAR=200) ===