s2dverification (version 2.5.0)

Season: Computes Seasonal Means

Description

Computes seasonal means on timeseries organized in a matrix 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

Matrix containing the timeseries along one of its dimensions.

posdim

Dimension along which to compute seasonal means = Time dimension

monini

First month of the time-series: 1 to 12.

moninf

Month when to start the seasonal means: 1 to 12.

monsup

Month when to stop the seasonal means: 1 to 12.

Value

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

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')