systemfunc批量读取站点文件转为nc文件出现问题

又又又加零

我用NCL中的systemfunc函数批量读取1961-2016的数据，出现了只能读取2007-2016年的数据，并且报错说数组超过上限（已经解决）                            



代码参考的是《NCL图形分析语言入门到精通》PDF文件的站点数据处理为nc文件，我把主程序读取的代码放在下面了。

1. ;-----------------------------------------主程序---------------------------------------------------
   
2. 
   
3. begin
   
4. 
   
5. 
   
6. ;----选取最大的文件获取所有站点信息
   
7. bigfile = "E:/DATA/196101-201612/SURF_CLI_CHN_MUL_DAY-PRE-13011-201606.TXT"
   
8.   stn_info = asciiread(bigfile,(/25170,13/),"integer")           ;25170行x13字段
   
9.   stn=stn_info(::30,0)                                           ;各站间隔30日(6月的数据)
   
10.   stnnum =dimsizes(stn)                                          ;站数(实际为839站)
    
11. 
    
12. 
    
13. ;----经纬度数据格式为:度分(需转换为度)
    
14. lats=stn_info(::30,1)/100+mod(stn_info(::30,1),100)/60.
    
15. lons=stn_info(::30,2)/100+mod(stn_info(::30,2),100)/60.
    
16. 
    
17. 
    
18. ;----高度单位0.1m和估测站点高度处理
    
19. height=stn_info(::30,3)*0.1
    
20. height=where(height.gt.10000.,height-10000.,height)
    
21. delete(stn_info)                                                        ;释放内存
    
22. sdate = 1951
    
23. edate = 2016
    
24. ndate=yyyymmdd_time(sdate,edate,"integer")
    
25. 
    
26. 
    
27. ;----生成时间段内每一天的数组,如 year*10000+month*100+day
    
28. ndays     = dimsizes(ndate)                                             ;总天数
    
29. data      = new((/stnnum,ndays/),"integer")                             ;定义变量数组
    
30. data      = 32766                                                       ;数组初始赋值缺省值
    
31. 
    
32. 
    
33. ;----------------------------------读取数据文件操作------------------------------------------------
    
34. 
    
35. files     = systemfunc("ls E:/AAAJDJS/AAA196101-201612/*.TXT")                         ;所有数据文件名列表
    
36. nf        = dimsizes(files)                                                            ;列表中文件数
    
37. cols      = 13                                                                         ;13个数据字段
    
38. pos       = 0
    
39. 
    
40. do i = 0,nf-1
    
41.    print(i+"Nowreading..."+files(i))                                 ;过程监控
    
42.    dummy                      = asciiread(files(i),-1,"integer")
    
43.    rows                       = dimsizes(dummy)/cols
    
44.    nfile                      = onedtond(dummy,(/rows,cols/))        ;转换为二维行x列
    
45.    ntm                        = max(nfile(:,6))                      ;天数
    
46.    stnf                       = nfile(::ntm,0)                       ;站号
    
47.    nstnf                      = dimsizes(stnf)                       ;站数
    
48.    idxstn                     = get1Dindex(stn,stnf)                 ;站索引位置
    
49.    temp                       = reshape(nfile(:,9),(/nstnf,ntm/))    ;匹配数组:(站,时间)
    
50.    data(idxstn,pos:pos+ntm-1) = temp(:,:)                            ;分段赋值累计降水量
    
51.    pos                        = pos + ntm                            ;移动时间维位置指针
    
52.    delete(dummy)                                                     ;注意清除循环内重复赋值变量
    
53.    delete(nfile)
    
54.    delete(rows)
    
55.    delete(stnf)
    
56.    delete(temp)
    
57.    delete(idxstn)
    
58. 
    
59. end do

复制代码

又又又加零

解决了，是因为07年以前的数据后缀是小写