登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册 
 
 
 
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,否则就会出现如上情况 
 
 
 
 
 
 
 
 
 
 
 
 
 |