Plot Raw Or Smoothed Anomalies

Share:

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

1
2
3
4
5
6
7
8
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.

Author(s)

History:
0.1 - 2011-03 (V. Guemas, virginie.guemas at ic3.cat) - Original code
1.0 - 2013-09 (N. Manubens, nicolau.manubens at ic3.cat) - Formatting to CRAN

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
# 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')

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.