关于画散度场出现的问题

清凉月

想画个散度的填色图，ncl一直不报错不出图，删减版程序如下：
（程序里涉及的文件都在附件里）

1. load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/calendar_decode2.ncl"
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
5. 
   
6. begin
   
7. 
   
8. 
   
9. ;-------------计算300hPa位势高度的气候态和距平--------------------------------------------------
   
10. 
    
11. f_data1       = addfile("500hPa_Global_1980_2021_monthly.nc","r")
    
12. lat           = f_data1->latitude
    
13. lon           = f_data1->longitude
    
14. 
    
15. f_data2       = addfile("div.nc","r")     ;改类型
    
16. div         = f_data2->div
    
17. div!0   = "lat"
    
18. div&lat = lat
    
19. div!1   = "lon"
    
20. div&lon = lon
    
21. div@units = "m/s^2"
    
22. div@units = "m/s^2"
    
23. 
    
24. data      = div(:,:)
    
25. 
    
26. wks = gsn_open_wks("png","xxx")
    
27. 
    
28. ;-------------散度--------------------------------------------------
    
29. gsn_define_colormap( wks ,"GreenMagenta16")
    
30. res2              = True
    
31. res2@gsnMaximize  = True
    
32. res2@gsnAddCyclic = False   ;有循环点（全球数据)
    
33. res2@cnLinesOn    = False
    
34. res2@cnFillOn     = True
    
35. ; ;设置地图
    
36. res2@mpFillOn     = False
    
37. res2@mpMinLonF    = -60
    
38. res2@mpMaxLonF    = 180
    
39. res2@mpMinLatF    = 15
    
40. res2@mpMaxLatF    = 90
    
41. res2@pmTickMarkDisplayMode = "Always"
    
42. ;地图数据
    
43. res2@mpDataSetName         = "./database/Earth..4"
    
44. res2@mpDataBaseVersion     = "MediumRes" ; or "Ncarg4_1"
    
45. res2@mpGeophysicalLineColor= "black"
    
46. ;设置地图经纬线
    
47. res2@mpGridAndLimbOn       = True
    
48. res2@mpGridLineThicknessF  = 0.5
    
49. res2@mpGridLineDashPattern = 2
    
50. res2@mpGridSpacingF        = 15
    
51. ;色标垂直放
    
52. res2@lbOrientation         = "Vertical"  
    
53. 
    
54. res2@cnLevelSelectionMode = "ExplicitLevels"
    
55. res2@cnLevels     = (/-7,-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7/)
    
56. res2@cnFillColors = (/2,3,4,5,6,7,8,0,0,11,12,13,14,15,16/)
    
57. 
    
58. data@_FillValue = 32766
    
59. data@missing_value = 32766
    
60. 
    
61. ; printVarSummary(data)
    
62. plot_2 = gsn_csm_contour_map(wks,data,res2)
    
63. 
    
64. end

复制代码

奇怪的是，当把关于mp的设置都删了，把这句话改成这样：plot_2 = gsn_csm_contour(wks,data,res2)之后，出来这么一张图【图片1】
但是用Panoply看了下散度算的也没有问题，别人拿python画了图也没有问题，如【图2】。

所以想问一下这个程序问题出在哪里，感谢~

下载附件居然要扣钱，放网盘了
链接：https://pan.baidu.com/s/1JOfNHKxGSHD9V88P8IT_Wg?pwd=u658
提取码：u658
--来自百度网盘超级会员V6的分享

cookie-o-o

问题可能出在，你的div数据里面包含了nan缺测值，你在读入div后，加一行，replace_ieeenan(div,32767.0,0)试试。另外，div本身有经纬度坐标，后面再给div坐标重新定义无必要。