Changes between Version 5 and Version 6 of ImplementAnalysis_3DVar

Timestamp:

Dec 9, 2021, 7:59:51 AM (2 years ago)

Author:

lnerger

Comment:

--

ImplementAnalysis_3DVar

@@ -17 +17 @@
 <li><a href="ImplementAnalysis_3DEnVar">Implementation for 3D-EnVar</a></li>
 <li><a href="ImplementAnalysis_Hyb3DVar">Implementation for hybrid 3D-Var</a></li>
-<li><a href="ImplementAnalysislenkfOmi">Implementation for LEnKF</a></li>
 <li><a href="PDAF_OMI_Overview">PDAF-OMI Overview</a></li>
 </ol>
@@ -186 +185 @@
 == Execution order of user-supplied routines ==
 
-The user-supplied routines are essentially executed in the order they are listed in the interface to `PDAFomi_assimilate_global`. The order can be important as some routines can perform preparatory work for later routines. For example, `U_init_dim_obs_pdafomi` prepares an index array that provides the information for executing the observation operator in `U_obs_op_pdafomi`. How this information is initialized is described in the documentation of OMI.
+The user-supplied routines are essentially executed in the order they are listed in the interface to `PDAFomi_assimilate_3dvar`. The order can be important as some routines can perform preparatory work for later routines. For example, `U_init_dim_obs_pdafomi` prepares an index array that provides the information for executing the observation operator in `U_obs_op_pdafomi`. How this information is initialized is described in the documentation of OMI.
 
 Before the analysis step is called the following routine is executed:
@@ -195 +194 @@
  1. [#U_init_dim_obs_pdafomicallback_obs_pdafomi.F90 U_init_dim_obs_pdafomi]
  1. [#U_obs_op_pdafomicallback_obs_pdafomi.F90 U_obs_op_pdafomi] (multiple calls, one for each ensemble member)
+
+Inside the analysis step the interative optimization is computed. This involves the repeated call of the routines:
+ 1. [#U_cvtcvt_pdaf.F90 U_cvt] 
+ 1. [#U_obs_op_linpdafomicallback_obs_pdafomi.F90 U_obs_op_lin_pdafomi]
+ 1. [#U_obs_op_adjpdafomicallback_obs_pdafomi.F90 U_obs_op_adj_pdafomi]
+ 1. [#U_cvt_adjcvt_adj_pdaf.F90 U_cvt_adj] 
+
+After the iterative optimization the following routines are executes to complte the analysis step:
+ 1. [#U_cvtcvt_pdaf.F90 U_cvt] (Call to the control vector transform to compute the final state vector increment
  1. [#U_prepoststepprepoststep_ens_pdaf.F90 U_prepoststep] (Call to act on the analysis ensemble, called with (positive) value of the time step)