利用read_grib2读ERA-TIGGE数据

aiai1992

头一回发帖分享~可能比较粗糙，不过应该是干货~

最近在用TIGGE数据，格式是Grib2。开始是一头雾水！在论坛上找了好久，发现有matlab读取GFS的grib2文件的示例，还有二爷的Matlab读取气象海洋数据教程（二爷原创），于是尝试改了一下，写了一个读ERA-TIGGE的小函数。想着论坛里还没有读欧洲格点的示例，就发一下，感谢论坛给我的帮助！


这个函数比较简单，适用性不好，需要跟据下载的数据适当修改，但思路是一样的~


注：附件里一些数据格式的解释说明。

1. function [ x ] =read_tigge(Data_File , Para)
   
2. %-------------load data----------------------------------------------------
   
3. temp = read_grib2(Data_File,{Para});
   
4. %-------------print info --------------------------------------------------
   
5. disp(temp(1).parameter)
   
6. disp(temp(1).unit)
   
7. disp(temp(1).gds)
   
8. 
   
9. if size(temp,1) == 41*8
   
10.     for ti = 1:41
    
11.         for k = 1:8
    
12.             x(:,:,k,ti) = temp(k+(ti-1)*8).data;
    
13.         end
    
14.     end
    
15.     disp('Done!')
    
16. else
    
17.     disp('Error!')
    
18. end
    
19. 
    
20. end
    

复制代码