QUIKSCAT/NCEP数据处理（已解决）

言深深

夏天爱裤衩 发表于 2012-7-6 16:04

                               
登录/注册后可看大图

着实不会写程序啊，只是想把数据中任意经纬度范围、时间范围的数据读出来就可以，刚刚把end=99删掉，编译 ...

这里的意思应该是read到标记99行的地方结束，但是标记99的地方居然是一个close（11）表示关闭文件，很奇怪，一般都是end，enddo，continue啥的，所以不太明白程序编写的人到底想干的啥啊

夏天爱裤衩

上传一个ncl处理QUIKSCAT/NCEP数据的脚本，只要改动do i = 0, 10
      ip1 = 20
      ip2 = 24
      jp  = 19这几个位置就可以

1. ; readqnblendwind.ncl
   
2. ; NCL program to read binary QSCAT/NCEP blended wind data files.
   
3. ; This routine can be adapted to read either the zonal and
   
4. ; meridional wind data ("uv") or the windstress curl data ("curl").
   
5. ;
   
6. ; Contributed by: Tianyi Fan, University of Colorado, Boulder; 2007-07-12
   
7. ;
   
8. 
   
9. begin
   
10. nlat = 353
    
11. nlon = 720
    
12. 
    
13. readuv = True   ; to read either uv or curlfiles
    
14. ;  readuv = False
    
15. 
    
16. ; ===============================
    
17. ;  open file
    
18. ; ===============================
    
19. if (readuv) then
    
20.    fils = "./uv.200001.bln"
    
21. else
    
22.    fils = "./curl.200001.bln"
    
23. end if   
    
24. print("intput" + fils)
    
25. 
    
26. 
    
27.   numRec = 30*4  ; number of record to read in, mamybe 31*4 or
    
28. 
    
29. day = new( numRec, "float")
    
30. u = new((/ numRec, nlat, nlon/), "float")
    
31. v = new((/ numRec, nlat, nlon/), "float")
    
32. curl = new((/ numRec, nlat, nlon/), "float")
    
33. 
    
34. if (readuv) then
    
35.    recl = 7+2*nlon*nlat      ; recl length
    
36. else
    
37.    recl = 5+nlon*nlat
    
38. end if  
    
39. dum_all = new(recl, "float")
    
40. 
    
41. print("Assigning coordinate varibal information")
    
42. 
    
43. u!0 = "time"
    
44. u&time = ispan(1, numRec, 1)
    
45. u!1 = "lat"
    
46. u&lat = fspan(-88., 88., nlat)
    
47. u&lat@units = "degrees_north"
    
48. u!2 = "lon"
    
49. u&lon = fspan(0.5,360.,nlon)
    
50. u&lon@units = "degrees_east"
    
51. 
    
52. v!0 = "time"
    
53. v&time = ispan(1, numRec, 1)
    
54. v!1 = "lat"
    
55. v&lat = fspan(-88., 88., nlat)
    
56. v&lat@units = "degrees_north"
    
57. v!2 = "lon"
    
58. v&lon = fspan(0.5,360.,nlon)
    
59. v&lon@units = "degrees_east"
    
60. 
    
61. 
    
62. curl!0 = "time"
    
63. curl&time = ispan(1, numRec, 1)
    
64. curl!1 = "lat"
    
65. curl&lat = fspan(-88, 88, nlat)
    
66. curl&lat@units = "degrees_north"
    
67. curl!2 = "lon"
    
68. curl&lon = fspan(0.5,360.,nlon)
    
69. curl&lon@units = "degrees_east"
    
70. 
    
71. ; ===============================
    
72. ;  read data
    
73. ; ===============================
    
74.   if(readuv) then
    
75.     do i = 0, numRec -1   
    
76.       dum_all    = fbindirread(fils, i, recl, "float")      
    
77. ;      idum1      = dum_all(0)
    
78.       day(i)     = dum_all(1)
    
79. ;      idum2      = dum_all(2)
    
80. ;      idum3      = dum_all(3)
    
81.       u(i, :,:)  = onedtond(dum_all(4:4+nlon*nlat-1),(/nlat, nlon/) )
    
82. ;      idum4      = dum_all(4+nlon*nlat)
    
83. ;      idum5      = dum_all(5+nlon*nlat)
    
84.       v(i, :,:)  = onedtond(dum_all(6+nlon*nlat:6+nlon*nlat*2-1),(/nlat, nlon/) )
    
85. ;      idum6      = dum_all(6+nlon*nlat*2)
    
