- 积分
- 3635
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-6-12
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
自己最近在学习如何用matlab绘图,参考网上很多语句,利用气象家园的帖子来记录一下
【今天主要用nc月平均高度场和风场资料简单绘制了7月份多年平均的500hPa环流形势场】
代码如下:
clear;
clc;
%ncdisp(strcat('F:\001\hgt.mon.mean.nc'));查看文件信息
lon=ncread('F:\001\hgt.mon.mean.nc','lon');
lat=ncread('F:\001\hgt.mon.mean.nc','lat');
h=ncread('F:\001\hgt.mon.mean.nc','hgt');
u=ncread('F:\001\uwnd.mon.mean.nc','uwnd');
v=ncread('F:\001\vwnd.mon.mean.nc','vwnd');
x=lon(29:57);y=lat(14:37);
%绘制7月份500hPa环流形势
h_500=squeeze(mean(h(29:57,14:37,6,7:12:475),4));%7月平均位势高度
u_500=squeeze(mean(u(29:57,14:37,6,7:12:475),4));
v_500=squeeze(mean(v(29:57,14:37,6,7:12:475),4));
load coast;
plot(long,lat,'-k','LineWidth',1.8);%画世界陆地界线
hold on;
addpath 'D:\Matlab\地图\m_map1'
%geoshow('D:\Matlab\地图\bou1_4l.shp');调用中国地图
guojie=shaperead('D:\Matlab\地图\chinamap\bou2_4l.shp');
bou2_4lx=[guojie(:).X];%提取经度信息
bou2_4ly=[guojie(:).Y];%提取纬度信息
plot(bou2_4lx,bou2_4ly,'-k','LineWidth',1.2)%绘国界
axis([70 140 15 55]);%设置坐标轴范围
set(gca,'XTick',70:10:140)%设置坐标轴样式
set(gca,'XTickLabel',{'70°E','80°E','90°E','100°E','110°E','120°E','130°E','140°E'})
set(gca,'YTick',15:10:55)
set(gca,'YTickLabel',{'15°N','25°N','35°N','45°N','55°N'});
hold on
quiver(x,y,u_500',v_500');
[C,b]=contour(x,y,h_500');
set(b,'ShowText','on','Textstep',get(b,'LevelStep'));
%LevelStep是显示值为m的倍数的等高线(M可以手动也可以默认),TextStep为n就是显示值为n的倍数的等高线的值
colormap jet
|
-
|