grads 批处理文件的问题

笨笨兔

      我的数据文件是2014122101.grd，2014122102.grd，2014122103.grd------2014122123.grd  .一共23个文件。我想把它们批处理合并到一个grd中
      我的ctl文件如下
DSET E:\data\testd3.20141221%h2.grd
TITLE NAQMPS output for domaind3
OPTIONS byteswapped  template
UNDEF     -1.E20
XDEF 204 LINEAR 1. 1.0
YDEF 204 LINEAR 1. 1.0
ZDEF   12 LEVELS
   50.00000       160.0000       300.0000       470.0000       680.0001   
   929.9999       1220.000       1570.000       1990.000       2500.000   
   3110.000       3840.000   
TDEF      23 LINEAR 01:00Z21DEC2014 60MN
VARS     63
dz         12  99 Layer Depth m
u          12  99 U         m/s
v          12  99 V         m/s
w          12  99 V         m/s
temp       12  99 temp      K
............


gs文件：

'reinit'
'open E:\data\testd3.2014122101.ctl'
'set gxout fwrite'
'set fwrite E:\data\outdata.grd'
'set grads off'
hour=1
t=1
while(hour<24)
'set x 1 204'
'set y 1 204'
'set t 't
'set z 1'
'd rh'
hour=hour+1
t=t+1
endwhile
'disable fwrite'
;

只能生成一个时次的图（只读入了t=1一个时次的数据），t=2以后都报错。。

如果固定了x,y就可以输出结果。。
不知道是哪里出了问题，想请教大家帮忙看看。{:soso_e134:}

river

你的循环有问题，时间根本没有循环······

笨笨兔

river 发表于 2015-1-6 07:55
你的循环有问题，时间根本没有循环······

是t的问题么？能不能大概说一下我要怎么改时间呢，试了很多次都没弄明白 。。

维克多·蓝晶

你只最开始打开了一个ctl，而且z明明有12层，'set t 't是什么东西？你原始grd里面都只有1个时次对吧？那除了set t 1还能选别的？


while(hour<10)
'open E:\data\testd3.201412210'hour'.ctl'
'set x 1 204'
'set y 1 204'
'set t 1'
'set z 1 12'
'd rh'
hour=hour+1
'close 1'
endwhile
while(hour<23)
'open E:\data\testd3.20141221'hour'.ctl'
'set x 1 204'
'set y 1 204'
'set t 1'
'set z 1 12'
'd rh'
hour=hour+1
'close 1'
endwhile

ha33ha33

太好了，就是需要这个

笨笨兔

维克多·蓝晶 发表于 2015-1-6 09:07
你只最开始打开了一个ctl，而且z明明有12层，'set t 't是什么东西？你原始grd里面都只有1个时次对吧？那除 ...

因为z我是想着只提取一层出来，所以当时设的‘set z 1’，这种批处理是不可以只提取一层么？那我再去把时间那块改改试试，grads不太熟练，谢谢你的指正

笨笨兔

维克多·蓝晶 发表于 2015-1-6 09:07
你只最开始打开了一个ctl，而且z明明有12层，'set t 't是什么东西？你原始grd里面都只有1个时次对吧？那除 ...

我的打开的ctl文件应该不只一个呀？那个ctl就是我上面贴出来得描述文件，里面DSET E:\data\testd3.20141221%h2.grd不是应该就表示一直读入所有23个grd么？我是所有grd都用这一个ctl来读取，虽然名字是testd3.2014122101.ctl，但是不是只读的testd3.2014122101.grd。

笨笨兔

river 发表于 2015-1-6 07:55
你的循环有问题，时间根本没有循环······

我把时间hour=1
t=1
while(hour<24)
'set x 1 204'
'set y 1 204'
'set t 't
'set z 1'
'd rh'
hour=hour+1
t=t+1
endwhile
改为

hour=1
while(hour<24)
'set x 1 204'
'set y 1 204'
'set t 1'
'set z 1'
'd rh'
hour=hour+1
endwhile
还是出错。。是hour不对么？
报错提示：Low level I/O Error: read error on data file
                 data filename=E:\data\outdata.grd
                Error reading 204 bytes at location 3995136
               Data Request Error:Error for variable'rh'
                 Error ocurred at column 1
               DISPLAY error: Invalid expression
                 Expression=rh
不太懂这是哪里的问题。。。

river

笨笨兔 发表于 2015-1-6 09:53
我把时间hour=1
t=1
while(hour

没有改到点子上，而且t是grads里面默认的时间，所以不能随便拿来当变量用，像z x y z也是一样的
hour=1
while(hour<24)
'set x 1 204'
'set y 1 204'
'set t 'hour''
'set z 1'
'd rh'
hour=hour+1
endwhile

路易十_W8L8f

像上面那样改，问题就没有了吗？