data(NCEP_Iberia_tas)
range(NCEP_Iberia_tas$Dates)
## [1] "1990-12-01 GMT" "2000-03-01 GMT"

sub.time <- subsetGrid(NCEP_Iberia_tas, years = 1995:1998)
range(sub.time$Dates)
## [1] "1994-12-01 GMT" "1998-03-01 GMT"

data(NCEP_Iberia_tas)
plotClimatology(climatology(NCEP_Iberia_tas, clim.fun = list(FUN = "mean", na.rm = T)),
                backdrop.theme = "countries", 
                scales = list(draw = T))

sub.geo <- subsetGrid(NCEP_Iberia_tas, lonLim = c(-1, 3), latLim = c(39,43))
plotClimatology(climatology(sub.geo, clim.fun = list(FUN = "mean", na.rm = T)),
                backdrop.theme = "countries", 
                scales = list(draw = T))

data(NCEP_Iberia_tp)
multi <- makeMultiGrid(NCEP_Iberia_tas, NCEP_Iberia_tp)
plotMeanGrid(multi)

sub.geo.multi <- subsetGrid(multi, lonLim = c(-1, 3), latLim = c(39,43))
plotMeanGrid(sub.geo.multi)

sub.Grid <- subsetGrid(multi, var = "tp")
plotClimatology(climatology(sub.Grid, list(FUN = "mean", na.rm = T)), 
                backdrop.theme = "countries", 
                main = attr(sub.Grid$Variable, "longname"))

data(S4_Iberia_tas)
sub.mem <- subsetGrid(S4_Iberia_tas, members = c(1, 3, 5,))
plotClimatology(climatology(sub.mem, list(FUN = "mean", na.rm = T)))