CSIndicators (version 1.2.0)

MultiVarExceedingThreshold: Binary probabilities based on climate thresholds (internal function)

Description

This function computes a binary array based on the thresholds imposed to one or multiple climate variables by the user (e.g. thresholds for temperature, precipitation and relative humidity). Values of 1 imply that the conditions are met, while values of 0 imply that the conditions are not met.

Usage

MultiVarExceedingThreshold(
  data,
  thresholds,
  operators,
  dates = NULL,
  start = NULL,
  end = NULL,
  time_dim = "time",
  fun = max,
  ncores = NULL,
  na.rm = FALSE,
  ...
)

Value

An array containing the binary outputs for whether all arrays in data pass the required thresholds.

Arguments

data

A list of arrays containing the climate variable(s)

thresholds

A list of list of scalar/vectors/arrays for each s2dv_cube object in data. If the threshold is an array and some dimensions are the same as in data, the dimension must be equal.

operators

A list of scalar/vectors containing the operators to be applied to each climate variable in data.

dates

A multidimensional array of dates with named dimensions matching the temporal dimensions on parameter 'data'. By default it is NULL, to select aperiod this parameter must be provided.

start

An optional parameter to define the initial date of the period to select from the data by providing a list of two elements: the initial date of the period and the initial month of the period. By default it is set to NULL and the indicator is computed using all the data provided in data.

end

An optional parameter to defined the final date of the period to select from the data by providing a list of two elements: the final day of the period and the final month of the period. By default it is set to NULL and the indicator is computed using all the data provided in data.

time_dim

A character string indicating the name of the dimension to compute the indicator. By default, it is set to 'time'.

fun

A function to be applied along the time dimension. By default, the max is computed over time_dim. To keep the original temporal resolution (i.e., return values for each time step), set `fun = NULL`, in which case no aggregation is applied.

ncores

An integer indicating the number of cores to use in parallel computation.

na.rm

A logical value indicating whether to ignore NA values (TRUE) or not (FALSE). Default value is FALSE.

...

Additional arguments for parameter 'fun'.

Author

Javier Corvillo & Alberto Bojaly,

Examples

# Simple arrays with single thresholds
tas <- array(rnorm(100, mean = 25, sd = 5), dim = c(time = 10, lat = 5, lon = 2))
prlr <- array(abs(rnorm(100, mean = 50, sd = 20)), dim = c(time = 10, lat = 5, lon = 2))

# Check where temperature > 25 AND precipitation > 40
result <- MultiVarExceedingThreshold(
  data = list(tas, prlr),
  operators = list(">", ">"),
  thresholds = list(25, 40)
)