Changes between Version 4 and Version 5 of PDAF_diag_stddev_nompi
- Timestamp:
- Mar 26, 2025, 2:13:46 PM (6 days ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
PDAF_diag_stddev_nompi
v4 v5 3 3 This page documents the routine `PDAF_diag_stddev_nompi` of PDAF, which was introduced with PDAF V3.0. 4 4 5 This routine computes mean ensemble standard deviation .5 This routine computes mean ensemble standard deviation as the square root of the spatial mean variance. 6 6 7 The routine can be called in the pre/poststep routine of PDAF both before and after the analysis step to compute the ensemble statistics. This variant only compute the standard deviation for the provided ensemble array `ens`. Thus, with a domain-decomposed parallel model, the value of of the standard deviation doe snot consider the full domain. The variant taking into account parallelization is [wiki:PDAF_diag_stddev]. 7 The routine can be called in the pre/poststep routine of PDAF both before and after the analysis step to compute the ensemble statistics. 8 9 This variant only computes the standard deviation for the provided ensemble array `ens`. Thus, with a domain-decomposed parallel model, the value of of the standard deviation does not consider the full domain. The variant taking into account parallelization is [wiki:PDAF_diag_stddev]. 8 10 9 11 The interface is the following: … … 11 13 SUBROUTINE PDAF_diag_stddev_nompi(dim, dim_ens, & 12 14 state, ens, stddev, do_mean, status) 13 }}} 14 with the following arguments: 15 {{{ 15 16 16 INTEGER, INTENT(in) :: dim !< state dimension 17 17 INTEGER, INTENT(in) :: dim_ens !< Ensemble size … … 26 26 * The ensemble standard deviation is a common measure of the estimate model root mean square error. It is typically computed in `prepoststep_pdaf` to monitor the assimilation process. 27 27 * We recommend to use this routine only in cases where the MPI parallelization is not initialized, e.g. in separate post-analysis programs. In assimilation programs, where the parallelization is initialized for PDAF, we recommend to use the routine [wiki:PDAF_diag_stddev] for better overall compatibility. 28 * The option `do_mean` exists mainly for performance reasons. If `state_p` contains is ensmeble mean state it does not need to be computed again. 28 * The option `do_mean` exists mainly for performance reasons. If `state_p` contains is ensemble mean state it does not need to be computed again. 29 * The routine does not perform any MPI operations and can also be used if PDAF was not initialized by calling `PDAF_init`.