Provides simplified access to OPeNDAP Real-time Ocean Forecast System (rtofs) data served by the NOMADS OPeNDAP server.
Using devtools
devtools::install_github("BigelowLab/rtofs")
Here's an example the accesses yesterday's surface current forecasts and draws then using rasterVis.
library(rtofs)
library(raster)
library(rasterVis)
uri = daily_glo2ds_uri()
NC = rtofs_daily_glo2ds(uri)
U = rtofs_glo2ds_get_raster(NC, bb = c(-72, -63, 39, 46), param = 'u_velocity')
V = rtofs_glo2ds_get_raster(NC, bb = c(-72, -63, 39, 46), param = 'v_velocity')
day3 = raster::stack(U[[3]], V[[3]])
names(day3) <- c("U", "V")
day5 = raster::stack(U[[5]], V[[5]])
names(day5) <- c("U", "V")
p3 <- vectorplot(day3, par.settings=PuOrTheme(), isField = 'dXY', margin=FALSE, main = 'Day 3 currents')
p5 <- vectorplot(day5, par.settings=PuOrTheme(), isField = 'dXY', margin=FALSE, main = 'Day 5 currents')
print(p3, split=c(1, 1, 1, 2), more=TRUE)
print(p5, split=c(1, 2, 1, 2))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.