留言板
- xueqixiang
- 谢谢 ，y长度 和MHGT纬度都没对上。已解决。
ljjtc: [C,H]=m_contour(x(17:46,8:85),y(19:45,8:85),MHGT(:),50,'ShowText','on',...
这里应该是y的长度错了
- xueqixiang
- 您好,能帮我看下2016年200高度场距平的程序吗?图就是出不来
hgt1 = ncread('D:\matlab\MATLAB\wavelet\level\hgt.mon.mean.nc','hgt');
lon = ncread('D:\matlab\MATLAB\wavelet\level\hgt.mon.mean.nc','lon');
lat = ncread('D:\matlab\MATLAB\wavelet\level\hgt.mon.mean.nc','lat');
%hgt2 = hgt(8:85,17:46,10,97:456);
ran = 360;
[x,y] = meshgrid(lon,lat);
hgt_mn=zeros(ran,2340);
for i=1:ran
squ=squeeze(hgt1(8:85,17:46,10,i));
hgt_mn(i,:)=reshape(squ,1,2340);
end
hgt_nan=isnan(hgt_mn);
k=0;
for j=1:2340
if sum(hgt_nan(:,j)==0)
k=k+1;
hgt_(:,k)=hgt_mn(:,j);
end
end
HGT=zeros(size(hgt_));
% 求每月的距平
for i=1:12
%
HGT(i:12:ran-12+i,:)=hgt_(i:12:end,:)-repmat(mean(hgt_(i:12:end,:),1),size(hgt_(i:12:end,:),1),1);
end
MHGT = mean(HGT(349:360,1:2340));
m_proj('Equidistant cylindrical','long',[9.5 212],'lat',[-24.5 52]);
m_coast('patch',[0.98 0.98 0.98],'edgecolor','k');
m_grid('box','on','xtick',10,'tickdir','in','yaxislocation','left','fontsize',10);
xlabel('经度','fontsize',12);
ylabel('纬度','fontsize',12);
title('200hpa位势高度距平','fontsize',15);
hold on
[C,H]=m_contour(x(17:46,8:85),y(19:45,8:85),MHGT(:,:),50,'ShowText','on',...
'Color','r','LineWidth',0.8);
clabel(C,H,'labelspacing',1000);
运行后的错误提示:
Error in
heightanomaly3
(line 38)
[C,H]=m_contour(x(17:46,8:85),y(19:45,8:85),MHGT(:),50,'ShowText','on',...