Fixup 4f9b29c , fixing bug in PointCloud constructor regarding global…

… index computation with nb_parts==1
ecmwf · Sep 22, 2023 · 184134b · 184134b
1 parent d02f399
commit 184134b
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/atlas/functionspace/PointCloud.cc b/src/atlas/functionspace/PointCloud.cc
@@ -158,13 +158,14 @@ PointCloud::PointCloud(const Grid& grid, const grid::Partitioner& _partitioner,
         idx_t j{0};
         gidx_t g{0};
         for (auto p : grid.lonlat()) {
-            if( distribution.partition(g++) == part_ ) {
+            if( distribution.partition(g) == part_ ) {
                 gidx(j) = g+1;
                 ridx(j) = j;
                 lonlat(j, 0) = p.lon();
                 lonlat(j, 1) = p.lat();
                 ++j;
             }
+            ++g;
         }
     }
     else {