# CST_AdamontAnalog finds analogous data in the reference dataset to experiment data based on weather types 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. ```r 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. ## Returns 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. ## Examples ```r 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) ``` ## Author(s) Paola Marson, paola.marson@meteo.fr for PROSNOW version Lauriane Batté, lauriane.batte@meteo.fr for CSTools adaptation

CST_AdamontAnalog function