- 积分
- 1473
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-28
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
用m_contour绘制了高度场,等值线的标注太密了,能不能隔几条线标注一下?谢谢各位老师;
程序如下:
lat= ncread('vwnd.mon.mean.nc','lat');
lon= ncread('vwnd.mon.mean.nc','lon');
v= ncread('vwnd.mon.mean.nc','vwnd');
u1= ncread('uwnd.mon.mean.nc','uwnd');
H = ncread('hgt.mon.mean.nc','hgt');
h = H(:,:,10,450);
[x,y] = meshgrid(lon,lat);
v1 = v(:,:,10,450);
v1(v1==-9999) = nan;
u2 = u1(:,:,10,450);
u2(u2==-9999) = nan;
u3 = u2.*cosd(y');
d = 2 ;dd = 2;
uv = sqrt(u3.^2+v1.^2);
m_proj('Equidistant cylindrical','long',[18 212],'lat',[-22 47]);
m_coast('patch',[0.98 0.98 0.98],'edgecolor','k');
m_grid('box','on','xtick',10,'tickdir','in','yaxislocation','left','fontsize',10);
xlabel('longitude');
ylabel('latitudel');
title('风场和位势高度场','fontsize',15);
%m_ruler([0 0.5],1);
hold on;
m_quiver(x(19:45,9:85),y(19:45,9:85),u3(9:85,19:45)'./dd,v1(9:85,19:45)'./dd,1);
% m_pcolor(x,y,double(uv'));
hold on;
m_contour(x(19:45,9:85),y(19:45,9:85),h(9:85,19:45)',25,'ShowText','on',...
'Color','r','LineWidth',0.8);
图片如下:
|
-
|