Changes between Version 32 and Version 33 of ImplementAnalysisseik


Ignore:
Timestamp:
Sep 15, 2011, 2:28:28 PM (13 years ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • ImplementAnalysisseik

    v32 v33  
    6262
    6363
    64 === `U_collect_state` (collect_state.F90) ===
     64=== `U_collect_state` (collect_state_pdaf.F90) ===
    6565
    6666This routine is independent of the filter algorithm used.
     
    6868
    6969
    70 === `U_init_dim_obs` (init_dim_obs.F90) ===
     70=== `U_init_dim_obs` (init_dim_obs_pdaf.F90) ===
    7171
    7272This routine is used by all global filter algorithms (SEEK, SEIK, EnKF, ETKF).
     
    8686
    8787
    88 === `U_obs_op` (obs_op.F90) ===
     88=== `U_obs_op` (obs_op_pdaf.F90) ===
    8989
    9090This routine is used by all global filter algorithms (SEEK, SEIK, EnKF, ETKF).
     
    109109
    110110
    111 === `U_init_obs` (init_obs.F90) ===
     111=== `U_init_obs` (init_obs_pdaf.F90) ===
    112112
    113113This routine is used by all global filter algorithms (SEEK, SEIK, EnKF, ETKF).
     
    128128
    129129
    130 === `U_prepoststep` (prepoststep_seik.F90) ===
     130=== `U_prepoststep` (prepoststep_seik_pdaf.F90) ===
    131131
    132132The routine has already been described on the [ModifyModelforEnsembleIntegration#U_prepoststepprepoststep_seik.F90 page on modifying the model code for the ensemble integration]. For completeness, the description is repeated:
     
    163163
    164164
    165 === `U_prodRinvA` (prodrinva.F90) ===
     165=== `U_prodRinvA` (prodrinva_pdaf.F90) ===
    166166
    167167This routine is used by all filter algorithms that use the inverse of the observation error covariance matrix (SEEK, SEIK, and ETKF).
     
    188188 * The interface has a difference for SEIK and ETKF: For ETKF the third argument is the ensemble size (`dim_ens`), while for SEIK it is the rank (`rank`) of the covariance matrix (usually ensemble size minus one). In addition, the second dimension of `A_p` and `C_p` has size `dim_ens` for ETKF, while it is `rank` for the SEIK filter.  (Practically, one can usually ignore this difference as the fourth argument of the interface can be named arbitrarily in the routine.)
    189189
    190 === `U_init_obsvar` (init_obsvar.F90) ===
     190=== `U_init_obsvar` (init_obsvar_pdaf.F90) ===
    191191
    192192This routine is used by the global filter algorithms SEIK and  ETKF as well as the local filters LSEIK and LETKF. The routine is only called if the adaptive forgetting factor is used (`type_forget=1` in the example impementation).