关于disgrid插值函数

奇迹一哥

这是我写的代码，不知道为什魔输出结果是这样

1. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
4. 
   
5. begin
   
6. 
   
7. filepath = "/home/zxy/work/20100101.txt"
   
8. lines = asciiread(filepath,-1,"string")
   
9. 
   
10. ;delete the line which R is missed
    
11. e=dimsizes(lines)
    
12. 
    
13. 
    
14. station = stringtofloat(str_get_cols(lines,1,5))
    
15. lon = stringtofloat(str_get_cols(lines,16,21))
    
16. lat = stringtofloat(str_get_cols(lines,9,13))
    
17. pre  = stringtofloat(str_get_cols(lines,38,42))
    
18. 
    
19. olon = new((/491/),float)
    
20. olat = new((/411/),float)
    
21. data1 = new((/491,411/),float)
    
22. 
    
23. 
    
24. 
    
25. olon = fspan(67.25,137.37,491)
    
26. olat = fspan(15.19,55.23,411)
    
27. 
    
28. olon!0 = "lon"
    
29. olon@long_name  = "lon"
    
30. olon@units      = "degrees-east"
    
31. olon&lon        = olon
    
32. olat!0          = "lat"
    
33. olat@long_name  = "lat"
    
34. olat@units      = "degrees_north"
    
35. olat&lat        = olat
    
36. 
    
37. ;R@_FillValue = 999.0
    
38. dssetp("exp",0.5)
    
39. data1 = dsgrid2(lon,lat,pre,olon,olat)
    
40. print(data1)
    
41. ;print(data1(0,2))
    
42. opt = True
    
43. opt@fout = "/home/zxy/work/out_dsgrid2_2.txt"
    
44. ;print(data1)
    
45. ;print(olat)
    
46. ;write_matrix(data1,"lf10.2",opt)
    
47. write_table(opt@fout,"w",[/olon,olat,data1/],"%3i %3i %7.4f")
    
48. end
    

复制代码