1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56
|
figure(indfile)
subplot(2,2,[1 2])
colormap(cmap);
imagesc(thelon, thelat, rot90(GOES_image, 5), [TBmin TBmax]);
hold on
plot(long, lat, 'y-')
xlim([plotLONMIN plotLONMAX]);
ylim([plotLATMIN plotLATMAX]);
title([struc_GOESinf.SATIMG ' - ch' struc_GOESinf.ID ' - ' ...
struc_GOESinf.YYYYMMDD ' - ' struc_GOESinf.HHMN],'Fontsize',15);
scatter(cllon1(indfile,1),cllat1(indfile,1),'r','linewidth',3);
scatter(cllon2(indfile,1),cllat2(indfile,1),'b','linewidth',3);
scatter(cllon3(indfile,1),cllat3(indfile,1),'y','linewidth',3);
scatter(cllon4(indfile,1),cllat4(indfile,1),'c','linewidth',3);
scatter(cllon5(indfile,1),cllat5(indfile,1),'g','linewidth',3);
for ak=1:size(goodok12000,1)
rectangle('Position',[goodok12000(ak,2,1) goodok12000(ak,1,1) 3 2],'EdgeColor','r')
end
hold on
for aq=1:length(vec2)
rectangle('Position',[goodok3000(aq,2,1) goodok3000(aq,1,1) 3 2],'EdgeColor','b')
end
hold on
if length(mat)~=0
for pind=1:size(mat,1)
rectangle('Position',[mat(pind,2) mat(pind,1) 3 2],'edgecolor','g')
end
end
set(gca, 'ydir', 'normal')
subplot(2,2,3)
pcolor(xx,1:20,profvertruc2_lat(:,:)');
axis([-120 -30 0 20])
xlabel('Longitude','fontsize',12)
ylabel('Simulation vertical layer','fontsize',12)
set(gca,'xtick',[-135 -120 -105 -90 -75 -60 -45 -30 -15 0])
set(gca,'xticklabel',{'-135','-120','-105','-90','-75','-60','-45','-30','-15','0'},'fontsize',6)
set(gca,'ytick',plot_lev)
set(gca,'yticklabel',{['',num2str(bb(1:20)),'']},'fontsize',10)
subplot(2,2,4)
pcolor(yy,1:20,profvertruc2_lon(:,:)');
axis([0 65 0 20])
xlabel('Latitude','fontsize',12)
ylabel('Simulation vertical layer','fontsize',12)
set(gca,'xtick',[0 10 20 30 40 50 60 70 80])
set(gca,'xticklabel',{'0','10','20','30','40','50','60','70','80'},'fontsize',6)
set(gca,'ytick',plot_lev)
set(gca,'yticklabel',{['',num2str(bb(1:20)),'']},'fontsize',10) |
Partager