- 积分
- 27
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-4-1
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
clear;
file1='wind-1979-2013-850mb.nc';
ncid1=netcdf(file1);
lon=ncid1{'LON47_55'}(:,:,:,:);
lat=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;
m_proj('mercator','lon',[40,60],'lat',[40,60]);
m_contourf(lon,lat,v2,30,'linestyle','none');
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');
|
|