| | 1 | = PDAF_iau_init_inc = |
| | 2 | |
| | 3 | This page documents the routine `PDAF_iau_init_inc` of PDAF, which was introduced with PDAF V3.0. |
| | 4 | |
| | 5 | This routine can be used by a user to fill the process-local increment array. |
| | 6 | A common use is when the IAU should be applied from the initial time of a run, for example if the increment was computed in a previous assimilation run and then stored for restarting. Since PDAF can only compute an increment in an analysis step, the user needs to provide the increment. |
| | 7 | |
| | 8 | The routine is usually called in init_pdaf after the initialization of IAU with `PDAF_iau_init`. It has to be called by all processes that are model processes and one needs to provide the task-local ensemble (i.e. with local ensmeble size `dim_ens_l=1` for the fully parallel mode, and usually `dim_ens_l>1` for the flexible parallelization mode) |
| | 9 | |
| | 10 | The interface is the following: |
| | 11 | {{{ |
| | 12 | SUBROUTINE PDAF_iau_init_inc(dim_p, dim_ens_l, ens_inc, flag) |
| | 13 | }}} |
| | 14 | with the following arguments: |
| | 15 | {{{ |
| | 16 | INTEGER, INTENT(in) :: dim_p !< PE-local dimension of model state |
| | 17 | INTEGER, INTENT(in) :: dim_ens_l !< Task-local size of ensemble |
| | 18 | REAL, INTENT(in) :: ens_inc(dim_p, dim_ens_l) !< PE-local increment ensemble |
| | 19 | INTEGER, INTENT(out) :: flag !< Status flag |
| | 20 | }}} |
| | 21 | |
