Computes Bias Corrected Climatologies
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:
per-pair method (Garcia-Serrano and Doblas-Reyes, CD, 2012)
Kharin method (Karin et al, GRL, 2012)
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.
Clim(var_exp, var_obs, memb = TRUE, kharin = FALSE, NDV = FALSE)
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. |
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. |
History:
0.9 - 2011-03 (V. Guemas, virginie.guemas@ic3.cat) - Original code
1.0 - 2013-09 (N. Manubens, nicolau.manubens@ic3.cat) - Formatting to R CRAN
# 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')Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.