Skip to content

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.