= PDAF_diag_ensmean = 

This page documents the routine `PDAF_diag_ensmean` of PDAF, which was introduced with PDAF V3.0.

This routine computes the ensemble mean state vector.

The routine can be called in the pre/poststep routine of PDAF both before and after the analysis step to compute the ensemble statistics.

The interface is: 
{{{
SUBROUTINE PDAF_diag_ensmean(dim_p, dim_ens, state_p, ens_p, status)

  INTEGER, INTENT(in) :: dim_p               ! PE-local state dimension
  INTEGER, INTENT(in) :: dim_ens             ! Ensemble size
  REAL, INTENT(inout) :: state_p(dim_p)      ! ensemble mean state vector
  REAL, INTENT(in)    :: ens(dim_p, dim_ens) ! State ensemble
  INTEGER, INTENT(out) :: status             ! Status flag (0=success)
}}}

**Notes:**
 * The routine compute the ensemble mean for the given ensemble array `ens_p`. This works for both a parallelized model and without parallelization. 
 * The routine does not perform any MPI operations and can also be used if PDAF was not initialized by calling `PDAF_init`.