GRADS绘制登陆台风路径图gs文件.docVIP

*绘制登陆台风路径图 *具体显示日期和气压等内容 year=1989 filename=yingxiang *******************以上为年内台风循环 *绘制地图 *reinit sdfopen f:/air.mon.ltm.nc set grads off set grads off set parea 1 9 1 7 set lev 1000 set cmin 1000 set t 1 set lon 90 160 set lat 5 35 set xlevs 90 100 110 120 130 140 150 160 set ylevs 5 10 15 20 25 30 35 set cint 1000 d air llp=1 while(llp=147) *读取第一行 parap=read(e:/tfnj/filename.txt) parap=sublin(parap,2) lno=subwrd(parap,3) name=subwrd(parap,8) tno=subwrd(parap,4) say lno say name lp=1 while(lp=lno) parap=read(e:/tfnj/filename.txt) parap=sublin(parap,2) *逐行读取 *get day time latitude lontitude* * ymd1=subwrd(parap,1) pres=subwrd(parap,5) *获取时间信息 year1=substr(ymd1,1,4) day1=substr(ymd1,7,2) *获取日期 time1=substr(ymd1,9,2) *获取时次 say day1 say time1 *time1=subwrd(parap,4) lat1=subwrd(parap,3) *分解得到经纬度值 lon1=subwrd(parap,4) latv=lat1/10 lonv=lon1/10 say lonv say latv q w2xy lonv latv *转成直角坐标 lonx=subwrd(result,3) laty=subwrd(result,6) say laty say lonx *draw wxsym 40 lonx laty 0.1 2 5 *在当前点位置画台风符号 if(time1=00) *如果是00时标出日期 set strsiz 0.08 0.08 set string 1 bc *draw string lonx laty-0.1 day1.08h *draw string lonx laty+0.1 pres endif if(lp=2) *如果不是第一个点，当前点与前一点连线 draw line lonxb latyb lonx laty *draw title year No. tno typhoon name Endif lp=lp+1 latyb=laty lonxb=lonx endwhile *******************以下为年内台风循环 llp=llp+1 endwhile printim e:/filenameyear1tnoname.gif x1024 y768 white enable print e:/filename.gmf print disable reinit ***年际循环