grads读不同的文件，生成一个grd数据文件的方法

笛声悠悠

楼主下了海洋垂直速度dzdt资料，但不是打包下载的，每一年一个单独的文件，现在要用grads的fwrite命令把每个文件中的dzdt提取出来，相信很多同学都遇到过类似的问题，下面交流一下我的经验：
'reinit'
'set fwrite d:/dzdt.grd'
t1=1979                               *开始的年份
t2=1                                     *记录打开文件的设备号
while(t1<=2000)                   *时间的循环
'sdfopen e:\data\dzdt\dzdt.'%t1%'.nc'  *%里是文件名的时间
'set gxout fwrite'
tt=1                                           *月份的循环
while(tt<13)
zz=1                                        *截取的层数
while(zz<41)                            *层数循环
'set lon 40.5 120.5'
'set lat -29.5 29.5'
'set t 'tt
'set z 'zz
'd dzdt.'t2'(t='tt',z='zz')'
zz=zz+1
endwhile
tt=tt+1
endwhile
t1=t1+1
t2=t2+1
endwhile
'disable fwrite'
;

日絮宝贝

十分感谢，看你的帖子，搞定问题了