Changes between Version 3 and Version 4 of PDAFomi_set_searchtype

Timestamp:

Apr 17, 2026, 12:14:41 PM (5 weeks ago)

Author:

lnerger

Comment:

--

PDAFomi_set_searchtype

@@ -17 +17 @@
            ! 0: Search routine of PDAF 3.0
            ! 1: Optimized search code
-           ! 2: Optimized search code with only one search loop
+           ! 2: Optimized search code with only one search loop (default)
            ! 11: Search type using sorted observations
            ! 12: Search type using sorted observations with only one search loop
@@ -24 +24 @@
            !       and only one search loop
            !<  (2, 12 and 22 need more memory as arrays are allocated with full observation dimension)
-    INTEGER, INTENT(in) :: sortdir         !< Sorted doordinate direction
+    INTEGER, INTENT(in) :: sortdir         !< Sorted doordinate direction (default=2)
 }}}
 
@@ -33 +33 @@
  * In experiments, we saw that the search performance of the different types depends on the distribution of the observations. Generally type 0 is slowest, but the performance of the other cases is case-specific.
  * search types 11-22 can give slightly different assimilation results compared 0-2, because the sorting of observations change the order of computations. 
- * Usually, `sortdir` should be the highest index of coordinate direction along which a large variability of coordinates exist, e.g. the latitude in geographic coordinates
+ * Usually, `sortdir` should be the highest index of coordinate direction along which a large variability of coordinates exist, e.g. the latitude in geographic coordinates; However, depending on the specific case a different direction might lead to a better performance of the search.
  * `sortdir` cannot be a direction with periodicity (e.g. longitude in global Earth models is periodic and not a proper choice)