Advertisement
Guest User

Untitled

a guest
Oct 26th, 2016
77
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.83 KB | None | 0 0
  1. require(ncdf4)
  2.  
  3. # to get all the future filenames:
  4. filenames <- list.files(path="G:/future",pattern="_.....nc", full.names=TRUE,recursive=T)
  5. # just change to path for the "present" ones
  6.  
  7. # let's run this for just one file:
  8. fn <- filenames[1]
  9.  
  10. # load one of the roms files
  11. roms <- nc_open(fn)
  12.  
  13. # to see what's in the roms file:
  14. print(roms)
  15.  
  16. # load temperature
  17. temp <- ncvar_get(roms,"temp")
  18. dim(temp)
  19.  
  20. # load lat & long
  21. lon <- ncvar_get(roms,"lon_rho")
  22. lat <- ncvar_get(roms,"lat_rho")
  23.  
  24. # load depth
  25. depth <- ncvar_get(roms,'h')
  26.  
  27. # load mask (land or sea)
  28. mask <- ncvar_get(roms,"mask_rho")
  29.  
  30. # load time
  31. ocean_time <- ncvar_get(roms,varid = 'ocean_time')
  32. base_time <- "2060-01-01" # use this for "future" ROMS files
  33. # base_time <- "1990-01-01" # use this for "present" ROMS files
  34. ocean_time <- as.POSIXct(unlist(ocean_time), origin = base_time)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement