Clim: Computes Bias Corrected Climatologies

In s2dverification: Set of Common Tools for Forecast Verification

Computes Bias Corrected Climatologies

Description

This function computes only per-pair climatologies from the experimental and observational matrices output from Load(). To compute plain climatologies from only experimental or observational data from Load(), the following code can be used:
clim <- array(apply(obs_data, c(1, 4, 5, 6), mean),
dim = dim(obs_datta)[-c(2, 3)])
The function Clim() computes per-pair climatologies 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() 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	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.

Author(s)

History:
0.9 - 2011-03 (V. Guemas) - Original code
1.0 - 2013-09 (N. Manubens) - Formatting to R CRAN

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')
 

s2dverification documentation built on April 20, 2022, 9:06 a.m.