求助grads取台风中心经纬度

一颗梅子在南京

求助一个问题，如何用grads自动识别台风中心经纬度？
我用NCEP的FNL1×1资料（总共37个时次），想要求出每一个时次的台风中心的经纬度，然后画出台风路径图，
利用论坛里Forcast分享的《GRADS绘制台风(热带气旋)路径-自动识别中心-位置和最低气压保存》方法求经纬度：
'reinit'  ;*第1部分--从海平面气压场（或1000hPa高度场）中找出台风中心位置，输出到文件
'!rm weizhi-w.txt'   ;*使用追加写入，所以要先删除旧数据文件 !是GRADS调用操作系统命令
'open 928s.ctl'   ;*WRF预报结果
i=1
while (i<29) ;*本来应该用<=,但是最后一个位置不对，不用
'set t ' i
t=subwrd(result,4)  ;*取得日期
tc=substr(t,8,5)  ;*取得部分日期串
'd slvl'  ;*海平面气压
'define locx=minloc(min(slvl,lat=10,lat=22),lon=106,lon=121)' ;*搜索最低气压经度
'define locy=minloc(min(slvl,lon=106,lon=121),lat=10,lat=22)' ;*搜索最低气压纬度
'd locx'
locxx=subwrd(result,4)
'd locy'
locyy=subwrd(result,4)
'q gr2w 'locxx' 'locyy ;*GRADS坐标转换成经纬度坐标
x=subwrd(result,3)
y=subwrd(result,6)
rc=write('weizhi-w.txt',x ' 'y ' ' tc,append)  ;*将经纬度位置和时间存入文件
*say rc
i=i+4    ;*只输出08和20时，数据间隔3小时
endwhile
close('weizhi-w.txt')  ;*关闭输出文件
'!cp weizhi-w.txt weizhi-r.txt'  ;*在GS里同一个文件不能既读又写，所以做1个拷贝

遇到一个问题，minloc语句无法求出正确的经纬度，每个时次所求的经纬度一样，且所求的经纬度明显不是台风中心的经纬度。
求大神们帮忙解决！

一颗梅子在南京

我不是故意要点分享资料的，我是想选脚本编辑求助的，可是没法选

wuyu

挺好的，不错