Rx2day, Rx5day resulting in some infinite values
Two extreme precipitation indices, rx2day and rx5day, produces some infinite values on e.g. GridClim and PTHBV datasets. For both these datasets there are no issues with rx1day.
Example (on bi):
import iris
from climix.metadata import load_metadata
from dask.distributed import Client
client = Client()
fname = "/nobackup/rossby26/users/sm_erhol/extremeEventAttribution/gavle2021/pr_gavle2021_SMHIGridClim_day_19610101-20181230.nc"
# This will be an masked array
cube = iris.load_cube(fname)
index_catalogue = load_metadata()
index = index_catalogue.prepare_indices(["rx2day"])[0]
index_cube = index([cube], client)
assert not np.any(np.isinf(index_cube.data))
When realising the index_cube some overflow warnings are thrown. If the spatial dimensions of the cube are collapsed (average, max), no infinite values are produced.