s2dverification (version 2.5.0)

PlotAno: Plot Raw Or Smoothed Anomalies

Description

Plots timeseries of raw or smoothed anomalies of any index output from Load() or Ano() or or Ano_CrossValid() or Smoothing() and organized in matrices with dimensions: c(nmod/nexp, nmemb/nparam, nsdates, nltime) for the model data c(nobs, nmemb, nsdates, nltime) for the observational data

Usage

PlotAno(exp_ano, obs_ano = NULL, sdates, toptitle = c("", "", "", "", "", "",  "", "", "", "", "", "", "", "", ""), ytitle = c("", "", "", "", "",  "", "", "", "", "", "", "", "", "", ""), limits = NULL,  legends = NULL, freq = 12, biglab = FALSE, fill = TRUE, memb = TRUE,  ensmean = TRUE, linezero = FALSE, points = FALSE, vlines = NULL,  fileout = c("output1_plotano.eps", "output2_plotano.eps",  "output3_plotano.eps", "output4_plotano.eps", "output5_plotano.eps"),  sizetit = 1)

Arguments

exp_ano

Array containing the experimental data: c(nmod/nexp, nmemb/nparam, nsdates, nltime).

obs_ano

Optional matrix containing the observational data: c(nobs, nmemb, nsdates, nltime)

sdates

List of starting dates: c('YYYYMMDD','YYYYMMDD').

toptitle

Main title for each experiment: c('',''), optional.

ytitle

Title of Y-axis for each experiment: c('',''), optional.

limits

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

legends

List of observational dataset names, optional.

freq

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

biglab

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

fill

TRUE/FALSE if the spread between members should be filled. Default = TRUE.

memb

TRUE/FALSE if all members/only the ensemble-mean should be plotted. Default = TRUE.

ensmean

TRUE/FALSE if the ensemble-mean should be plotted. Default = TRUE.

linezero

TRUE/FALSE if a line at y=0 should be added. Default = FALSE.

points

TRUE/FALSE if points instead of lines should be shown. Default = FALSE.

vlines

List of x location where to add vertical black lines, optional.

fileout

Name of the output eps file for each experiment: c('','').

sizetit

Multiplicative factor to scale title size, optional.

Examples

# 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)
runmean_nb_months <- 12
dim_to_smooth <- 4  # Smooth along lead-times
smooth_ano_exp <- Smoothing(ano_exp, runmean_nb_months, dim_to_smooth)
smooth_ano_obs <- Smoothing(ano_obs, runmean_nb_months, dim_to_smooth)
PlotAno(smooth_ano_exp, smooth_ano_obs, startDates, 
        toptitle = paste('smoothed anomalies'), ytitle = c('K', 'K', 'K'), 
        legends = 'ERSST', biglab = FALSE, fileout = 'tos_ano.eps')