CST_AdamontAnalog {CSTools}

CST_AdamontAnalog finds analogous data in the reference dataset to experiment data based on weather types

Description

This function searches for analogs in a reference dataset for experiment data, based on corresponding weather types. The experiment data is typically a hindcast, observations are typically provided by reanalysis data.

This function searches for analogs in a reference dataset for experiment data, based on corresponding weather types. The experiment data is typically a hindcast, observations are typically provided by reanalysis data.

Usage

CST_AdamontAnalog(
  exp,
  obs,
  wt_exp,
  wt_obs,
  nanalogs,
  method = "pattcorr",
  thres = NULL,
  search_obsdims = c("member", "sdate", "ftime"),
  londim = "lon",
  latdim = "lat"
)

AdamontAnalog(
  exp,
  obs,
  wt_exp,
  wt_obs,
  nanalogs = 5,
  method = "pattcorr",
  thres = NULL,
  search_obsdims = c("member", "sdate", "ftime"),
  londim = "lon",
  latdim = "lat"
)

Arguments

exp

A multidimensional array with named dimensions containing the experiment data.

obs

A multidimensional array with named dimensions containing the reference data. Note that lat/lon dimensions need to be the same as exp.

wt_exp

Corresponding weather types (same dimensions as exp$data but lat/lon).

wt_obs

Corresponding weather types (same dimensions as obs$data but lat/lon).

nanalogs

Integer defining the number of analog values to return (default: 5).

method

A character string indicating the method used for analog definition. It can be:

• 'pattcorr': pattern correlation.

• 'rain1' (for precip patterns): rain occurrence consistency.

• 'rain01' (for precip patterns): rain occurrence/non occurrence consistency

thres

Real number indicating the threshold to define rain occurrence/non occurrence in rain (0)1.

search_obsdims

List of dimensions in obs along which analogs are searched for.

londim

Name of longitude dimension.

latdim

Name of latitude dimension.

Value

analog_vals An object of class s2dv_cube containing nanalogs analog values for each value of exp input data.

analog_vals An array containing nanalogs analog values.

Author(s)

Paola Marson, [e-mail] for PROSNOW version

Lauriane Batté, [e-mail] for CSTools adaptation

Examples

wt_exp <- sample(1:3, 15*6*3, replace = TRUE)
dim(wt_exp) <- c(dataset = 1, member = 15, sdate = 6, ftime = 3)
wt_obs <- sample(1:3, 6*3, replace = TRUE)
dim(wt_obs) <- c(dataset = 1, member = 1, sdate = 6, ftime = 3)
exp <- NULL
exp$data <- 1 : c(1 * 15 * 6 * 3 * 8 * 8)
dim(exp$data) <- c(dataset = 1, member = 15, sdate = 6, ftime = 3,
                  lat = 8, lon = 8)
class(exp) <- 's2dv_cube'
obs <- NULL
obs$data <- 101 : c(100 + 1 * 1 * 6 * 3 * 8 * 8)
dim(obs$data) <- c(dataset = 1, member = 1, sdate = 6, ftime = 3,
                  lat = 8, lon = 8)
class(obs) <- 's2dv_cube'
analog_vals <- CST_AdamontAnalog(exp = exp, obs = obs, wt_exp = wt_exp, 
                                wt_obs = wt_obs, nanalogs = 2)
wt_exp <- sample(1:3, 15*6*3, replace = TRUE)
dim(wt_exp) <- c(dataset = 1, member = 15, sdate = 6, ftime = 3)
wt_obs <- sample(1:3, 6*3, replace = TRUE)
dim(wt_obs) <- c(dataset = 1, member = 1, sdate = 6, ftime = 3)
exp <- 1 : c(1 * 15 * 6 * 3 * 8 * 8)
dim(exp) <- c(dataset = 1, member = 15, sdate = 6, ftime = 3, lat = 8, lon = 8)
obs <- 101 : c(100 + 1 * 1 * 6 * 3 * 8 * 8)
dim(obs) <- c(dataset = 1, member = 1, sdate = 6, ftime = 3, lat = 8, lon = 8)
analog_vals <- AdamontAnalog(exp = exp, obs = obs, wt_exp = wt_exp, 
                            wt_obs = wt_obs, nanalogs = 2)