Skip to content

Data Access

Cristian Lussana edited this page Jan 30, 2021 · 3 revisions

R-Examples

seNorge_2018 Latest (see also the manual of the R-package ncdf4)

library(ncdf4)
f<-"https://thredds.met.no/thredds/dodsC/senorge/seNorge_2018/Latest/seNorge2018_20210130.nc"
nc<-nc_open(f)
for (i in 1:nc$nvar) print( paste(i,nc$var[[i]]$name))
tg <- nc$var[[6]]
data <- ncvar_get( nc, tg )
print(data)
nc_close(nc)

seNorge_2018 Archive (see also the manual of the R-package ncdf4):

library(ncdf4)
f<-"https://thredds.met.no/thredds/dodsC/senorge/seNorge_2018/Archive/seNorge2018_2021.nc"
nc<-nc_open(f)
for (i in 1:nc$nvar) print( paste(i,nc$var[[i]]$name))
tg <- nc$var[[2]]
varsize <- tg$varsize
ndims <- tg$ndims
nt <- varsize[ndims]
for( i in 1:nt ) {
  start <- rep(1,ndims)
  start[ndims]
  count <- varsize
  count[ndims] <- 1
  data <- ncvar_get( nc, tg, start=start, count=count )
  timeval <- ncvar_get( nc, tg$dim[[ndims]]$name, start=i, count=1 )
  print(paste("Data for variable",tg$name,"at timestep",i,
  " (time value=",timeval,tg$dim[[ndims]]$units,"):"))
  print(data)
}
nc_close(nc)

Clone this wiki locally