NCEP上下载u、v风速nc文件后如何通过matlab编程实现风向的显示

maxwellll

ncep上下载u、v风速nc文件后，该如何通过matlab编程实现风向在地图上（10°N-40°N,110°E-150°E）的显示呢。在网上找了代码，但是调试不成功，毕业设计要用到，但我不是本专业的，matlab不是很熟悉，希望有大大帮忙答疑解惑，感谢。

maxwellll

还有想问一下 在matlab中如何读取nc文件并提取相关数据呢

ゼヅ欢チニ

读取nc文件，并将津京冀地面温度取出来放到tem里
该nc数据为0.125°*0.125°，月均值数据
clear all;
clc
track='f:\DATA\ECMWF\ECMWF-Levels\';
file=dir([track,'*.nc']);
lat1=90/0.125-(35/0.125);  %找到对应经纬度的位置
lat0=90/0.125-(45/0.125);  
lon0=110/0.125;
lon1=120/0.125;
for i=1:length(file)
    str=file(i).name;      
    str1=strcat(track,str);
    ncdisp(strcat(track,str)); %可以显示nc中有什么变量及其维度等，
    lat = double(ncread(str1,'latitude'));
    lon = double(ncread(str1,'longitude'));
% time = ncread(str1,'time');
    t = ncread(str1,'t')-272.15;          %K-272.15=℃
    tem=t([lon0:lon1],[lat0:lat1],16,1);     %第16层高度为1000hpa，1代表1月            
end
laty=lat(lat0:lat1);
lonx=lon(lon0:lon1);
pcolor(lonx,laty,tem)
shading flat; %  加上这句即可显示
[x,y]=meshgrid(lonx,laty);
phandle=pcolor(x,y,tem');
colorbar

ENSO爱好者

ゼヅ欢チニ 发表于 2019-6-26 08:51
读取nc文件，并将津京冀地面温度取出来放到tem里
该nc数据为0.125°*0.125°，月均值数据
clear all;

你用的是ECWMF数据是吧？

ゼヅ欢チニ

ENSO爱好者 发表于 2019-6-27 21:01
你用的是ECWMF数据是吧？

恩恩  对 不是NCEP数据    忘标题了

江湖小鱼儿

请问楼主您问题解决了吗？我也遇到类似的问题，能否分享一下解决的办法呢？感激不尽