grads计算整层水汽输送通量

大果子

各位大神，这是我编写的计算整层水汽输送通量的脚本，运行时总出现Data request warning：Request is completly outside file limits。谁能帮我看看怎么回事？不胜感激！


'reinit'

y=1981
while(y<=2014)
'sdfopen k:\nc\air.'y'.nc'
'sdfopen k:\nc\rhum.'y'.nc'
'sdfopen k:\nc\uwnd.'y'.nc'
'sdfopen k:\nc\pres.sfc.'y'.nc'

'set fwrite k:\qj\qu'y'.grd'

'set gxout fwrite'

nt=1
while(nt<=365)
'set t 'nt''
'set lon 0 160'
'set lat -20 70'

nz=1
while(nz<=8)
'set z 'nz''
'define es=6.112*exp(17.269*(air-273.15)/(air-35.86))'
'define qs=0.622*es/(lev-0.378*es)'
'define q=rhum.2*qs/100'
'define qu=uwnd.3*q/9.8'

nz=nz+1
endwhile

'set lon 0 160'
'set lat -20 70'
'set z 1'
'set t 'nt''
'aqu=vint(pres.4,qu,300)'
'd aqu'

nt=nt+1
endwhile
‘close 4’
'close 3'
'close 2'
'close 1'
'disable fwrite'
y=y+1
endwhile
;

Eleven

可能是nc文件的时间范围不在设置的范围内。