86.     end do
    
87.   else
    
88.     do i = 0, numRec -1
    
89.       dum_all    = fbindirread(fils, i, recl, "float")  
    
90.       idum1      = dum_all(0)
    
91.       day(i)     = dum_all(1)
    
92.       idum2      = dum_all(2)
    
93.       idum3      = dum_all(3)
    
94.       curl(i,:,:)= onedtond(dum_all(4:4+nlon*nlat-1),(/nlat, nlon/) )
    
95.       idum4      = dum_all(4+nlon*nlat)        
    
96.     end do
    
97.   end if
    
98.   
    
99. 
    
100. ;============================
     
101. ;  print
     
102. ; ===========================
     
103.   if (readuv) then
     
104.     do i = 0, 10
     
105.       ip1 = 20
     
106.       ip2 = 24
     
107.       jp  = 19
     
108.       print("day = " + day(i) + "  print data at i = " + (ip1+1) + "  to  " + (ip2+1) + ", j = " + (jp+1)  )
     
109.       print(" u = " + u(i, jp, ip1:ip2))
     
110.       print(" v = " + v(i, jp, ip1:ip2))
     
111.       
     
112.       jp  = 299
     
113.       print("day = " + day(i) + "  print data at i = " + (ip1+1) + "  to  " + (ip2+1) + ", j = " + (jp+1)  )
     
114.       print(" u = " + u(i, jp, ip1:ip2))
     
115.       print(" v = " + v(i, jp, ip1:ip2))
     
116.     end do  
     
117.   else
     
118.     do i = 0, 5
     
119.       ip1 = 20
     
120.       ip2 = 24
     
121.       jp  = 19
     
122.       print("day = " + day(i) + "  print data at i = " + (ip1+1) + "  to  " + (ip2+1) + ", j = " + (jp+1)  )
     
123.       print(" c = " + curl(i, jp, ip1:ip2))
     
124.       
     
125.       jp  = 299
     
126.       print("day = " + day(i) + "  print data at i = " + (ip1+1) + "  to  " + (ip2+1) + ", j = " + (jp+1)  )
     
127.       print(" c = " + curl(i, jp, ip1:ip2))
     
128.     end do   
     
129.   end if
     
130. 
     
131.   
     
132. end
     

复制代码

mofangbao

夏天爱裤衩 发表于 2012-7-6 19:28

                               
登录/注册后可看大图

上传一个ncl处理QUIKSCAT/NCEP数据的脚本，只要改动do i = 0, 10
      ip1 = 20
      ip2 = 24

哈哈  又从哪儿要了一个啊

夏天爱裤衩

mofangbao 发表于 2012-7-6 19:52

                               
登录/注册后可看大图

哈哈  又从哪儿要了一个啊

哦哈哈，网上找啊，我还有个matlab的呢，不过到处都是错误，改了半天没搞出结果来，如果有人需要，我也可以上传

小肥虾619

请问楼主知不知道QuikSCAT资料的原始格式啊？以及到底是6小时一次还是12小时一次具体的时间间隔是啥样的？谢啦谢啦

夏天爱裤衩

小肥虾619 发表于 2012-10-30 09:01

                               
登录/注册后可看大图

请问楼主知不知道QuikSCAT资料的原始格式啊？以及到底是6小时一次还是12小时一次具体的时间间隔是啥样的？谢 ...

QuikSCAT和QuikSCAT/NCEP还是有区别的，你仔细找找看，后者是QuikSCAT卫星数据融合了NCEP再分析数据而成的混合风场，为6h分辨率，前者没用过，应该是按轨道走的吧

小肥虾619

夏天爱裤衩 发表于 2012-10-30 13:50

                               
登录/注册后可看大图

QuikSCAT和QuikSCAT/NCEP还是有区别的，你仔细找找看，后者是QuikSCAT卫星数据融合了NCEP再分析数据而成的 ...

多谢楼主了哈

jeffggg

你好 这个QUickscat 数据有没有NC格式的 下载地址是哪个啊？

geyx

正着急怎么处理呢   正需要帮助啊

老A不加V

夏天爱裤衩 发表于 2012-10-30 13:50
QuikSCAT和QuikSCAT/NCEP还是有区别的，你仔细找找看，后者是QuikSCAT卫星数据融合了NCEP再分析数据而成 ...

你好，想问下，你最后是怎么处理的这个数据，我现在也正在使用这个数据，不知如何处理？能否处理成Grads可识别的格式？