Exmaple5——contour plots over maps.docx

Exmaple5——contour plots over maps这个示例从netCDF文件中读取数据，并且展示如何命名维度。如何获得现有的颜色图并且改变它的值，如何在各种地图投影上叠加等值线图，如何填充特定的等值线，如何掩蔽，如何在你任意想要的位置绘制文本字符串。这个示例从netCDF文件中读取数据并且在不同的地图投影上创建了4个等值线图。从等值线和地图两个方面改变了resources。运行这个示例，必须下载以下文件：gsun05n.ncl，然后键入：ncl gsun05n.ncl示例5代码及解释 1. load $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl 2. 3. begin4. data_dir = ncargpath(data) ; Get ready to open three netCDF files.5. cdf_file1 = addfile(data_dir + /cdf/941110_P.cdf,r)6. cdf_file2 = addfile(data_dir + /cdf/sstdata_netcdf.nc,r)7. cdf_file3 = addfile(data_dir + /cdf/Pstorm.cdf,r)打开三个netCDF文件，cdf_file1、cdf_file2和cdf_file3是在addfile函数中指定的三个文件。 8.9. psl = cdf_file1-Psl ; Store some data from the three netCDF10. sst = cdf_file2-sst ; files to local variables.11. pf = cdf_file3-p从刚刚打开的三个netCDF中读取一些数据，并且将它们存储在本地的NCL变量中。psl是一个在纬度/经度格点测量的压力值的2维数组，sst是一个在纬度/经度和12个不同时间步长上测量的海洋表面温度的3维数组。pf是在纬度/经度格点和64个不同时间步长上测量的压力值的3维数组。请注意：本地NCL变量和netCDF文件变量中所使用的不同的变量名字（如，pf和cdf_file-p）。这是可以选择的，此示例中是为了展示你可以给你想要命名的变量随意命名。 12.13. psl@nlon = dimsizes(psllon) ; Store the sizes of the lat/lon14. psl@nlat = dimsizes(psllat) ; arrays as attributes of the psl15. pf@nlon = dimsizes(pflon) ; and pf variables.16. pf@nlat = dimsizes(pflat)Psl和pf都包含叫做“lat”和“lon”的坐标变量，因此获得这些属性变量的常量并且把它们的值以psf和pf的属性存储起来。 17.18. sst!1 = lat ; Name dimensions 0 and 1 of sst19. sst!2 = lon ; lat and lon.20. sstlat = cdf_file2-lat ; Create coordinate variables21. sstlon = cdf_file2-lon ; for sst and store the sizes22. sst@nlon = dimsizes(sstlon) ; of the arrays as attributes 23. sst@nlat = dimsizes(sstlat) ; of sst.sst唯一的坐标变量是“time”，和它的第一个维度相对应（这是在先前的文件中事先确定的）。使用iscoord函数来确定变量是否含有该特定的坐标变量。在18-21行中，因为sst只有第一个维度是命名过的，你可以给sst创建named dimensions和叫做“lat”和“lon”的第二个、第三个坐标变量。在22-23行，存储lat和lon数组的长度作为sst变量。请记住：为了创建另一个变量的坐标变量，它必须与变量的一个named dimension有相同的名字。 24. 25. wks = gsn_open_wks(x11,gsun05n) ; Open a workstation. 26.27. ;----------- Begin first plot