s2dverification (version 2.5.0)

Clim: Computes Per-pair/Kharin/Fuckar Climatologies

Description

This function computes climatologies from the experimental and observational matrices output from Load() using one of the following methods:

• 1) per-pair method (Garcia-Serrano and Doblas-Reyes, CD, 2012)
• 2) Kharin method (Karin et al, GRL, 2012)
• 3) Fuckar method (Fuckar et al, GRL, 2014)

Clim function computes climatologies using the startdates covered by the whole experiments/observational data sets. The startdates not available for all the data (model and obs) are excluded when computing the climatologies.

Usage

Clim(var_exp, var_obs, memb = TRUE, kharin = FALSE, NDV = FALSE)

Arguments

var_exp

Model data: c(nmod/nexp, nmemb/nparam, nsdates, nltime) up to c(nmod/nexp, nmemb/nparam, nsdates, nltime, nlevel, nlat, nlon)

var_obs

Observational data: c(nobs, nmemb, nsdates, nltime) up to c(nobs, nmemb, nsdates, nltime, nlevel, nlat, nlon)

memb

memb: TRUE/FALSE (1 climatology for each member). Default = TRUE.

kharin

TRUE/FALSE (if Kharin method is applied or not). Default = FALSE.

NDV

TRUE/FALSE (if Fuckar method is applied or not). Default = FALSE.

Value

clim_exp

Array with same dimensions as var_exp except the third (starting dates) and, depending on the parameters, the second (members), which disappear.

clim_obs

Array with same dimensions as var_obs except the third (starting dates) and, depending on the parameters, the second (members), which disappear.

Examples

# Load sample data as in Load() example:
example(Load)
clim <- Clim(sampleData$mod, sampleData$obs)
PlotClim(clim$clim_exp, clim$clim_obs, 
         toptitle = paste('sea surface temperature climatologies'), 
         ytitle = 'K', monini = 11, listexp = c('CMIP5 IC3'), 
         listobs = c('ERSST'), biglab = FALSE, fileout = 'tos_clim.eps')