Consolidation with CDI-PIO (develop)

Generally, this is !11 (merged) rebased onto the develop branch.

This is the first develop-based version accepted by the Buildbot of ICON: https://buildbot.dkrz.de/grid?branch=m300488/cdi-pio-2-integration

src/pio_dist_grid.c

@@ -450,6 +450,20 @@ cdiPioDistGridInit(grid_t *gridptr, int gridtype, int size, int xsize, int ysize
       }
     else
       {
+        /* check if hints are in range */
+        static const char dimName[2] = { 'X', 'Y' }, caller[] = "cdiPioDistGridCreate";
+        for (size_t aRank = 0; aRank < saMaxRank; ++aRank)
+          {
+            int32_t dimUB = sub_arrays[cdiPioGDsaArea].rect[aRank].first + sub_arrays[cdiPioGDsaArea].rect[aRank].size;
+            if (xy_decomposition[aRank][0] < 0 || xy_decomposition[aRank][0] > dimUB)
+              Errorc("decomposition hint for %c start out "
+                     "of range [0,%d]\n",
+                     dimName[aRank], (int) dimUB - 1);
+            if (xy_decomposition[aRank][1] < 0 || (xy_decomposition[aRank][0] + xy_decomposition[aRank][1] > dimUB))
+              Errorc("decomposition hint value %d for %c size out "
+                     "of range [0,%d]\n",
+                     (int) xy_decomposition[aRank][1], dimName[aRank], (int) dimUB - xy_decomposition[aRank][0]);
+          }
         if (irregular) /* for irregular grids, all arrays are decomposed identically */
           for (size_t j = 0; j < cdiPioGDsaNum; ++j)
             for (size_t i = 0; i < saMaxRank; ++i)