%% 公众号:海洋与大气科学
%% 作者mapm
clear;clc;close all;
%% 下面是fnl
filename='fnl_20150816_00_00.nc';area=[0 100;0 50];inter2=10;inter1=10;
%% 读取风速和位势高度;500hpa;
lon = double(ncread(filename,'longitude'));
lat = double(ncread(filename,'latitude'));
GHT_500=double(ncread(filename,'HGT_500mb'));
U_500=double(ncread(filename,'UGRD_500mb'));
V_500=double(ncread(filename,'VGRD_500mb'));
%% 区域;
ln=find(lon>=area(1,1)&lon<=area(1,2));
la=find(lat>=area(2,1)&lat<=area(2,2));
lonx=lon(ln);laty=lat(la);
[x,y]=meshgrid(lonx,laty);
x=x';y=y';
%% 把区域变量找出来;
HGT=GHT_500(ln,la);
U=U_500(ln,la);
V=V_500(ln,la);
%%
% 读取颜色;
% colormore_1-49系列获取,回复”morecolor“
color=colormore_46(2000);
num=46;
fig_h=figure;
set(gcf,'position',[50,50,1000,800],'color','w')
m_proj('miller','lon',[area(1,1) area(1,2)],'lat',[area(2,1) area(2,2)]);
inter=2;
m_windbarb(x(1:inter:end,1:inter:end),y(1:inter:end,1:inter:end),U(1:inter:end,1:inter:end),V(1:inter:end,1:inter:end));
hold on
m_coast('linewidth',1.5,'color',[0.3,0.3,0.3]);
m_grid('ytick',[area(2,1):inter2:area(2,2)],'fontsize',16,'fontweight','bold','linewidth',3);
set(gca,'xtick',[area(1,1):inter1:area(1,2)],'xticklabel',[area(1,1):inter1:area(1,2)],'linewidth',2);
export_fig(fig_h,['colormore_',num2str(num),'大气环流彩图bar'],'-dpng','-r300');
%% ,'color','b','linewidth',2
color=colormore_46(2000);
num=46;
fig_h=figure;
set(gcf,'position',[50,50,1000,800],'color','w')
m_proj('miller','lon',[area(1,1) area(1,2)],'lat',[area(2,1) area(2,2)]);
inter=2;
m_windbarb(x(1:inter:end,1:inter:end),y(1:inter:end,1:inter:end),U(1:inter:end,1:inter:end),V(1:inter:end,1:inter:end),'color','b','linewidth',2);
hold on
m_coast('linewidth',1.5,'color',[0.3,0.3,0.3]);
m_grid('ytick',[area(2,1):inter2:area(2,2)],'fontsize',16,'fontweight','bold','linewidth',3);
set(gca,'xtick',[area(1,1):inter1:area(1,2)],'xticklabel',[area(1,1):inter1:area(1,2)],'linewidth',2);
export_fig(fig_h,['colormore_',num2str(num),'大气环流彩图ba1r'],'-dpng','-r300');
%%