求助grads提取多时间，多层次

天道酬勤

'reinit'
'sdfopen e:\osci\vapor\wind\uwnd.1976.nc'
'set fwrite e:\osci\vapor\1976uu.grd'
'set gxout fwrite'
‘set y 1 73’
'set x 1 145'
tt=122
while(tt<=274)
z=1
while (z<=12)
'set z 'z''
'set t 'tt''
'd uwnd'
z=z+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
;

dset e:\osci\vapor\1976uu.grdtitle paomian
undef -9.99e+08
XDEF 145 linear 0 2.5
YDEF 73 linear -90 2.5
ZDEF 12 levels 1000 925 850 700 600 500 400 300 250 200 150 100
TDEF 153 LINEAR 01may1976 1dy
VARS 1
ji  0 99 paomian
endvars
上面分别是提取的gs文件和ctl。我用某个特定的时间，层次和原nc文件做了核对，发现提取出来的图和原始的图不一致。到底是哪里出了问题呢， 求教。。。

曾记否

tt=122
while(tt<=274)

TDEF 153 LINEAR 01may1976 1dy
你的gs和ctl时间怎么都没对应？让tt从1开始循环试试！

jeffggg

u场不是17层吗 你只提取了其中的12层 是不是是这里出了问题？而且时间貌似不太对应啊

天道酬勤

曾记否 发表于 2013-4-14 15:05

                               
登录/注册后可看大图

tt=122
while(tt

我是从5月1号开始提取的。总共153天。

天道酬勤

jeffggg 发表于 2013-4-14 15:05

                               
登录/注册后可看大图

u场不是17层吗 你只提取了其中的12层 是不是是这里出了问题？而且时间貌似不太对应啊

因为另一个变量只有12层，所以为了以后的工作，我就只提了12层的数据。时间是对的啊，不知道哪里出了问题。

曾记否

天道酬勤 发表于 2013-4-14 15:22

                               
登录/注册后可看大图

我是从5月1号开始提取的。总共153天。

但你的ctl中时间变量只有153，你的gs中t却到了274，这样写应该不对啊

天道酬勤

找到问题了，很简单的一个问题，ctl倒数第二行0应该变成12！

river

VARS 1
ji  0 99 paomian
endvars
红色那一行有问题，你提取了12层，怎么写了0 99，那就是告诉grads这个变量只有地面一层，应该写成12 99

天道酬勤

river 发表于 2013-4-14 18:54

                               
登录/注册后可看大图

VARS 1
ji  0 99 paomian
endvars

恩，我已经找出来了。浪费了整整一天。。。谢谢啊！

river

天道酬勤 发表于 2013-4-14 19:29

                               
登录/注册后可看大图

恩，我已经找出来了。浪费了整整一天。。。谢谢啊！

找到错误就好，我也是回复完了才发现楼主已经解决了呢。楼主编辑一下帖子的名字吧，添加【已解决】