s2dv (version 0.1.1)

Ano: Compute forecast or observation anomalies

Description

This function computes anomalies from a multidimensional data array and a climatology array.

Usage

Ano(data, clim, ncores = NULL)

Arguments

data

A numeric array with named dimensions, representing the model or observational data to be calculated the anomalies. It should involve all the dimensions in parameter 'clim', and it can have more other dimensions.

clim

A numeric array with named dimensions, representing the climatologies to be deducted from parameter 'data'. It can be generated by Clim(). The dimensions should all be involved in parameter 'data' with the same length.

ncores

An integer indicating the number of cores to use for parallel computation. The default value is NULL.

Value

An array with same dimensions as parameter 'data' but with different dimension order. The dimensions in parameter 'clim' are ordered first.

Examples

# NOT RUN {
# Load sample data as in Load() example:
example(Load)
clim <- Clim(sampleData$mod, sampleData$obs)
ano_exp <- Ano(sampleData$mod, clim$clim_exp)
ano_obs <- Ano(sampleData$obs, clim$clim_obs)
ano_exp <- Reorder(ano_exp, c(1, 2, 4, 3))
ano_obs <- Reorder(ano_obs, c(1, 2, 4, 3))
# }
# NOT RUN {
PlotAno(ano_exp, ano_obs, startDates, 
       toptitle = 'Anomaly', ytitle = c('K', 'K', 'K'), 
       legends = 'ERSST', biglab = FALSE, fileout = 'tos_ano.png')
# }