循环报错求解；求大家分享一下t检验思路

WOT

我使用的是ncep的hgt月数据（hgt.mon.mean.nc），想计算一下1988~2001年7-8月的位势高度的平方差，下面的循环是计算平方差的分子，运行时grads报错了
这个循环只是我脚本内容的一部分，其他部分在没有这个循环之前都是没问题的（我敢保证），求求各位老师帮我看一下究竟哪里出错了
还有就是怎么弄t检验，我要做的t检验和下图这个一样，孩子不会弄，都折腾了好久了，要好心人可以分享一下t检验的思路吗
i=487
ssum=0
while i<=643
'set t ' i
'define hg7=hgt'
'set t ' i+1
'define hg8=hgt'
'define ssum=pow(hg7+hg8,2)+ssum'
i=i+12
endwhile

孙玻玻

请问您解决了吗？