grads 批处理1个月的数据，却只读入了第一个数据。恳请大家帮忙指正！

白兔糖0217

文件夹下有201401整整一个月的ctl 和 dat文件。文件如图：

现在我先对这些文件进行批处理读取，生成grd文件。下面是我的脚本。结果：只有第一个文件读取到了，后面的没有读取到，想知道是怎么回事？该做怎么样的修改？？坐等大家的指正！
'reinit'
dd=1
while(dd<=31)
if(dd<10)
    aaa=substr(20140101,1,7)
else
    aaa=substr(20140101,1,6)
endif
'open D:/Fortran/ARWpost/201401/'aaa''dd'.ctl'
'set gxout fwrite'
*'set fwrite /global/research/iumhouyy/grads/201401'dd'.grd'
'set lat 32.0'
'set lon 118.8'
it=1
while(it<5)
'set t 'it
'q time'
m=substr(result,13,3)
d=substr(result,11,2)
h=substr(result,8,2)
'set fwrite D:/Fortran/ARWpost/201401/'aaa''dd''h'.grd'
i=1
***** here,we need to change i ,this time ,we set max i is 58
while(i<35)
'set z 'i
'd height'
'd pressure'
'd tk'
'd tc'
'd td'
'd rh'
i=i+1
endwhile
'disable fwrite'
it=it+1
endwhile
dd=dd+1
endwhile
**grads2ascii result D:\Fortran\ARWpost\test3.txt %10.3f 7 1
'reinit'
;

白兔糖0217

啊  倒是提点意见啊

river

有报错信息吗？

白兔糖0217

river 发表于 2014-11-14 16:50
'open D:/Fortran/ARWpost/201401/'aaa''dd'.ctl'
'open D:/Fortran/ARWpost/201401'aaa''dd'.ct ...

D:/Fortran/ARWpost/201401/  是我的ctl路径。去掉"/",恐怕不对哦

river

白兔糖0217 发表于 2014-11-14 16:52
D:/Fortran/ARWpost/201401/  是我的ctl路径。去掉"/",恐怕不对哦

你提取出来的数据的ctl 和原始数据的ctl ？

白兔糖0217

river 发表于 2014-11-14 17:03
你提取出来的数据的ctl 和原始数据的ctl ？

嗯啊   都在../201401/目录下。我感觉问题在这边：
dd=1
while(dd<=31)
if(dd<10)   
aaa=substr(20140101,1,7)
else
     aaa=substr(20140101,1,6)
endif
'open D:/Fortran/ARWpost/201401/'aaa''dd'.ctl'

river

白兔糖0217 发表于 2014-11-14 17:07
嗯啊   都在../201401/目录下。我感觉问题在这边：
dd=1
while(dd

我觉得用笨办法就可以解决，不需要像substr那么高深的功能，可以这样
dd=1
while(dd<=31)
if(dd<10)
'open D:/Fortran/ARWpost/201401/2014010'dd'.ctl'

`````````
else
'open D:/Fortran/ARWpost/201401/201401'dd'.ctl'

`````````

endif

白兔糖0217

river 发表于 2014-11-14 21:24
我觉得用笨办法就可以解决，不需要像substr那么高深的功能，可以这样
dd=1
while(dd

我试试去 ，多谢您啦