登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 灰色节能君 于 2020-3-20 21:46 编辑
数据是ECMWF下载的u,v风场nc数据,然后用m_quiver画不出风场,只在90°S上显示了一行箭头,数据里没有缺测值,代码如下,数据我也传上去了,恳请大神帮我看一下哪里出问题了clear all;clc;close all;
filename='_grib2netcdf-atls04-a562cefde8a29a7288fa0b8b7f9413f7-4DYffR.nc';
lon=ncread(filename,'longitude')+180;
lat=ncread(filename,'latitude');
u0=ncread(filename,'u');
v0=ncread(filename,'v');
u=u0(:, :, 1);
v=v0(:, :, 1);
[x,y]=meshgrid(lon,lat);
m_proj('Equidistant cylindrical','lat',[-90,90],'lon',[0,360]);
m_grid('linestyle','none','xtick',6,'fontsize',14,'tickdir','out','ytick',5);
m_coast('color','k');
hold on
m_quiver(x,y,u',v',10);
%%%%已解决
风场的数据不能包括90°N和90°S,否则就会出现如上情况
|