如何提取年份并写入一个grd文件中？请问我写的gs和ctl哪里出错了？

LemonTree007

提取年份的gs 如下：

1. 'reinit'
   
2. 'set gxout fwrite'
   
3. 'set fwrite D:\matlabht\tpdsi-spei\yt.grd'
   
4. 'open D:\matlabht\pdsi\pdsi1901-2014.ctl'
   
5. tt=1
   
6. while(tt<115)
   
7.     'set t 'tt
   
8.     'q time'
   
9.     res=subwrd(result,3)
   
10.     res=substr(res,9,12)          **00z01jan1901-00z01jan2014
    
11.     'd 'res
    
12.     tt = tt+1
    
13. endwhile
    
14. 'disable fwrite'
    
15. 'reinit'

复制代码

运行之后顺利得出grd文件，但是写了ctl却打不开，不知道哪里错了，ctl如下：

1. dset d:\matlabht\tpdsi-spei\yt.grd
   
2. undef 9999.9
   
3. title  1901-2014 year time
   
4. tdef 114 linear 1901 12mo
   
5. vars  1
   
6. res  0 99 year time
   
7. endvars

复制代码

LemonTree007

请忽略‘set gxout fwrite’,这个命令是用来写入grd文件的，对此问题没用，忘记删除了