Changes between Version 3 and Version 4 of PDAFomi_gather_obs

Timestamp:

Jun 3, 2024, 4:43:31 PM (6 months ago)

Author:

lnerger

Comment:

--

PDAFomi_gather_obs

@@ -10 +10 @@
 {{{
   SUBROUTINE PDAFomi_gather_obs(thisobs, dim_obs_p, obs_p, ivar_obs_p, ocoord_p, &
-       ncoord, lradius, dim_obs)
+       ncoord, cradius, dim_obs)
 
     TYPE(obs_f), INTENT(inout) :: thisobs   ! Data type with full observation information
@@ -18 +18 @@
     REAL, INTENT(in) :: ocoord_p(:,:)       ! Array of process-local observation coordinates
     INTEGER, INTENT(in) :: ncoord           ! Number of rows (first index) of coordinate array
-    REAL, INTENT(in) :: lradius             ! Typical localization radius (the maximum radius used in this process domain) 
+    REAL, INTENT(in) :: cradius             ! Typical localization radius (the maximum radius used in this process domain) 
     INTEGER, INTENT(out) :: dim_obs         ! Full number of observations
 
@@ -25 +25 @@
 '''Notes:'''
 * The only direect output is `dim_obs`, which is then returned to `init_dim_obs_pdafomi` in `callback_obs_pdafomi.F90`.
-* The value of `lradius` can be some maximum value. It is only utilized if domain-limited observations are activated ([wiki:OMI_use_global_obs see page on using the option `use_global_obs`])
+* The value of `cradius` can be some maximum value. It is only utilized if domain-limited observations are activated ([wiki:OMI_use_global_obs see page on using the option `use_global_obs`])
+* Also in case of non-isotropic localization `cradius` is a single value. It should be set of the largest radius used in the parallel process domain.