s2dverification (version 2.10.3)

Season: 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.

Examples

# NOT RUN {
# 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)
 
# }
# NOT RUN {
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')
 
# }