Version 5 (modified by juaco, 7 years ago) (diff) |
---|
Global domain selections
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).
> ex.global <- loadECOMS(dataset = "WFDEI", var = "ps", lonLim = NULL, latLim = NULL, season = c(12,1,2), years = 2001) [2014-06-17 13:18:14] Defining homogeneization parameters for variable "ps" [2014-06-17 13:19:45] Defining geo-location parameters [2014-06-17 13:19:45] Defining time selection parameters [2014-06-17 13:19:58] Done > print(object.size(ex.global), units = "Mb") 178 Mb
Note that, unlike the 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. Therefore, these arguments are always ignored when loading WFDEI or any other non-forecast dataset (they originate a note on screen if not omitted, though).
> ps.mean <- apply(ex.global$Data, FUN = mean, MAR = c(1,2)) > x <- ex.global$xyCoords$x > y <- ex.global$xyCoords$y > # Requires library 'fields'. Please install if not present to reproduce the example > library(fields) > image.plot(x, y, ps.mean, asp = 1, main = "Mean DJF surface pressure 2001" )
Attachments (3)
- ethiopia_drift.png (63.9 KB) - added by juaco 7 years ago.
- ethiopia_mean.png (31.2 KB) - added by juaco 7 years ago.
- image-20140617-134514.png (67.1 KB) - added by juaco 7 years ago.
Download all attachments as: .zip