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