提取垂直剖面数据问题【已解决】

aa22556688

资料是2.5*2.5的omega数据， z:1~12; t:1~795, 固定x=35，提取经向剖面图（y:21~53），然后对时间序列的剖面进行eof分析。

gs文件如下：
'reinit'
'sdfopen E:\Data\omega.mon.mean.nc'
'set fwrite E:\ZC_test\omega1.grd'
'set gxout fwrite'
tt=1
while(tt<=795)
'set t 'tt''
zz=1
while(zz<=12)
'set z 'zz''
'set x 35'
'set y 21 53'
'd omega'
zz=zz+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
；
提取数据的ctl：
dset E:\ZC_test\omega1.grd
undef -9.99e+33
title omega
xdef 1 levels 85.5
ydef 33 linear -40.0 2.5
zdef 12 levels 1000 925 850 700 600 500 400 300 250 200 150 100
tdef 795 linear 00JAN1948 1mo
vars 1
omega 12 99 Vertical velocity  *** 此处维数信息有误，现在是正确的。
endvars

但是用grads画图后，提示z维有问题，不知道哪里错了， 按照这个顺序提取数据，应该没有问题啊，求大神给指点迷津，非常感谢！！！

river

xdef 1 levels 85.5
ydef 33 linear -40.0 2.5
红色部分为什么是那样的呢？明明它的格式应该和ydef那一句是一样的。

lqouc

修改ctl如下
vars 1
omega 12 99 ertical velocity
endvars

aa22556688

river 发表于 2014-10-20 07:23
xdef 1 levels 85.5
ydef 33 linear -40.0 2.5
红色部分为什么是那样的呢？明明它的格式应该和ydef那一 ...

因为是沿着y方向的剖面图，所以y只有一层

aa22556688

lqouc 发表于 2014-10-20 08:07
修改ctl如下
vars 1
omega 12 99 ertical velocity

非常感谢，确实是这的问题，忽略了层次，再次感谢。

小预报员

真的吗？！！！！！！！！！！！！！！

aa22556688

小预报员 发表于 2014-10-20 21:28
真的吗？！！！！！！！！！！！！！！

上面的都是正确的，试衣镜修改好了的。