Changes between Version 1 and Version 2 of PDAFomi_init_dim_obs_l

Timestamp:

Apr 23, 2024, 3:39:56 PM (7 months ago)

Author:

lnerger

Comment:

--

PDAFomi_init_dim_obs_l

@@ -16 +16 @@
     REAL, INTENT(in) :: coords_l(:)          ! Coordinates of current local analysis domain
     INTEGER, INTENT(in) :: locweight         ! Type of localization function
-    REAL, INTENT(in) :: lradius              ! Localization radius
+    REAL, INTENT(in) :: cradius              ! Localization cut-off radius
     REAL, INTENT(in) :: sradius              ! Support radius of localization function
     INTEGER, INTENT(inout) :: dim_obs_l      ! Local dimension of current observation vector
@@ -25 +25 @@
  * The coordinate vector `coords_l` has be set consistently with the observation coordinates specified in `init_dim_obs_OBSTYPE`.
  * a convenient place to initialize `coords_l` is in `init_dim_l`. The templates and tutorial codes show this approach. (see the [wiki:ImplementAnalysisLocal page on implementing the local analysis step])
- * `locweight`, `lradius` and `sradius` can be set to the same values for all local analysis domains. However, the interface allows the user to specify different localization radii and weight functions for each single local analysis domain depending on the index `domain_p` or the coordinates `coords_l`.
+ * `locweight`, `cradius` and `sradius` can be set to the same values for all local analysis domains. However, the interface allows the user to specify different localization radii and weight functions for each single local analysis domain depending on the index `domain_p` or the coordinates `coords_l`.
+
+== Settings for non-isotropic localization ==
+
+Introduced by PDAF V2.2 a non-isotropic localization is possible for this use vectors for `cradius` and `lradius`. This in the interface of `PDAFomi_init_dim_obs_l` one can use
+{{{
+    REAL, INTENT(in) :: cradius(:)              ! Localization radius
+    REAL, INTENT(in) :: sradius(:)              ! Support radius of localization function
+
+}}}
+Here, the size `cradius` and `sradius` has to be equal to `thisobs%ncoords`.