用grads读取txt中经纬度画站点图

用grads读取txt中经纬度，用于画站点图
　　
　　txt中的站点信息格式如下：
　　
　　stn_num lon lat data1 data2 .... （每个变量用空格格开一行为一个记录）
　　
　　****************draw grid**************************************
　　
　　 'open ******.ctl' *打开一个格点资料用于画底图的值，譬如全球或者区域海温分布，有值就好
　　……
　　
　　……
　　
　　 'd var' *画出底图
　　
　　***********************draw station********************************
　　
　　cnt=1000
　　
　　i=1
　　
　　while(i　　
　　aa=read('westpacific.txt')
　　 aa1=sublin(aa,2)
　　 lat=subwrd(aa1,5)
　　 lon=subwrd(aa1,4)
　　
　　 data＝subwrd(aa1,6)
　　
　　'q w2xy 'lon' 'lat
　　
　　x = subwrd(result,3)
　　y = subwrd(result,6)
　　
　　if(data<=……) *判读控制mark的值，从而体现data的差别
　　'draw mark 8 'x' 'y' 0.01'
　　
　　endif
　　
　　if(data<=……)
　　
　　'draw mark 8 'x' 'y' 0.03'
　　
　　endif
　　
　　if(data<=……)
　　
　　'draw mark 8 'x' 'y' 0.06'
　　
　　endif
　　
　　………
　　i=i+1
　　endwhile
　　ret=close(‘westpacific.txt')
　　'cbarn 1 0 4.9 0.3'
　　'print'
　　'disable print'
　　'reinit'