Extraction of specified indices from file with additional threshold file supplied - HELP WANTED

[joelzeder]

Hi @corviday ,
I was wondering whether you could help me with the following issue:
Whenever I try to extract a spcified index (e.g. rx5day) from a file (infile)

infile <- "pr_day_UKESM1-0-LL_ssp585_r1i1p1f2_gn_20150101-20181230.nc"
threshfile <- "UKESM1-0-LL_historical_r1i1p1f2_gn.thresh_1981-2010.nc"

create.indices.from.files(infile, 'output_path', output.filename.template=infile, author.data, climdex.vars.subset = "rx5day", climdex.time.resolution = c("annual"), axis.to.split.on = "Y", fclimdex.compatible = TRUE, base.range = c(1981,2010), parallel = 6, verbose = TRUE, thresholds.files = threshfile, max.vals.millions = 20, cluster.type = "SOCK")

the process fails with the following error message:

Error in check.quantile.validity(quantiles, present.var.list, days.in.base) :
Quantiles must be present for all variables provided.

However, why are quantiles needed in the first place, since Rx5day does not rely on quantiles? Second, all quantiles as computed by create.thresholds.from.file are part of the threshold file threshfile : tx10thresh tx90thresh r95thresh r99thresh. Thus it should actually have all the necessary input right?

However, if I extract the indices the very same way, but from the input file infile <- "pr_day_UKESM1-0-LL_historical_r1i1p1f2_gn_19710101-20101230.nc" the original thresholds were calculated with, then it works with specifying only one index.

Thank you very much!