Changes between Initial Version and Version 1 of PDAFomi_observation_diagnostics


Ignore:
Timestamp:
Mar 26, 2025, 5:07:07 PM (7 days ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • PDAFomi_observation_diagnostics

    v1 v1  
     1= PDAF-OMI Observation Diagnostics =
     2
     3{{{
     4#!html
     5<div class="wiki-toc">
     6<h4>PDAF-OMI Guide</h4>
     7<ol><li><a href="PDAF_OMI_Overview">Overview</a></li>
     8<li><a href="OMI_Callback_obs_pdafomi">callback_obs_pdafomi.F90</a></li>
     9<li><a href="OMI_observation_modules">Observation Modules</a></li>
     10<li><a href="OMI_observation_operators">Observation operators</a></li>
     11<li><a href="OMI_error_checking">Checking error status</a></li>
     12<li><a href="OMI_debugging">Debugging functionality</a></li>
     13<li><a href="OMI_ImplementationofAnalysisStep">Implementing the analysis step with OMI</a></li>
     14<ol>
     15<li> <a href="ImplementFilterAnalysisOverview"> General overview for ensemble filters</a></li>
     16<ol>
     17<li><a href="ImplementAnalysisGlobal">Implementation for Global Filters</a></li>
     18<li><a href="ImplementAnalysisLocal">Implementation for Local Filters</a></li>
     19<li><a href="ImplementAnalysislenkfOmi">Implementation for LEnKF</a></li>
     20</ol>
     21<li> <a href="Implement3DVarAnalysisOverview"> General overview for 3D-Var methods</a></li>
     22<ol>
     23<li><a href="ImplementAnalysis_3DVar">Implementation for 3D-Var</a></li>
     24<li><a href="ImplementAnalysis_3DEnVar">Implementation for 3D Ensemble Var</a></li>
     25<li><a href="ImplementAnalysis_Hyb3DVar">Implementation for Hybrid 3D-Var</a></li>
     26</ol>
     27</ol>
     28<li><a href="OMI_nondiagonal_observation_error_covariance_matrices">Using nondiagonal R-matrices</a></li>
     29<li><a href="Porting_to_OMI">Porting an existing implemention to OMI</a></li>
     30<li><a href="PDAFomi_additional_functionality">Additional OMI Functionality</a></li>
     31<li>Observation diagnostics</li>
     32</ol>
     33</div>
     34}}}
     35
     36[[PageOutline(2-3,Contents of this page)]]
     37
     38The PDAF-OMI observation diagnostics module provides functionality to obtain statistics about the differences between observations and the observed model state. In addition, there are routines that provide the user access to the observations and observed quantities like the observed ensemble mean state.
     39
     40Here, we describe the functionalities of the observation diagnostics routines.
     41
     42A common place to call the PDAFomi_diag diagnostics routines is in `prepoststep_pdaf`, which is the usual place to also analyze the ensemble. Usually, PDAF initializes the observations after `prepoststep_pdaf` was executed after the forecast. To be able to compare the observations and the forecast ensemble, one has to switch the place at which observations are initialized. This is done with
     43{{{
     44  CALL PDAF_set_iparam(9, 0)
     45}}}
     46which can be called in `init_pdaf` after the initialization of PDAF with `PDAF_init`.
     47
     48The routines for observation diagnostics can be organized in four groups
     49 * Activation of observation diagnostics
     50   * [PDAFomi_observation_diagnostics#PDAFomi_set_obs_diag PDAFomi_set_obs_diag]
     51 * Statistics
     52   * PDAFomi_diag_rmsd
     53   * PDAFomi_diag_stats
     54 * Access to observation dimensions
     55   * PDAFomi_diag_nobs
     56   * PDAFomi_diag_dimobs
     57 * Acces to observation arrays
     58   * PDAFomi_diag_get_obs
     59   * PDAFomi_diag_get_HXmean
     60   * PDAFomi_diag_get_HX
     61   * PDAFomi_diag_get_ivar
     62
     63== Activation of observation diagnostics ==
     64
     65=== PDAFomi_set_obs_diag ===
     66
     67The routine is used to activate or deactivate the observation diagnostics.
     68
     69The routine can be called by all processes, but it is sufficient to call it for those processes that handle observations, which usually are the filter processes. A common place is to call the routine in `init_pdaf` afer the initialization of PDAF in `PDAF_init`.
     70
     71The interface is:
     72{{{
     73  SUBROUTINE PDAFomi_set_obs_diag(diag)
     74
     75    INTEGER, INTENT(in) :: diag   ! Value for observation diagnostics mode
     76                                  ! >0 activates observation diagnostics
     77}}}
     78
     79
     80
     81