关于NC文件循环求平均的问题

蓉蓉不着急

求助，想把下面这段程序写成循环
'reinit'
'set gxout shaded'
'sdfopen d:\ozone\go3.2003.nc'
'set z 3'
'define a=ave(go3,t=60,t=151)'
'close 1'
'sdfopen d:\ozone\go3.2004.nc'
'set z 3'
'define b=ave(go3,t=61,t=152)'
'close 1'
'sdfopen d:\ozone\go3.2005.nc'
'set z 3'
'define c=ave(go3,t=60,t=151)'
'close 1'
'sdfopen d:\ozone\go3.2006.nc'
'set z 3'
'define d=ave(go3,t=60,t=151)'
'close 1'
'sdfopen d:\ozone\go3.2007.nc'
'set z 3'
'define e=ave(go3,t=60,t=151)'
'close 1'
'sdfopen d:\ozone\go3.2008.nc'
'set z 3'
'define f=ave(go3,t=61,t=152)'
'close 1'
'sdfopen d:\ozone\go3.2009.nc'
'set z 3'
'define g=ave(go3,t=60,t=151)'
'close 1'
'sdfopen d:\ozone\go3.2010.nc'
'set z 3'
'define h=ave(go3,t=60,t=151)'
'close 1'
'sdfopen d:\ozone\go3.2011.nc'
'set z 3'
'define m=ave(go3,t=60,t=151)'
'close 1'
'sdfopen d:\ozone\go3.2012.nc'
'set z 3'
'define n=ave(go3,t=61,t=152)'
'define s=(a+b+c+d+e+f+g+h+m+n)*1000000/10'
'd s'
'cbar 1.0 0'
'disable contour'


然后我编了一个
'reinit'
'set gxout contour'
yy=2003
i=1
'a=0'
while(yy<=2012)
'sdfopen d:\ozone\go3.'yy'.nc'
'set z 1'
if ((math_mod(yy,4)=0&math_mod(yy,100)!=0)|math_mod(yy,400)=0)
'define a'i'=ave(go3,t=61,t=152)'
else
'define a'i'=ave(go3,t=60,t=151)'
endif
'a=a+a'i''
i=i+1
yy=yy+1
endwhile
'd a/10'
'cbar 1.0 0'
'disable contour'

但是运行出错，求助怎么改呀？

蓉蓉不着急

{:cry:}求大佬帮忙看看