- 积分
- 76
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-1-8
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
请大家帮忙看看 程序运行没有错误 这个程序是为了在一张图上画出30个台风路径 然而我设置了先画后画的顺序 可是出来的图地图部分把路径盖住了 不知哪里有错 请大神指点
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/csm/contributed.ncl"
wks=gsn_open_wks("png","track")
number=(/"01","02","03","04","05","06","07","08","09","10","11","12","13","14","15","16","17","18","19","20","21","22","23","24","25","26","27","28","29","30"/)
row=(/38,38,14,34,16,27,9,51,47,18,41,17,39,24,6,24,21,30,46,28,28,42,27,44,44,18,9,49,20,14/)
res=True
res@mpDataBaseVersion = "MediumRes"
;res@pmTickMarkDisplayMode = "Always" ; better map tickmarks
res@gsnDraw=False
res@gsnFrame=False ;don't advance frame
res@mpMinLatF=0
res@mpMaxLatF=40
res@mpMinLonF=100
res@mpMaxLonF=160
res@mpFillDrawOrder="PreDraw"
res@mpLandFillColor="grey90"
res@tfPolyDrawOrder="PostDraw"
plot=gsn_csm_map(wks,res)
resn=True
resn@gsLineColor=60
resp=True
do i=0,29
typhoon:=asciiread("/cygdrive/c/Users/Administrator/Desktop/1997 best track/"+number(i)+".txt",-1,"float")
a=row(i)
use:=onedtond(typhoon(8:),(/a,6/))
y:=use(:,2)*0.1
x:=use(:,3)*0.1
c=a-1
do b=0,c
x1:=use(b,3)*0.1
y1:=use(b,2)*0.1
gsn_polyline(wks,plot,x,y,resn)
resp@gsMarkerColor=use(b,1)*10
resp@gsMarkerIndex=1
resp@gsMarkerSizeF=0.01
resp@tfPolyDrawOrder="PostDraw"
gsn_polymarker(wks,plot,x1,y1,resp)
end do
end do
ress=True
ress@txFontHeightF=0.015
ress@txPerimOn=True
ress@tfPolyDrawOrder="PostDraw"
xue=gsn_add_text(wks,plot,"",x(c),y(c),ress)
draw(plot)
frame(wks)
|
-
|