Changes between Version 1 and Version 2 of PDAF_diag_ensmean

Timestamp:

Mar 26, 2025, 11:09:04 AM (7 days ago)

Author:

lnerger

Comment:

--

PDAF_diag_ensmean

@@ -7 +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.
 
-The interface is the following: 
+The interface is: 
 {{{
-SUBROUTINE PDAF_diag_ensmean(dim, dim_ens, state, ens, status)
+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)
 }}}
-with the following arguments:
-{{{
-  INTEGER, INTENT(in) :: dim               ! state dimension
-  INTEGER, INTENT(in) :: dim_ens           ! Ensemble size
-  REAL, INTENT(inout) :: state(dim)        ! ensemble mean state vector
-  REAL, INTENT(in)    :: ens(dim, 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`.