cfsr matlab,NCEP CFSR数据读取

一、

NCEP CFSR再分析数据，时间分辨率是1小时。

1、整体读取数据情况

clear all

setup_nctoolbox

tic

%% 读取数据文件

wind= ncgeodataset('wnd10m.cdas1.201612.grb2');

wind_list = wind.variables;%文件的列表情况

wind_list包含了

2、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);

3、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);

4、

高度只有一个值，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米

5、纬度

纬度有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);

6、经度

同纬度，从西经到东经，经度的分辨率为 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);

7、时间分量

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);

8、风速

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

二、带入时间和坐标数据读取特定地点的风速。