Currently our code will hang if the spatial field is not evenly divisible by the spatial decomposition. We should raise an exception during the call to get_local_slices to alert the user early to this type of unsupported configuration.
Once this issue is resolved, we can consider opening an issue to support this type of configuration.