while loop problem

[Mike Ray]

I am trying to estimate the average of the square of the anomalies,
multiplied by cos(theta). theta=latitude.
can someone find the mistake!

'open u.ctl'
'open v.ctl'
'define ubar=ave(u,t=1,t=977)'
'define vbar=ave(v.2,t=1,t=977)'
t1=1
sumud2=0
while(t1<978)
'set t 't1
'set lon 40 180'
'set lat -10 10'
'set lev 850'
'define ud2=(u-ubar)*(u-ubar)'
'define sumud2=sumud2+ud2'
t1=t1+1
endwhile
'd  (sumud2/977)*cos(lat*22/7/180)'

_________________________________________________________________
Is your PC infected? Get a FREE online computer virus scan from McAfee®
Security. http://clinic.mcafee.com/clinic/ibuy/campaign.asp?cid=3963