ncl 非全球contour显示为全球

jstzysq

请教各位大神。我只是在ncl里用自己编的数据，画一个非全球的contour，可是最后显示的图是全球范围的(等值线充满全球)，不是我想要的100~140度，20~40度的范围。我觉得的可能是坐标数组分配时出了什么问题。以前画wrfout数据时，数据是区域性的，导入数据后不用分配坐标数据，就能正确显示范围
begin
lon=fspan(100.,140.,41);东经100到140度
lat=fspan(20.,40.,21);北纬20到40度
lon@units = "degrees_east"
lat@units = "degrees_north"
u=new((/41,21/),"float")
do i=0,40
  do j=0,20
    u(i,j)=sqrt((lon(i)-120)^2+(lat(j)-30)^2);自编数据，以120度、30度为中心的同心圆
  end do
end do
u!0="lon"
u!1="lat"
u&lon=lon
u&lat=lat
wks=gsn_open_wks("x11","aa")
res=True
res@gsnFrame=False
res@gsnDraw=False
res@gsnAddCyclic=False
res@tfDoNDCOverlay=True
res@cnFillOn=True
res@cnLinesOn=False
plot=gsn_csm_contour_map(wks,u,res)
;如果用gsn_csm_contour()，画出的图范围是我所想要的，只是没有地图而已
;另外无论怎么调整(mpMinLonF,...)这四个属性，等值线总是充满整个图
draw(plot)
frame(wks)
end

liyf

res@mpMaxLonF=140
res@mpMaxLatF=40
res@mpMinLonF=100
res@mpMinLatF=20

jstzysq

liyf 发表于 2016-11-23 12:51
res@mpMaxLonF=140
res@mpMaxLatF=40
res@mpMinLonF=100

谢谢。不过会遇到这种情况：比如某数据集只在南北纬80度之间有值，老板仍让我画出全球地图，俩80度之间是等值线，80到90度之间为空白。设置mpMaxLonF......就不管用了

jstzysq

问题解决了。tfDoNDCOverlay要设置为False。还有二维数组要设置为lat x lon，而不是lon x lat

xugk

好东西啊，感谢

圈圈圆圆圈圈

谢楼主！正好我也遇到了这样的问题，参考了您的程序后终于完美的解决了