Make smaller chunks to improve parallelisation at cost of scheduling

midas/doy_statistics.py

@@ -385,7 +385,7 @@ class DOYStatistics:
             data.append(batch[..., self.inds[idx_d]])
         data = da.stack(data, axis=-3)
 
-        chunks = (1000,) * (data.ndim - 3) + ("auto", -1, -1)
+        chunks = (1000,) * (data.ndim - 3) + (10, -1, -1)
         data = data.rechunk(chunks)
 
         logger.info("Starting sorting")