Changes between Version 4 and Version 5 of OMI_debugging_PDAF3

Timestamp:

May 9, 2026, 2:36:44 PM (6 days ago)

Author:

lnerger

Comment:

--

OMI_debugging_PDAF3

@@ -86 +86 @@
      REAL :: inno_omit_ivar=1.0e-12       !< Value of inverse variance to omit observation
                                           !<     (should be much larger than actual observation error variance)
+     CHARACTER(len=20) :: name=""         !< Name of observation type (optional)
 
      ! ----  The following variables are set in the routine PDAFomi_gather_obs ---
@@ -99 +100 @@
      INTEGER, ALLOCATABLE :: id_obs_f_lim(:) !< Indices of domain-relevant full obs. in global vector of obs.
   END TYPE obs_f
+     INTEGER, ALLOCATABLE :: idx_sort(:)  !< Re-sorted indices for search_type>0
+
+     ! ----  The following variables are used internally for observation diagnostics ---
+     REAL, ALLOCATABLE :: obs_diag_p(:)        !< Full observed field
+     REAL, ALLOCATABLE :: ocoord_diag_p(:,:)   !< Coordinates of full observation vector
+     REAL, ALLOCATABLE :: ivar_obs_diag_p(:)   !< Inverse variance of full observations
+     REAL, ALLOCATABLE :: HX_diag_p(:,:)       !< Ensemble of observed model states
+     REAL, ALLOCATABLE :: HXmean_diag_p(:)     !< Observed ensemble meanstate
+
+     ! ----  The following variables can be set in the routine PDAFomi_set_localize_covar ---
+     INTEGER :: nradii                    !< Length of CRADIUS and SRADIUS
+     INTEGER, ALLOCATABLE :: locweight(:) !< localization weight type (single value or vector for factorized localization)
+     REAL, ALLOCATABLE :: cradius(:)      !< Localization cut-off radius (single value or vector)
+     REAL, ALLOCATABLE :: sradius(:)      !< support radius for localization function (single value or vector)
+
+     ! ----  Other internal variables ---
+     INTEGER :: locweight_v               !< Type of localization function in vertical direction (for disttype>=10)
 }}}
 
@@ -110 +128 @@
      REAL, ALLOCATABLE :: sradius_l(:)    !< directional support radii of local observations
      REAL, ALLOCATABLE :: ivar_obs_l(:)   !< Inverse variance of local observations
+     REAL, ALLOCATABLE :: dist_l_v(:)     !< Vertical distances of local observations
      INTEGER :: locweight                 !< Specify localization function
+     INTEGER :: locweight_v=0             !< Specify localization function for vertical direction
+     INTEGER :: nradii                    !< Length of CRADIUS and SRADIUS
      REAL, ALLOCATABLE :: cradius(:)      !< Localization cut-off radius (single value or vector)
      REAL, ALLOCATABLE :: sradius(:)      !< support radius for localization function (single value or vector)