rx5day index fails with NGCD input
Running the following command line:
climix -v -e -x rx5day /home/rossby/joint_exp/climix/1/testdata/NGCD/input_data/SverigeScaled/NGCD_RR_type2_*.nc
results in the following traceback:
Traceback (most recent call last):
File "/home/sm_joalo/.conda/envs/climix-devel/bin/climix", line 8, in <module>
sys.exit(main())
File "/home/sm_joalo/dev/repos/climix/climix/main.py", line 293, in main
do_main(
File "/home/sm_joalo/dev/repos/climix/climix/main.py", line 271, in do_main
result = index(input_data, client=scheduler.client, sliced_mode=sliced_mode)
File "/home/sm_joalo/dev/repos/climix/climix/index.py", line 52, in __call__
aggregated = multicube_aggregated_by(
File "/home/sm_joalo/dev/repos/climix/climix/iris.py", line 126, in multicube_aggregated_by
result = list(map(agg, groupby_subcubes))
File "/home/sm_joalo/dev/repos/climix/climix/iris.py", line 109, in agg
result = aggregate(data, axis=dimension_to_groupby, cube=cube, **kwargs)
File "/home/sm_joalo/.conda/envs/climix-devel/lib/python3.10/site-packages/iris/analysis/__init__.py", line 549, in lazy_aggregate
return self.lazy_func(data, axis=axis, **kwargs)
File "/home/sm_joalo/dev/repos/climix/climix/index_functions/index_functions.py", line 449, in lazy_func
return super().call_func(thresholded_data, axis, **kwargs)
File "/home/sm_joalo/dev/repos/climix/climix/index_functions/index_functions.py", line 354, in call_func
rolling_view = np.lib.stride_tricks.sliding_window_view(
File "<__array_function__ internals>", line 180, in sliding_window_view
File "/home/sm_joalo/.conda/envs/climix-devel/lib/python3.10/site-packages/dask/array/core.py", line 1760, in __array_function__
return da_func(*args, **kwargs)
File "/home/sm_joalo/.conda/envs/climix-devel/lib/python3.10/site-packages/dask/array/overlap.py", line 815, in sliding_window_view
safe_chunks = tuple(
File "/home/sm_joalo/.conda/envs/climix-devel/lib/python3.10/site-packages/dask/array/overlap.py", line 816, in <genexpr>
ensure_minimum_chunksize(d + 1, c) for d, c in zip(depths, x.chunks)
File "/home/sm_joalo/.conda/envs/climix-devel/lib/python3.10/site-packages/dask/array/overlap.py", line 353, in ensure_minimum_chunksize
raise ValueError(
ValueError: The overlapping depth 5 is larger than your array 1.