| 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. |