Changes between Version 1 and Version 2 of PDAFomi_init_dim_obs_l


Ignore:
Timestamp:
Apr 23, 2024, 3:39:56 PM (8 months ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • PDAFomi_init_dim_obs_l

    v1 v2  
    1616    REAL, INTENT(in) :: coords_l(:)          ! Coordinates of current local analysis domain
    1717    INTEGER, INTENT(in) :: locweight         ! Type of localization function
    18     REAL, INTENT(in) :: lradius              ! Localization radius
     18    REAL, INTENT(in) :: cradius              ! Localization cut-off radius
    1919    REAL, INTENT(in) :: sradius              ! Support radius of localization function
    2020    INTEGER, INTENT(inout) :: dim_obs_l      ! Local dimension of current observation vector
     
    2525 * The coordinate vector `coords_l` has be set consistently with the observation coordinates specified in `init_dim_obs_OBSTYPE`.
    2626 * 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])
    27  * `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`.
     27 * `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`.
     28
     29== Settings for non-isotropic localization ==
     30
     31Introduced 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
     32{{{
     33    REAL, INTENT(in) :: cradius(:)              ! Localization radius
     34    REAL, INTENT(in) :: sradius(:)              ! Support radius of localization function
     35
     36}}}
     37Here, the size `cradius` and `sradius` has to be equal to `thisobs%ncoords`.