Grads程序求助

曜城

        提取的1979-2012年逐日的850hPa经向风资料。现在要对其做5点平滑，去掉5天内的信号，从1979年01月01日到2012年12月31日，共12419天。可用Grads程序编出来的，发现只有12418天。程序如下，请各位高手指教！
'reinit'
'open D:/Data/Ncep2/1xdaily/v850/v850.1979-2012.ctl'
'set gxout fwrite'
'set fwrite D:/Data/Ncep2/1xdaily/v850/v850.1979-2012.5sm.grd'


'set x 1 144'
'set y 1 73'

'd v850(t=1)'
'd v850(t=2)'

i=3
while(i<=12417)

j1=i-2
j2=i-1
j4=i+1
j5=i+2

'define v8502=(-3*v850(t='j1')+12*v850(t='j2')+17*v850(t='i')+12*v850(t='j4')-3*v850(t='j5'))/35.0'
'd v8502'
say i
i=i+1

endwhile

'd v850(t=12418)'
'd v850(t=12419)'

destroyer

确定只有12418天？用fortran读一下层数试试？

destroyer

确定只有12418天？用fortran读一下层数试试？

曜城

destroyer 发表于 2013-12-26 11:04

                               
登录/注册后可看大图

确定只有12418天？用fortran读一下层数试试？

谢谢，已经解决了，最后一行没有打回车，grads运行的时候就没有运行最后一行，少了最后一个时次。

若兰

grads里面最后最后加上reinit命令，然后打几个回车。