ubscripts must be integer when not using coordinate indexing

sususue

数据包括全国很多站点，按一个站点一个经纬度，Nstp个数据读取，现在想要把得到的二维数组clod(站点号，Nstp)存在一个三维数组里面，来得到全国的分布

1. do i =0,ndim1-1               
   
2.                  
   
3.   t     = asciiread( path1+files1(i),-1,"string")                
   
4.   nrows = dimsizes(t)               
   
5.   yr    = stringtofloat(str_get_field(t(nrows-1), 9,   "        "))
   
6. 
   
7.   if ( yr .eq. 2014) then
   
8.   delim = str_get_tab()
   
9.   id(mm)   = stringtofloat(str_get_field(t(nrows-1), 1,  delim))
   
10.   lat(mm)   = stringtofloat(str_get_field(t(nrows-1), 3, delim))
    
11.   lon(mm)   = stringtofloat(str_get_field(t(nrows-1), 4, delim))
    
12.   clod(mm,:) = stringtofloat(str_get_field(t(nrows-Nstp:nrows-1),27, delim))
    
13. 
    
14.         mm = mm+1
    
15. 
    
16.   end if               
    
17.                
    
18. delete(t)               
    
19. delete(nrows)               
    
20. delete(yr)               
    
21. 
    
22. end do
    
23. 
    
24.        
    
25. 
    
26. kk=addfile("/disk02/hg-ws/data40G/1.nc","r")
    
27. aa=kk->ts
    
28. printVarSummary(aa)
    
29. 
    
30. 
    
31. M=new((/360,180,Nstp/),"float",aa@_FillValue)
    
32. M!0 ="lat"
    
33. M@lat= aa@lat
    
34. M!1 ="lon"
    
35. M@lon= aa@lon
    
36. M!2 ="lat"
    
37. M@time=  ispan(1,Nstp,1)
    
38. 
    
39. 
    
40. lat=round(lat,0)
    
41. lon=round(lon,0)
    
42. 
    
43. 
    
44. do tt=0,Nstp-1
    
45. do i=0,mm-1       
    
46.         <font color="Red">M(lat(i),lon(i),tt) = clod(i,tt)</font>
    
47. end do
    
48. end do
    
49. kk1=addfile("/data40G/test.nc","c")
    
50. kk1 ->clod =M
    

复制代码

但是提示出错为
我查了一下好像是因为我的经纬度没有按照单调递增这样有序存入，是零散分布在全国无序的。
请问大神们这个要怎么处理？