循环提取冬季的数据

luoziwuhui

想利用循环提取1979-2010年冬季的数据，原始的数据是每年放在一个文件里的（1979-2011年）
这就涉及到一个跨年飞问题，下面是自己编的程序。
前一年的12月能提取，没问题。到了后一年的1-2月份，就出错，说超出了范围{:soso_e127:}。
请教一下各位！
谢谢{:soso_e113:}！


附程序：
************************************
'reinit'
iyr=1979
while(iyr<=2010)
iyr2=iyr+1
'open E:\tj\pre\nizi\rain'iyr'.ctl'（前一年）
'open E:\tj\pre\nizi\rain'iyr2'.ctl'（后一年）

'set gxout fwrite'
'set fwrite E:\tj\pre\nizi\winter.'iyr'.dat'
********前一年12月******************************************
'set time 01dec'iyr' 31dec'iyr''
'set x 1 142'
'set y 1 82'
'set z 1'
'd rain'

'set dfile 2'
***********后一年1-2月****************************************
'set time 01djan'iyr2' 28feb'iyr2''
'set z 1'
'set x 1 142'
'set y 1 82'
'd rain'
iyr=iyr+1
'disable fwrite'
endwhile

lotina

我没用过SET DFILE的命令
这个不需要关闭第一个文件吗？
一般来说再循环里打开一个文件就要关闭一个文件啊。

mofangbao

如果两个ctl的时间不一致还是一个个的打开 关闭比较保险

凡不凡

正在学grads