NCL出图呈现方块状，带角，很丑

Herminia

画一个SST梯度随时间的变化，同时做了该区域纬向平均，平均时将小值区域设为缺测，画出的图呈现方块状，脚本如下（只贴了比较重要的片段）：

begin
;*****读取数据。。。。
wtemp=fin[:]->water_temp(:,0,{LATA:LATB},{LONA:LONB})  ;;;1维time，2维depth
;;;;;;计算月平均
do md=0,11
wtav(md,:,:)=dim_avg_n(wtemp(md*6:md*6+5,:,:),0)
end do
;;;;通过计算x.y方向的中央差，计算梯度
dwdy=center_finite_diff_n(wtav(:,:,:),dy,False,0,1)  ;;;;dy为y方向的格距大小
do m=0,30
dwdx(:,m,:)=center_finite_diff_n(wtav(:,m,:),dx(m),False,0,1)  ;;;dx为x方向的格距大小
end do
;;;;求梯度
grade=10000*sqrt(dwdy^2+dwdx^2)   
;;;去除小值
dtdx=10000*0.05/1852   
grade=where((grade.gt.dtdx),grade,grade@_FillValue)
;;;;;;;;;纬向平均;;;;;;;;;;;;;;;;;;
latav=dim_avg_n_Wrap(grade,1)
;;;;;画图
res@gsnAddCyclic=False
res@cnLevelSelectionMode="ManualLevels"
res@cnMinLevelValF       = 0.35                 ; min level
res@cnMaxLevelValF       = 0.65                 ; max level;
res@cnLevelSpacingF      = 0.05                 
end