Changes between Version 6 and Version 7 of InitPdaf

Timestamp:

Aug 24, 2010, 3:06:27 PM (14 years ago)

Author:

lnerger

Comment:

--

InitPdaf

@@ -58 +58 @@
  * `status_pdaf`: An integer used as status flag of PDAF. If `status_pdaf` is zero upon exit from `PDAF_init` the initialization was successful. An error occurred for non-zero values. (The error codes are documented in the routine PDAF_init.) 
 
+It is recommended that the value of `status_pdaf` is checked in the program after PDAF_init is executed. Only if its value is 0 the initialization was successful.
+
 == Other variables for the assimilation ==
 
@@ -98 +100 @@
 The only mode-based filter supplied with PDAF is currenly the SEEK filter. For this filter the initialization bases on the decomposition of the state error covariance matrix in the form '''P''' = '''VUV^T^'''. According to this decomposition, the array `ens_p` has to be initialized to hold the modes from matrix '''V''' and `Uinv` holds the inverse of matrix '''U'''. In addition `state_p` has to be initialized with the initial state estimate.  If a model with domain decomposition is used, the part of all modes for the local sub-domain of the MPI process and the corresponding part of the state vector has to be initialized.
 
+== Testing the PDAF initialization ==
+
+The PDAF initialization can be tested by compiling the program and executing it. At this stage it will not be meaningful to perform an actual integration. However, one can test if the initialization in PDAF_init is sucessful and if the ensemble array is correctly initialized.
+
+Standard output from PDAF_init should look like the following:
+{{{
+  PDAF: Initialize SEIK Filter
+
+         SEIK configuration
+               filter sub-type = 0
+                 --> Standard SEIK
+                 --> Transform ensemble with deterministic Omega
+                 --> Use fixed forgetting factor: 1.00
+
+         Parallelization - Filter on model PEs:
+                      Total number of PEs:   4
+          Number of parallel model tasks:   1
+                           PEs for Filter:   4
+     # PEs per ensemble task and local ensemble sizes: 
+           Task   1
+           #PEs   4
+              N 300
+}}}
+
+The correctness of the ensemble initialization in `U_ensemble_init` should be checked by the user.