NCL色斑图及图例

代码如下：

begin
 f = addfile("/cygdrive/d/learn/ncl-test/data/hgt.sfc.nc","r")     ;再分析资料 二进制文件
 var = short2flt( f->hgt(0,:,:))          ;将短型数据转换为浮点型 就是第一个时间，所有经纬度的值
 wks = gsn_open_wks("png", "/cygdrive/d/learn/ncl-test/output/contourmap")
 gsn_define_colormap(wks, "rainbow")    ;GMT_gray灰色色板


 res                = True
 res@gsnLeftString  = "" ;不绘制左表标题
 res@gsnRightString = ""

 res@pmTickMarkDisplayMode = "Always"  ;坐标轴签上增加°符号

 ;;绘制的地图区域
 res@mpMinLatF  = 10 
 res@mpMaxLatF  = 80
 res@mpMinLonF  = 30
 res@mpMaxLonF  = 130

 res@mpFillOn              = False  
 res@mpGridAndLimbOn       = True  ;设置经纬度线
 res@mpGridLineThicknessF  = 0.5
 res@mpGridLineDashPattern = 2
 res@mpGridSpacingF        = 10

 res@cnFillOn  = True
 res@cnLinesOn = False

 res&#