| | 1 | = Implementation of the analysis step = |
| | 2 | |
| | 3 | {{{ |
| | 4 | #!html |
| | 5 | <div class="wiki-toc"> |
| | 6 | <h4>Implementation Guide</h4> |
| | 7 | <ol><li><a href="ImplementationGuide">Main page</a></li> |
| | 8 | <li><a href="AdaptParallelization">Adaptation of the parallelization</a></li> |
| | 9 | <li><a href="InitPdaf">Initialization of PDAF</a></li> |
| | 10 | <li><a href="ModifyModelforEnsembleIntegration">Modifications for ensemble integration</a></li> |
| | 11 | <li>Implementation of the analysis step</li> |
| | 12 | <ol> |
| | 13 | <li><a href="ImplementAnalysisestkf">Implementation for ESTKF</a></li> |
| | 14 | <li><a href="ImplementAnalysislestkf">Implementation for LESTKF</a></li> |
| | 15 | <li><a href="ImplementAnalysisetkf">Implementation for ETKF</a></li> |
| | 16 | <li><a href="ImplementAnalysisletkf">Implementation for LETKF</a></li> |
| | 17 | <li><a href="ImplementAnalysisseik">Implementation for SEIK</a></li> |
| | 18 | <li><a href="ImplementAnalysislseik">Implementation for LSEIK</a></li> |
| | 19 | <li><a href="ImplementAnalysisseek">Implementation for SEEK</a></li> |
| | 20 | <li><a href="ImplementAnalysisenkf">Implementation for EnKF</a></li> |
| | 21 | <li><a href="ImplementAnalysislenkf">Implementation for LEnKF</a></li> |
| | 22 | <li><a href="ImplementAnalysisnetf">Implementation for NETF</a></li> |
| | 23 | <li><a href="ImplementAnalysislnetf">Implementation for LNETF</a></li> |
| | 24 | <li><a href="ImplementAnalysislpf">Implementation for PF</a></li> |
| | 25 | </ol> |
| | 26 | <li><a href="AddingMemoryandTimingInformation">Memory and timing information</a></li> |
| | 27 | </ol> |
| | 28 | </div> |
| | 29 | }}} |
| | 30 | |
| | 31 | || This page describes the implementation of the analysis step without using PDAF-OMI. Please see the [wiki:ImplementationofAnalysisStep page on the analysis with OMI] for the more modern and efficient implementation variant using PDAF-OMI. || |
| | 32 | |
| | 33 | The implementation of the analysis step requires that the user-supplied routines whose names are given in the call to the filter-specific put-state routine are implemented with full functionality. |
| | 34 | |
| | 35 | Several of the routines are general and hence usable for all filter algorithms. There are some differences between the available local and global filters. The routines are described in separate pages for each filter algorithm. The general routines will be repeated, but they are marked as general. |
| | 36 | |
| | 37 | Filter-specific implementations for global filters: |
| | 38 | * [ImplementAnalysisestkf Implementation of the analysis for the ESTKF] |
| | 39 | * [ImplementAnalysisetkf Implementation of the analysis for the ETKF] |
| | 40 | * [ImplementAnalysisseik Implementation of the analysis for the SEIK filter] |
| | 41 | * [ImplementAnalysisseek Implementation of the analysis for the SEEK filter] |
| | 42 | * [ImplementAnalysisenkf Implementation of the analysis for the EnKF] |
| | 43 | * [ImplementAnalysisnetf Implementation of the analysis for the NETF] |
| | 44 | * [ImplementAnalysispf Implementation of the analysis for the PF] |
| | 45 | |
| | 46 | Filter-specific implementations for domain-local filters: |
| | 47 | * [ImplementAnalysislestkf Implementation of the analysis for the LESTKF] |
| | 48 | * [ImplementAnalysisletkf Implementation of the analysis for the LETKF] |
| | 49 | * [ImplementAnalysislseik Implementation of the analysis for the LSEIK filter] |
| | 50 | * [ImplementAnalysislnetf Implementation of the analysis for the LNETF] |
| | 51 | |
| | 52 | Filter-specific implementation for filter with covariance localization: |
| | 53 | * [ImplementAnalysislenkf Implementation of the analysis for the LEnKF] |
| | 54 | |
| | 55 | (The filters ESTKF and LESTKF have been added with PDAF version 1.8; the filters LEnKF, NETF and LNETF have been addid with PDAF version 1.12) |
| | 56 | |
| | 57 | == Testing the full implementation == |
| | 58 | |
| | 59 | After the implementation of the user-supplied routines for the analysis step, the data assimilation system is fully implemented. In order to test the program one has to compile it without setting 'PDAF_NO_UPDATE'. |
| | 60 | |
| | 61 | To speed up the testing it can be useful to comment out the time stepping part. In that case only the model initialization and the assimilation cycle would be executed, without the usually costly time integration of the ensemble. One should then check if all user-supplied routines perform as they should. |
| | 62 | |
| | 63 | |
