求南海夏季风指数程序出错

hgsndg

'reinit'
'sdfopen E:\aa\uwnd.1987.nc'
'sdfopen E:\aa\vwnd.1987.nc'
'set grads off'
'set lat 0 10'
'set lon 100 130'
'set lev 850'
'define u8=uwnd.1'
'define v8=vwnd.2'
'set dfile 1'
'set lat 0 10'
'set lon 100 130'
'set lev 200'
'define u2=uwnd'
'set dfile 2'
'set lat 0 10'
'set lon 100 130'
'set lev 200'
'define v2=vwnd'
'define aa=aave(hdivg(u8,v8)-hdivg(u2,v2),lon=100,lon=130,lat=0,lat=10)'
'define bb=hdivg(u8,v8)-hdivg(u2,v2)'
'define cc=pow(bb,2)'
'define dd=asumg(cc,x=1,x=13,y=1,y=5)'
'define ee=sqrt(dd/65)'
'define udif=(bb/65)/ee'
'd udif'
'define r=tcorr(udif*100000,v2,t=160,t=190)'
'set gxout shaded'
'd r'
;
错误应该是出在求南海夏季风指数那里。计算公式为：Id= ΔD/sqrt( ∑(ΔDj)2/ n)
数来研究南海夏季风的活动与大气环流的关
系,其中ΔD=(D850- D200),这里的“-”
表示空间平均,Di是某格点的散度,j是时间
序列。   但是不知道该怎么写，谢谢大家帮忙

chenxiaoyun

请问有改好吗？求分享