某一格式数据用grads批处理转化为二进制数据

单音和弦

各位好，我下载的逐日数据格式有点奇怪，但是网站配有ctl可以打开，每一个逐日数据是一个文件，数据是V1.0_0.50deg.lnx.19790501，V1.0_0.50deg.lnx.19790502，V1.0_0.50deg.lnx.19790503。。。我现在想用让它们共用一个crl文件把它们合并到一个grd文件中，（将19790501—19790930这153个数据合并）。但是合并后的数据画某一天的图和原始数据不一样，只有第一天是一样的，后面都不一样，不知道是什么问题？？请教各位！！！！谢谢！！！
ctl 如下：
dset e:\1979\V1.0_0.50deg.lnx.1979%m2%d2
options template
title  Pre
undef -999.0
xdef 720 linear   0.25 0.50
ydef 360 linear -89.75 0.50
zdef 1   linear 1 1
tdef 153 linear 01May1979 1dy
vars 2
r    1  00 the grid analysis (0.1mm/day)
gnum     1  00 the number of stn
ENDVARS

gs 如下：
'reinit'
'open E:\1979\Pre1979.ctl'
'set gxout fwrite'
'set fwrite E:\1979\1979_56789.grd'
'set grads off'
tn=1
while(tn<=153)
'set x 1 720'
'set y 1 360'
'set z 1'
'set t 'tn''
'd r'
tn=tn+1
endwhile
'disable fwrite'
''
;
合并后的数据的ctl:
dset e:\1979\1979_56789.grd
title Pre
undef -999.0
xdef 720 linear 0.25 0.50
ydef 360 linear -89.75 0.50
zdef 1   linear 1 1
tdef 153 linear 01May1979 1dy
vars 2
r    1  00 the grid analysis (0.1mm/day)
gnum     1  00 the number of stn  
ENDVARS
原始数据画有一个的图：

用合并后的ctl文件也画这天的图：

单音和弦

原始数据的图

合并后数据的图
终于上图了

单音和弦

为什么总是不能上图呢 反正两个图出来不一样，请教各位！！

river

都已经批量描述了，然后再合并成一个二进制文件不是多余的么

单音和弦

river 发表于 2015-4-23 17:32
都已经批量描述了，然后再合并成一个二进制文件不是多余的么

因为这个数据不是二进制的呀，想合并生成一个二进制的文件然后fortran可以读了。写这个合并二进制文件的ctl是为了和原始数据验证一下，不知道为什么会不一样呢？？？

单音和弦

没有人知道吗？请大家帮忙看看！！