Season: Computes Seasonal Means

In s2dverification: Set of Common Tools for Forecast Verification

Computes Seasonal Means

Description

Computes seasonal means on timeseries organized in a array of any number of dimensions up to 10 dimensions where the time dimension is one of those 10 dimensions.

Usage

Season(var, posdim = 4, monini, moninf, monsup)

Arguments

var	Array containing the timeseries along one of its dimensions.
posdim	Dimension along which to compute seasonal means = Time dimension.
monini	an integer indicating the first month of the time series: 1 to 12.
moninf	an integer indicating the month when to start the seasonal means: 1 to 12.
monsup	an integer indicating the month when to stop the seasonal means: 1 to 12.

Value

Array with the same dimensions as var except along the posdim dimension whose length corresponds to the number of seasons. Partial seasons are not accounted for.

Author(s)

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

Examples

# Load sample data as in Load() example:
example(Load)
leadtimes_dimension <- 4
initial_month <- 11
mean_start_month <- 12
mean_stop_month <- 2
season_means_mod <- Season(sampleData$mod, leadtimes_dimension, initial_month,
                          mean_start_month, mean_stop_month)
season_means_obs <- Season(sampleData$obs, leadtimes_dimension, initial_month,
                          mean_start_month, mean_stop_month)
 
PlotAno(season_means_mod, season_means_obs, startDates, 
       toptitle = paste('winter (DJF) temperatures'), ytitle = c('K'), 
       legends = 'ERSST', biglab = FALSE, fileout = 'tos_season_means.eps')
 

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