Exmaple6vectorplotsovermaps.docx

Exmaple6——vector plots over maps该示例从netCDF文件中读取数据，并且展示如何在各种地图投影上重复展示矢量图，如何细化矢量数据，如何增加图的大小。这个示例从一些netCDF文件中读取数据，并且在不同的地图投影上创建三种矢量图。为改变矢量和地图设置resources。运行这个示例，必须下载以下文件：gsun06n.ncl，然后键入：ncl gsun06n.ncl。示例6代码及解释 1. load $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl 2. 3. begin /Document/Manuals/Getting_Started/Examples/gsun06n.shtml4. dir = ncargpath(data) /Document/Manuals/Getting_Started/Examples/gsun06n.shtml5. uf = addfile(dir+/cdf/Ustorm.cdf,r) ; Open three netCDF files. /Document/Manuals/Getting_Started/Examples/gsun06n.shtml6. vf = addfile(dir+/cdf/Vstorm.cdf,r) /Document/Manuals/Getting_Started/Examples/gsun06n.shtml7. tf = addfile(dir+/cdf/Tstorm.cdf,r)打开含有风暴数据的三个netCDF文件。uf、vf和tf是在addfile函数中特定的三个引用。 8. /Document/Manuals/Getting_Started/Examples/gsun06n.shtml9. u = uf-u(0,:,:) ; Get u, v, lat, lon data. /Document/Manuals/Getting_Started/Examples/gsun06n.shtml10. v = vf-v(0,:,:) /Document/Manuals/Getting_Started/Examples/gsun06n.shtml11. lat = uf-ulat /Document/Manuals/Getting_Started/Examples/gsun06n.shtml12. lon = uf-ulon在2个netCDF文件引用中读取变量u和V以及坐标变量lat和lon，uf和vf并且将它们存储到本地的NCL变量中。uf-u和vf-v是风速矢量数据的三维数组（同样大小），第一个维度代表时间，第二个维度代表纬度，第三个维度代表经度（用print (/Document/Functions/Built-in/getfilevardims.shtmlgetfilevardims(uf,u))可以验证）。通过符号（0,:,:），选择了的所有经纬度范围内第一个时间步长内的矢量场。 13. 14. wks = gsn_open_wks(x11,gsun06n) ; Open a workstation. 15. /Document/Manuals/Getting_Started/Examples/gsun06n.shtml16. ;----------- Begin first plot -----------------------------------------在地图上绘制矢量场，使用默认的地图投影（圆柱等距）。为这个图设置的唯一的resources是为了确定在地图上的什么位置覆盖上矢量图。 17. 18. resources = True 19. /Document/Manuals/Getting_Started/Examples/gsun06n.shtml20. nlon = dimsizes(lon) /Document/Manuals/Getting_Started/Examples/gsun06n.shtml21. nlat = dimsizes(lat) /Document/Manuals/Getting_Started/Examples/gsun06n.shtml22. resources@/Document/Graphics/Resources/vf.shtmlvfXCStartV = lon(0) ; Define lat/lon corners /Document/Manuals/Getting_Started/Examples/gsun06n.shtml23. resources@/Document/Graphics/