| 1 | = Global domain selections |
| 2 | |
| 3 | For this example we will load the gridded observational dataset WFDEI (`dataset = "WFDEI"`). In particular, we will load the near-surface atmospheric pressure (`var = ps`) for all land areas globally in boreal winter (DJF, `season = c(12,1,2)`) for the year 2010 (`years = 2010`). Note that, for full spatial domain requests, arguments `lonLim` and `latLim` are set to `NULL` (the default, which is equivalent to just omitting both arguments). |
| 4 | |
| 5 | {{{ |
| 6 | > ex.global <- loadECOMS(dataset = "WFDEI", var = "ps", lonLim = NULL, latLim = NULL, season = c(12,1,2), years = 2001) |
| 7 | [2014-06-17 13:18:14] Defining homogeneization parameters for variable "ps" |
| 8 | [2014-06-17 13:19:45] Defining geo-location parameters |
| 9 | [2014-06-17 13:19:45] Defining time selection parameters |
| 10 | [2014-06-17 13:19:58] Done |
| 11 | > print(object.size(ex.global), units = "Mb") |
| 12 | 178 Mb |
| 13 | }}} |
| 14 | |
| 15 | Note that, unlike the [http://meteo.unican.es/ecoms-udg/RPackage/Examples previous examples], we have omitted the `leadMonth` and `members` arguments, as this is not a forecast dataset, but a gridded observational dataset lacking the initialization and ensemble dimensions. If any of this arguments is specified when loading `WFDEI` or any other non-forecast dataset, they will be simply ignored (originating a note on screen). |
| 16 | |
| 17 | {{{ |
| 18 | > ps.mean <- apply(ex.global$Data, FUN = mean, MAR = c(1,2)) |
| 19 | > x <- ex.global$xyCoords$x |
| 20 | > y <- ex.global$xyCoords$y |
| 21 | > # Requires library 'fields'. Please install if not present to reproduce the example |
| 22 | > library(fields) |
| 23 | > image.plot(x, y, ps.mean, asp = 1, main = "Mean DJF surface pressure 2001" ) |
| 24 | }}} |
| 25 | |
| 26 | [[Image(image-20140617-134514.png)]] |
| 27 | |
| 28 | {{{#!comment |