NCEP CFSR數據讀取

1、

NCEP CFSR再分析數據，時間分辨率是1小時。

一、總體讀取數據狀況

clear all
setup_nctoolbox
tic
%% 讀取數據文件
wind= ncgeodataset('wnd10m.cdas1.201612.grb2');
wind_list = wind.variables;%文件的列表狀況

wind_list包含了

二、u份量

size_of_u = wind.size(wind_list(1));%u份量的數據尺寸,777小時，1個高度，經緯度數據，4D數據
data_u=wind.geovariable(wind_list(1));%取得數據類型爲風速u的數據
u_1=data_u.data(1,1,:,:); %
u_2=squeeze(u_1);

三、v 份量

%% v 份量
size_of_v = wind.size(wind_list(2));%v份量的數據尺寸,777小時，1個高度，經緯度數據，4D數據
data_v=wind.geovariable(wind_list(2));%取得數據類型爲風速v的數據
v_1=data_v.data(1,1,:,:); %
v_2=squeeze(v_1);

四、

高度只有一個值，10米

%% 高度份量
size_of_h= wind.size(wind_list(5));%v份量的數據尺寸,777小時，1個高度，經緯度數據，4D數據
data_h=wind.geovariable(wind_list(5));%取得數據類型爲風速v的數據
v_1=data_h.data(1); %高度10米

五、緯度

 緯度有880 個值。也就是從北緯90度到南緯90度，分紅了880個值，其緯度分辨率爲 0.2045度。

%% 緯度
size_of_lat = wind.size(wind_list(3));%u份量的數據尺寸,777小時，1個高度，經緯度數據，4D數據
data_lat=wind.geovariable(wind_list(1));%取得數據類型爲風速u的數據
lat_1=data_lat.data(1,1,:,:); %
lat_2=squeeze(lat_1);

六、經度

同緯度，從西經到東經，經度的分辨率爲 360/1760 = 0.2045。

%% 經度
size_of_lon = wind.size(wind_list(3));%u份量的數據尺寸,777小時，1個高度，經緯度數據，4D數據
data_lon=wind.geovariable(wind_list(1));%取得數據類型爲風速u的數據
lon_1=data_lon.data(1,1,:,:); %
lon_2=squeeze(lon_1);

七、時間份量

NCEP 一小時採集一次數據，因此 12月份 共 31*24=744 個數據

%% 時間數據
size_of_time = wind.size(wind_list(6));%u份量的數據尺寸,777小時，1個高度，經緯度數據，4D數據
data_time=wind.geovariable(wind_list(6));%取得數據類型爲風速u的數據
time=data_time(1:size_of_time);

八、風速

wind_speed=sqrt(u_2.^2+v_2.^2);

2、帶入時間和座標數據讀取特定地點的風速。