s2dverification (version 2.5.0)

PlotClim: Plots Climatologies

Description

Plots climatologies as a function of the forecast time for any index output from Clim() and organized in matrix with dimensions: c(nmod/nexp, nmemb/nparam, nltime) or c(nmod/nexp, nltime) for the experiment data c(nobs, nmemb, nltime) or c(nobs, nltime) for the observational data

Usage

PlotClim(exp_clim, obs_clim = NULL, toptitle = "", ytitle = "", monini = 1,  freq = 12, limits = NULL, listexp = c("exp1", "exp2", "exp3"),  listobs = c("obs1", "obs2", "obs3"), biglab = FALSE, leg = TRUE,  fileout = "output_plotclim.eps", sizetit = 1)

Arguments

exp_clim

Matrix containing the experimental data with dimensions: c(nmod/nexp, nmemb/nparam, nltime) or c(nmod/nexp, nltime)

obs_clim

Matrix containing the observational data (optional) with dimensions: c(nobs, nmemb, nltime) or c(nobs, nltime)

toptitle

Main title, optional

ytitle

Title of Y-axis, optional.

monini

Starting month between 1 and 12. Default = 1.

freq

1 = yearly, 12 = monthly, 4 = seasonal, ... Default = 12.

limits

c(lower limit, upper limit): limits of the Y-axis, optional.

listexp

List of experiment names, optional.

listobs

List of observational dataset names, optional.

biglab

TRUE/FALSE for presentation/paper plot. Default = FALSE.

leg

TRUE/FALSE to plot the legend or not.

fileout

Name of the output eps file.

sizetit

Multiplicative factor to scale title size, optional.

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('climatologies'), 
         ytitle = 'K', monini = 11, listexp = c('CMIP5 IC3'), 
         listobs = c('ERSST'), biglab = FALSE, fileout = 'tos_clim.eps')