| | 1 | = PDAF_diag_stddev_nompi = |
| | 2 | |
| | 3 | This page documents the routine `PDAF_diag_stddev_nompi` of PDAF, which was introduced with PDAF V3.0. |
| | 4 | |
| | 5 | This routine computes mean ensemble standard deviation. |
| | 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]. |
| | 8 | |
| | 9 | The interface is the following: |
| | 10 | {{{ |
| | 11 | SUBROUTINE PDAF_diag_stddev_nompi(dim, dim_ens, & |
| | 12 | state, ens, stddev, do_mean, status) |
| | 13 | }}} |
| | 14 | with the following arguments: |
| | 15 | {{{ |
| | 16 | INTEGER, INTENT(in) :: dim !< state dimension |
| | 17 | INTEGER, INTENT(in) :: dim_ens !< Ensemble size |
| | 18 | REAL, INTENT(inout) :: state(dim) !< State vector |
| | 19 | REAL, INTENT(in) :: ens(dim, dim_ens) !< State ensemble |
| | 20 | REAL, INTENT(out) :: stddev !< Standard deviation of ensemble |
| | 21 | INTEGER, INTENT(in) :: do_mean !< Whether to compute ensemble mean |
| | 22 | INTEGER, INTENT(out) :: status !< Status flag (0=success) |
| | 23 | }}} |
| | 24 | |
| | 25 | **Note:** |
| | 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. |