- 积分
- 211
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-9-11
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
求教:
下面是我用NCL画的2m温度图,运行没有报错,能画出2m温度的填色图,但是后面用 gsn_add_polyline(swks,plot,xpts(i:i+1),ypts(i:i+1),resp) 命令画的方框却没有显示出来,请各位大神帮忙看看有什么问题~~
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
wrf=addfile("/home/sunhw/wrfout-data/wrfout3/wrfout_d01_2011-12-01_00:00:00","r")
tt=wrf->T2
swks=gsn_open_wks("eps","wrfd01-tt")
res = True
pltres = True
mpres = True
res@cnFillOn = True
res@cnMonoFillPattern = True
res@cnMonoFillColor = False
res@cnLineDrawOrder="Predraw"
res@cnFillDrawOrder="Predraw"
contour_tc=wrf_contour(wrf,swks,tt(0,:,:),res)
plot=wrf_map_overlays(wrf,swks,(/contour_tc/),pltres,mpres)
xpts=(/107.0,110.0,110.0,107.0,107.0/)
ypts=(/41.0,41.0,43,43.0,41.0/)
resp=True
resp@gsLineColor="black"
resp@gsLineThicknessF=2.0
dum=new(4,graphic)
do i=0,3
dum(i)=gsn_add_polyline(swks,plot,xpts(i:i+1),ypts(i:i+1),resp)
end do
end
|
|