3 | | In this example we will load data for Europe for the variable surface (2m) minimum temperature (`var = tasmin`), for the first two members (`members = 1:2`) of the CFSv2 hindcast (dataset = `CFSv2_seasonal_16`), considering the wintertime (DJF, `season = c(12,1,2)`) for the 10-year period 2001-2010 (`years = 2001:2010`), according to the forecast of previous September (`leadMonth = 3`). The original variable is stored as 6-hourly data for this particular dataset. Instead of loading the whole 6-hourly time series, or filtering by a particular time as in the [http://meteo.unican.es/ecoms-udg/RPackage/Examples/pointSelection previous example] we will retrieve the daily mean values, by setting the argument `time = "DD"`, that internally computes the daily mean from the 6-hourly instantaneous values. |
| 3 | In this example we will load data for Europe for the variable surface (2m) minimum temperature (`var = tasmin`), for the first two members (`members = 1:2`) of the CFSv2 hindcast (dataset = `CFSv2_seasonal_16`), considering the wintertime (DJF, `season = c(12,1,2)`) for the 10-year period 2001-2010 (`years = 2001:2010`), according to the forecast of previous September (`leadMonth = 3`). The original variable is stored as 6-hourly data for this particular dataset. We will retrieve the daily mean values, by setting the argument `time = "DD"`, that internally computes the daily minimum from the 6-hourly instantaneous values (see the NOTE below when executing the command): |
8 | | [2014-06-17 12:47:49] Defining homogeneization parameters for variable "tasmin" |
9 | | NOTE: daily mean will be calculated from the 6-h instantaneous model output |
10 | | [2014-06-17 12:47:49] Defining geo-location parameters |
11 | | [2014-06-17 12:47:49] Defining initialization time parameters |
12 | | [2014-06-17 12:47:54] Retrieving data subset ... |
13 | | [2014-06-17 12:54:33] Done |
| 8 | [2014-09-02 16:45:58] Defining homogeneization parameters for variable "tasmin" |
| 9 | NOTE: daily minimum will be calculated from the 6-h model output |
| 10 | [2014-09-02 16:45:58] Defining geo-location parameters |
| 11 | [2014-09-02 16:45:58] Defining initialization time parameters |
| 12 | [2014-09-02 16:46:03] Retrieving data subset ... |
| 13 | [2014-09-02 16:52:57] Done |
31 | | > library(fields) # Install if not available to reproduce the example |
32 | | > member1 <- apply(ex2$Data[,,,1], FUN = mean, MARGIN = c(2,3)) |
33 | | > member2 <- apply(ex2$Data[,,,2], FUN = mean, MARGIN = c(2,3)) |
34 | | > x <- ex2$xyCoords$x |
35 | | > y <- ex2$xyCoords$y |
36 | | > par(mfrow = c(1,2)) |
37 | | > image.plot(x,y,member1, asp = 1, main = "Member 1") |
38 | | > world(add = TRUE) |
39 | | > image.plot(x,y,member2, asp = 1, main = "Member 2") |
40 | | > world(add = TRUE) |
| 31 | plotMeanField(ex2, multi.member = TRUE) |