写了一个批处理ctl和gs

白兔糖0217

自己写了一个批处理的ctl和gs脚本。
一边请大家帮帮指正错误，顺便分享一下。
ps: 文件路径如图： 这些ctl只是做测试用。
批处理的ctl:
dset D:\Fortran\ARWpost\test\2014%m2%d2.dat
title data in grads format
options template
undef 1.e30
xdef  121 linear  116.33591   0.04054054
ydef  161 linear   28.73539   0.04054054
zdef   20 levels （详细高度：lve）
tdef    3 linear 01JUL2014     1day  

gs脚本：
'reinit'
'open D:/Fortran/ARWpost/test/all_ctl.ctl'                       *打开批处理ctl
year=2014                                                                  * 批处理开始的年份
month=07                                                                               *批处理开始的月份
t=1                                                                                           **从第几个时间序列开始
while(month<13)
mm=%month
'set gxout fwrite'
'set lat 35.0'
'set lon 118.8'
y=2014
it=1
while(it<5)
'set t 'it
'q time'
h=substr(result,8,2)
'set fwrite D:\Fortran\ARWpost\test\'%year''mm''h'.grd''
i=1
while(i<48)
'set z 'i
'd pressure'
'd tk'
'd tc'
'd td'
'd rh'
i=i+1
endwhile
'disable fwrite'
it=it+1
endwhile
month=month+1
t=t+1
endwhile
'reinit'
;


运行时，老说第二行有问题，ctl写错了？？？

白兔糖0217

草包问题，难倒寡人了

白兔糖0217

重新写了一个，'reinit'
yy=2014
mm=07
dd=01
while(dd<32)
'open D:\Fortran\ARWpost\test\'yy''mm''dd'.ctl'      


这个的写法不对？为什么呢