Changes between Version 2 and Version 3 of PDAFomi_diag_stats


Ignore:
Timestamp:
Mar 26, 2025, 4:43:48 PM (6 days ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • PDAFomi_diag_stats

    v2 v3  
    11= PDAFomi_diag_stats =
    22
    3 This page documents the routine `PDAFomi_diag_stats` of PDAF, which was introduced with PDAF V3.0. This is part of the PDAF-OMI observation diagnostics module.
     3This page documents the routine `PDAFomi_diag_stats` of PDAF, which was introduced with PDAF V3.0. This is part of the [wiki:PDAFomi_observation_diagnostics PDAF-OMI observation diagnostics module].
    44
    5 The routine returns
    6 a pointer to an array of a selection of 6 statistics comparing the observations and the observed ensemble mean for each observation type.
     5The routine returns a pointer to an array of a selection of 6 statistics comparing the observations and the observed ensemble mean for each observation type. The statistics can, for example, be used to plot a Taylor diagram.
    76
    8 Usually all PDAFomi_diag routines are called in `prepoststep_pdaf` where the observation information can be retrieved and analyzed. These statistics can, for exmaple, be used to plot a Taylor diagram.
     7Usually all PDAFomi_diag routines are called in `prepoststep_pdaf` where the observation information can be retrieved and analyzed.
    98
    109The interface is the following:
    1110{{{
    1211  SUBROUTINE PDAFomi_diag_stats(nobs, obsstats_ptr, verbose)
    13 }}}
    14 with the following arguments:
    15 {{{
     12
    1613    INTEGER, INTENT(inout) :: nobs                     ! Number of observation types
    1714    REAL, POINTER, INTENT(inout) :: obsstats_ptr(:,:)  ! Array of observation statistics
    1815          ! Included statistics are:
    1916          !  (1,:) correlations between observation and observed ensemble mean
    20           !  (2,:) centered RMSD between observation and observed ensemble mean
    21           !  (3,:) mean bias for observation minus observed ensemble mean
     17          !  (2,:) centered RMS difference between observation and observed ensemble mean
     18          !  (3,:) mean bias (observation minus observed ensemble mean)
    2219          !  (4,:) mean absolute difference between observation and observed ensemble mean
    2320          !  (5,:) variance of observations
    2421          !  (6,:) variance of observed ensemble mean
    25     INTEGER, INTENT(in) :: verbose                     ! Verbosity flag
     22    INTEGER, INTENT(in) :: verbose                     ! Verbosity flag, >0 to write output
    2623}}}
    2724
    2825**Note:**
    2926 * The computed statistics are for the global model domain. Thus, in case of a parallelized model, all process sub-domains are taken into account and calling `PDAFomi_diag_stats` will return the same value for all processes.
    30  * In Fortran user code the pointer should be declared in the form[[BR]] `REAL, POINTER :: obsstats_ptr(:)`[[BR]] It does not need to be allocated.. The target array has the size `(6, nobs)`.
     27 * In Fortran user code the pointer should be declared in the form[[BR]] `REAL, POINTER :: obsstats_ptr(:)`[[BR]] It does not need to be allocated. The target array has the size `(6, nobs)`.
    3128 * If the observation diagnostics have not be activated by using [wiki:PDAFomi_set_obs_diag] the pointer array will not be set and `nobs=0` is. One can check this value before assessing the pointer array
    3229 * The routine returns the centered RMSD as displayed in Taylor diagrams. The non-centered RMSD can be computed using [wiki:PDAFomi_diag_obs_rmsd].