ncep_hgt500_2000 | R Documentation |
warpGrid.R
This exemplary grid contains monthly mean data of geopotential height at 500hPa from NCEP renalysis, for the year 2000.
This R data object was obtained from the UDG server (http://www.meteo.unican.es/udg-tap,
log-in is requiered,
see loginUDG
) by means of function loadGridData
(package loadeR
) in the following manner:
loginUDG("username", "pasword")
ncep_hgt500_2000 <- loadGridData(dataset = "http://meteo.unican.es/tds5/dodsC/ncepReanalysis1/ncepReanalysis1_4xDaily.ncml", var = "hgt@500", latLim = c(20,90), years = 2000, time="DD", aggr.d="mean", aggr.m="mean")
data(ncep_hgt500_2000)
A grid object.
https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.html and http://www.meteo.unican.es/udg-tap.
loadGridData
data("ncep_hgt500_2000")
# Plot in longlat projection
visualizeR::spatialPlot(climatology(ncep_hgt500_2000), backdrop.theme = "coastline")
# Plot in polar stereographic projection
grid <- warpGrid(climatology(ncep_hgt500_2000))
library(visualizeR)
l1 <- get(load(paste0(find.package("visualizeR"), "/countries.rda"))) # world coastline
l1 <- sp::spTransform(l1[[2]], CRSobj = attr(grid$xyCoords, "projection"))
visualizeR::spatialPlot(grid, sp.layout = list(list(l1, first = FALSE)))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.