50年西风急流轴的变率线？

酒国

上图的程序用的是这个网站的。 http://bbs.06climate.com/forum.php?mod=viewthread&tid=17479

上图是显示的50年每一年西风急流轴的位置。这样每个经度上都有50年的急流轴位置的序列，请教一下怎么画出一条50年变率线？

我理解的是每个经度上50年的急流轴的变率，这样做出的折线图横轴是经度值，纵轴是变率值，不知道这样理解对吗？
如果是这样的话，需要输出50年每一年每一个经度上纬度值。请问这在gs怎样实现？用这行命令（ rc=write('d:/Max-U200hPa.txt', lon1，lat1)）试过，结果输出的是虚页尺寸。
gs的主体是：
it=3
while(it<=835)
lon1=17
while(lon1<=80)
    if(lon1=17)
    'set x 'lon1''
    'd maxloc(uwnd(t='it'),lat=0,lat=70)'
    a=sublin(result,2)
    lat1=subwrd(a,4)
    'q gr2xy 'lon1' 'lat1''
    x0=subwrd(result,3)
    y0=subwrd(result,6)          
else
    'set x 'lon1''
    'd maxloc(uwnd(t='it'),lat=0,lat=70)'
    a=sublin(result,2)
    lat1=subwrd(a,4)
    'q gr2xy 'lon1' 'lat1''
    x1=subwrd(result,3)
    y1=subwrd(result,6)
    'set line 1 1 1'
    'draw line 'x0' 'y0' 'x1' 'y1''
    x0=x1
    y0=y1
endif      
lon1=lon1+1
endwhile
it=it+12
endwhile