Guest User

Untitled

a guest
Oct 24th, 2017
93
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.78 KB | None | 0 0
  1. 7 variables (excluding dimension variables):
  2. double lat[rlon,rlat]
  3. bounds: lat_vertices
  4. long_name: latitude coordinate
  5. standard_name: latitude
  6. units: degrees_north
  7. double lat_vertices[vertices,rlon,rlat]
  8. units: degrees_north
  9. double lon[rlon,rlat]
  10. bounds: lon_vertices
  11. long_name: longitude coordinate
  12. standard_name: longitude
  13. units: degrees_east
  14. double lon_vertices[vertices,rlon,rlat]
  15. units: degrees_east
  16. double time_bnds[bnds,time]
  17. float hfss[rlon,rlat,time]
  18. standard_name: surface_upward_sensible_heat_flux
  19. long_name: Surface Upward Sensible Heat Flux
  20. units: W m-2
  21. original_name: HFX
  22. cell_methods: time: mean
  23. missing_value: 1.00000002004088e+20
  24. _FillValue: 1.00000002004088e+20
  25.  
  26. grid_mapping: rotated_latitude_longitude
  27. coordinates: lat lon
  28. positive: up
  29. int rotated_pole[]
  30. grid_mapping_name: rotated_latitude_longitude
  31. grid_north_pole_latitude: 39.25
  32. grid_north_pole_longitude: -162
  33.  
  34. library(raster)
  35. library(rasterVis)
  36. library(rworldmap)
  37.  
  38. data(wrld_simpl)
  39. colnames(mat) <-c("lon","lat","value")
  40.  
  41. x.min=min(mat$lon)
  42. x.max=max(mat$lon)
  43. y.min=min(mat$lat)
  44. y.max=max(mat$lat)
  45.  
  46. v.min=min(mat$value)
  47. v.max=max(mat$value)
  48.  
  49. breaks=seq(v.min,v.max,2)
  50. colour=rev(rainbow(length(breaks), alpha=1))
  51.  
  52. ras <- raster(nrows=453, ncols=453,xmn=x.min,xmx=x.max, ymn=y.min, ymx=y.max)
  53. rvm <- matrix(mat[,3],453,453,byrow=TRUE)
  54. ras[]=rvm[nrow(rvm):1,]
  55.  
  56. levelplot(ras,at=breaks,col.regions=colour, margin = FALSE) +
  57. layer(sp.polygons(wrld_simpl, fill = NA, col = "grey50"))
Add Comment
Please, Sign In to add comment