matlab风场图画不出矢量

zzw9701

用matlab画风场的图画不出来，就是没有风矢量，只有地形。。。。，下面是程序，有没有大神能够指点一下：
clear all; close all
%% Import data
ncdisp('D:\biyesheji\ECMWF\Jul_sst.nc');
lon=ncread('D:\biyesheji\ECMWF\Jul_sst.nc','longitude');
lat=ncread('D:\biyesheji\ECMWF\Jul_sst.nc','latitude');
u10=ncread('D:\biyesheji\ECMWF\Jul_sst.nc','u10');
v10=ncread('D:\biyesheji\ECMWF\Jul_sst.nc','v10');
u=u10(:, :, 1);
v=v10(:, :, 1);
[x,y]=meshgrid(lat,lon);
m_proj('Equidistant cylindrical','lat',[5,40],'lon',[140,170]);
h=m_quiver(x,y,u,v,0);
m_coast('color','k');
m_grid('linestyle','none');

二爷名声在外

1.m_quiver最后的参数请不要设为0；2.如果有报错信息，请放出来

zzw9701

二爷名声在外 发表于 2018-3-31 16:39
1.m_quiver最后的参数请不要设为0；2.如果有报错信息，请放出来

我改成1也没有风矢量，matlab没有报错。。。

二爷名声在外

zzw9701 发表于 2018-3-31 16:54
我改成1也没有风矢量，matlab没有报错。。。

上传个数据

zzw9701

二爷名声在外 发表于 2018-3-31 17:01
上传个数据

数据已上传，这个是只有2017年7月每天00时的，原文件太大了。。。。

二爷名声在外

clear all;clc;close all;
filename='_grib2netcdf-atls01-a82bacafb5c306db76464bc7e824bb75-7lCDTL.nc';
lon=ncread(filename,'longitude')-180;
lat=ncread(filename,'latitude');
u10=ncread(filename,'u10');
v10=ncread(filename,'v10');
u=u10(:, :, 1);
v=v10(:, :, 1);
[x,y]=meshgrid(lon,lat);
m_proj('Equidistant cylindrical','lat',[5,40],'lon',[140,170]);
hold on
m_grid
m_coast('color','k');
m_quiver(x,y,u',v');

和数据排列方式有关，多看看过程中的数据，楼主自己参悟吧

zzw9701

恩恩，好的，非常感谢！！！