计算水汽通量时运行错误，原因是？

好猫

按照默小陌的http://bbs.06climate.com/forum.php?mod=viewthread&tid=7290进行到第一步便出现了问题，地面气压数据莫名其妙输不出来，gs如下：
'reinit'
'sdfopen e:\ncep\uwnd.mon.mean.nc'
'sdfopen e:\ncep\vwnd.mon.mean.nc'
'sdfopen e:\ncep\shum.mon.mean.nc'
'sdfopen e:\ncep\pres.mon.mean.nc'
'set gxout fwrite'
'set fwrite e:\ncep\shuiqi\water_high.dat'

'set z 1 6'
'set t 1 12'
'have=ave(uwnd,t+360,t=744,12)'
'have2=ave(vwnd.2,t+360,t=744,12)'
'have3=ave(shum.3,t+360,t=744,12)'

'modify have seasonal'
'modify have2 seasonal'
'modify have3 seasonal'

'set t 1'
'define djfu=(have(t=3)+have(t=4)+have(t=5))/3'
'define djfv=(have2(t=3)+have2(t=4)+have2(t=5))/3'
'define djfs=(have3(t=3)+have3(t=4)+have3(t=5))/3'

'set t 1'
iz=1
while(iz<=6)
'set z 'iz''
'd ave(have,t=423,t=425)-djfu'
'd ave(have,t=471,t=473)-djfu'
'd ave(have,t=567,t=569)-djfu'
'd ave(have,t=603,t=605)-djfu'
'd ave(have,t=663,t=665)-djfu'
'd ave(have,t=711,t=713)-djfu'
iz=iz+1
endwhile
'set t 1'
iz=1
while(iz<=6)
'set z 'iz''
'd ave(have2,t=423,t=425)-djfv'
'd ave(have2,t=471,t=473)-djfv'
'd ave(have2,t=567,t=569)-djfv'
'd ave(have2,t=603,t=605)-djfv'
'd ave(have2,t=663,t=665)-djfv'
'd ave(have2,t=711,t=713)-djfv'
iz=iz+1
endwhile
'set t 1'
iz=1
while(iz<=6)
'set z 'iz''
'd ave(have3,t=423,t=425)-djfs'
'd ave(have3,t=471,t=473)-djfs'
'd ave(have3,t=567,t=569)-djfs'
'd ave(have3,t=603,t=605)-djfs'
'd ave(have3,t=663,t=665)-djfs'
'd ave(have3,t=711,t=713)-djfs'
iz=iz+1
endwhile

'set t 1 12'
'have4=ave(pres.4,t+360,t=744,12)'
'modify have4 seasonal'
'set t 1'
'define djfp=(have4(t=3)+have4(t=4)+have4(t=5))/3'
'd ave(have4,t=423,t=425)-djfp'
'd ave(have4,t=471,t=473)-djfp'
'd ave(have4,t=567,t=569)-djfp'
'd ave(have4,t=603,t=605)-djfp'
'd ave(have4,t=663,t=665)-djfp'
'd ave(have4,t=711,t=713)-djfp'

'disable fwrite'
'reinit'


只输出前面的比湿和风完全没有问题，再加上输出地面气压的时候，grads就显示了下面的错误，然后grads就自动关掉了，这还是我手快截下的图= =
希望知道原因的童鞋给解释一下！拜托~~

好猫

额，下午来到仔细检查了下，把
‘have4=ave(pres.4,t+360,t=744,12)'这句改成了
'have4=ave(pres.4(z=1),t+360,t=744,12)'
就OK了……粗心大意真害人啊，留给以后的童鞋们当反面教材啦

qucc

谢谢，来学习一下

Olajuwon

入门级选手   看着很晕

更大华

pre(z=1)是必须的，要仔细啊！