grads提取数据成grd文件时出现错误

旺旺小小酥

将模式的结果提取某几个变量到grd文件，写了grd文件的描述文件ctl之后发现提取出来的数据都是缺省值，不知道为什么

这是提取数据的脚本
'reinit'
'open /test/xkxu/daymean_data/50km/sub50ATMALL.ctl'
'set fwrite pre1Q1.grd'
'set gxout fwrite'
'define t=ta'
'define u=ua'
'define v=va'
'define w=omega*(-10)'
i=1
while(i<=92)
'set t 'i''
p=1
while(p<=12)
'set z 'p''
'set x 1 228'
'set y 1 119'
'd t'
'd u'
'd v'
'd w'
p=p+1
endwhile
i=i+1
endwhile
'disable fwrite'

这是编写的描述提取grd文件的ctl

dset   ^pre1Q1.grd
title  original data
undef  1e+20_FillValue
xdef 228 linear 58.48 0.45
ydef 119 linear 6.49  0.45
zdef 12 levels 1000.0  925.0  850.0  700.0  600.0  500.0  400.0  300.0  250.0  200.0 150.0  100.0
tdef 92 linear 01Jun2000  1dy
vars 4
t=>t 12 t,z,y,x  Air Temperature (K)
u=>u 12 t,z,y,x  Grid Eastward Wind (m s-1)
v=>v 12 t,z,y,x  Grid Northward Wind (m s-1)
w=>w 12 t,z,y,x  vertical velocity (m s-1)
ENDVARS

这是open pre1Q1.ctl之后随便d 一个变量的样子

[img]file:///C:\Users\47257\AppData\Roaming\Tencent\Users\472578644\QQ\WinTemp\RichOle\%5MD{UU(QWRE1PT~B0C)V@G.png[/img]

旺旺小小酥

问题已解决，最内层的循环不能全部写到一起，要拆开，一个气压层的循环，输出一个变量，下面是正确的脚本，希望对大家有帮助
'reinit'
'open /test/xkxu/daymean_data/50km/sub50ATMALL.ctl'
'set fwrite pre1Q1.grd'
'set gxout fwrite'
i=3288
while(i<=3379)
'set t 'i''
p=1
while(p<=12)
'set x 97.7111 148.822'
'set y 37.6889 71.0222'
'set z 'p''
'define t=ta'
'd t'
p=p+1
endwhile

n=1
while(n<=12)
'set x 97.7111 148.822'
'set y 37.6889 71.0222'
'set z 'n''
'define u=ua'
'd u'
n=n+1
endwhile

m=1
while(m<=12)
'set x 97.7111 148.822'
'set y 37.6889 71.0222'
'set z 'm''
'define v=va'
'd v'
m=m+1
endwhile

l=1
while(l<=12)
'set x 97.7111 148.822'
'set y 37.6889 71.0222'
'set z 'l''
'define w=omega*100'
'd w'
l=l+1
endwhile

i=i+1
endwhile
'disable fwrite'

bright

非常非常感谢