画轨迹报错!

对方的凤飞飞

1. 
   
2. begin
   
3. 
   
4.    diri="/media/wjj/data/FLEXPART/"
   
5.    files = systemfunc ("ls " + diri + "*.txt")
   
6.    ntime=dimsizes(files)
   
7.    loctn=new((/ntime,39,2/),float,-999)
   
8.    do i=0, ntime-1
   
9.       print(" "+files(i))
   
10.       ;data= asciiread(files(i),(/39,13/),"float")
    
11.       nrow=numAsciiRow(files(i));读取行
    
12.       ncol=numAsciiRow(files(i))
    
13.       lines=asciiread(files(i),(/nrow,ncol/),"float")
    
14.       loctn(i,0:nrow-1,0)=lines(0:nrow-1,2);经度
    
15.       loctn(i,0:nrow-1,1)=lines(0:nrow-1,3);纬度
    
16.       delete(lines)
    
17.    end do
    
18.    do i=0, ntime-1
    
19.    wks = gsn_open_wks("png","traj")               ; send graphics to PNG file
    
20.   
    
21. 
    
22.    
    
23.    res             = True                              ; plot mods desired
    
24.    res@gsnFrame    = False                             ; don't advance yet
    
25.    res@gsnMaximize = True
    
26.    res@gsnDraw          = False   ; so we can add poly stuff
    
27.    res@mpDataBaseVersion     =  "MediumRes" ;
    
28.    res@mpPerimOn             = False          ;-- don't draw the box around the plot
    
29.    res@gsnPolar="SH"
    
30.    res@mpCenterLonF          = 0              ;-- center at lon=10
    
31.    res@mpMaxLatF          = -30   
    
32.    res@mpGridLineColor       = "grey30"         ;-- set grid line color
    
33.    res@tiMainString    = ""  ;
    
34.    res@gsnCenterString = ""   ; center string
    
35.    res@mpFillOn               = False    ;
    
36.    res@mpOutlineOn            = True    ; turn on continental outline
    
37. 
    
38. 
    
39.    map = gsn_csm_map(wks,res)                        ; gsn_csm_map_ce
    
40.    pres                  = True               ; polyline resources
    
41.    pres@gsLineThicknessF = 2.0                ; line thickness
    
42.   
    
43.    pres@gsLineColor = "red"                        
    
44.    gsn_polyline(wks,map,loctn(i,:,0),loctn(i,:,1),pres)      ; draw the traj
    
45.    draw(map)                                          
    
46.    frame(wks)
    
47.    end do                                       
    
48.   
    
49. end
    

复制代码

代码写成这样,还是报错,求指教!
warning:asciiread: End of file reached and only (507) elements were read from the file, filling remaining elements with the default missing value for the requested type