- 积分
- 27
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-4-1
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
这是结果图
这是程序:
clear;
file1='wind-1979-2013-850mb.nc';
ncid1=netcdf(file1);
lont=ncid1{'LON47_55'}(:,:,:,:);
latt=ncid1{'LAT45_55'}(:,:,:,:);
lev=ncid1{'LEV15_15'}(:,:,:,:);
v=ncid1{'UWND'}(:,:,:,:);
missing=ncid1{'UWND'}.missing_value;
v(find(v==missing))=NaN;
v0=squeeze(v);
ind=0;
monind=[];
for year=1979:2013;
for mon=1:12
ind=ind+1;
monind(ind)=mon;
end
end
juv=find(monind==1);
v1=v0(juv,:,:);
v1mean=mean(v1,1);
v2=squeeze(v1mean);
figure;
hold on;
%[Lon,Lat]=meshgrid(lon,lat);
m_proj('mercator','lon',[40,60],'lat',[40,60]);
m_contourf(lont,latt,v2,[-4:2:10],'linestyle','none');caxis([-4 10]),colorbar;
m_gshhs_i('LineWidth',1,'Color',[0.5 0.5 0.5]);
%m_gshhs_i('patch',[.6 .6 .6],'edgecolor','b');
m_grid('box','fancy','tickdir','in');
hold on;
%[Lon,Lat]=meshgrid(lon,lat);
%xx=m_quiver(Lon,Lat,v,0,'b');
%m_gshhs_i('LineWidth',1,'Color',[0.5 0.5 0.5]);
%m_grid('box','fancy','YaxisLocation','right','xtick',lona:lonstep:lonb,'ytick',lata:latstep:latb,'tickdir','in','linestyle','none','fontsize',12);
xlabel('\fontsize{14}Longitude');
ylabel('\fontsize{14}Latitude');
问题出在哪啊??
file:///C:/Users/a/AppData/Roaming/Tencent/Users/596079566/QQ/WinTemp/RichOle/GPSPFJW0@4V(6)OR2DR8~)R.png
|
|