Computes Forecast or Observed Anomalies
This function computes anomalies from any experimental or observational
matrix output from Load() and their climatologies output from
Clim().
Ano(var, clim)
var |
Model or observational data: |
clim |
Climatologies from clim: c(nmod/nexp/nobs, nltime) |
Array with same dimensions as 'var'.
History:
0.1 - 2012-03 (V. Guemas, virginie.guemas@ic3.cat) - Original code
1.0 - 2013-09 (N. Manubens, nicolau.manubens@ic3.cat) - Formatting to R CRAN
# 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')Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.