grads绘制登陆台风路径图gs文件

《grads绘制登陆台风路径图gs文件》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、*绘制登陆台风路径图*具体显示日期和气压等内容year='1989'filename='yingxiang'*******************以上为年内台风循环*绘制地图*'reinit''sdfopenf:/air.mon.ltm.nc''setgradsoff''setgradsoff''setparea1917''setlev1000''setcmin1000''sett1''setlon90160''setlat535''setxlevs90100110120130140150160''setylevs5101520253035''set

2、cint1000''dair'llp=1while(llp<=147)*读取第一行parap=read('e:/tfnj/'filename'.txt')parap=sublin(parap,2)lno=subwrd(parap,3)name=subwrd(parap,8)tno=subwrd(parap,4)saylnosaynamelp=1while(lp<=lno)parap=read('e:/tfnj/'filename'.txt')parap=sublin(parap,2)*逐行读取*getdaytimelatitudelontitude*

3、*ymd1=subwrd(parap,1)pres=subwrd(parap,5)*获取时间信息year1=substr(ymd1,1,4)day1=substr(ymd1,7,2)*获取日期time1=substr(ymd1,9,2)*获取时次sayday1saytime1*time1=subwrd(parap,4)lat1=subwrd(parap,3)*分解得到经纬度值lon1=subwrd(parap,4)latv=lat1/10lonv=lon1/10saylonvsaylatv'qw2xy'lonv''latv''*转成直角坐标lonx=

4、subwrd(result,3)laty=subwrd(result,6)saylatysaylonx*'drawwxsym40'lonx''laty'0.125'*在当前点位置画台风符号if(time1=00)*如果是00时标出日期'setstrsiz0.080.08''setstring1bc'*'drawstring'lonx''laty-0.1''day1'.08h'*'drawstring'lonx''laty+0.1''presendifif(lp>=2)*如果不是第一个点，当前点与前一点连线'drawline'lonxb''latyb'

5、'lonx''laty*'drawtitle'year'No.'tno'typhoon'name''Endiflp=lp+1latyb=latylonxb=lonxendwhile*******************以下为年内台风循环llp=llp+1endwhile'printime:/'filename''year1''tno''name'.gifx1024y768white''enableprinte:/'filename'.gmf''print''disable''reinit'***年际循环