

clear; clc;

load('fig3.mat');

figure;
set(gca,'color',[1 1 1]);
set(gca,'defaultAxesFontSize',18);
subplot(411);
plot(fig3a.x,fig3a.y,'-r','linewidth',1.0); hold on;
datetick('x','HH:MM'); 
xlabel('Time (UT)'); ylabel('Flux (W/m^2)')
ylim([1e-7,1e-3]);
set(gca,'XLim',[fig3a.xind fig3a.xind+1]);
set(gca,'XTick',[fig3a.xind:1/6:fig3a.xind+1]);
set(gca,'XTickLabel',{'00:00' '04:00' '08:00' '12:00' '16:00' '20:00' '24:00'});
set(gca,'ytick',[1e-7 1e-6 1e-5 1e-4 1e-3])
set(gca,'yscale','log');
set(gca,'xminortick','on');
box on; grid on;
set(gcf,'Position',[300 250 750 300],'Unit','pixels');

subplot(412);
imagesc(fig3b.t,fig3b.wl,log10(fig3b.spectra)); axis xy;
set(gca,'XTick',[0:4:24]);
set(gca,'XTickLabel',{'00:00' '04:00' '08:00' '12:00' '16:00' '20:00' '24:00'});
xlabel('Time (UT)'); ylabel('Wavelength (nm)');
colormap('jet'); caxis([-8 -2]); 

subplot(413);
imagesc(fig3c.t,fig3c.alt,log10(fig3c.ion)); axis xy;
set(gca,'XTick',[1:600:length(fig3c.t)]); set(gca,'XTickLabel',{'17:20' '17:30' '17:40' '17:50' '18:00' '18:10' '18:20' '18:30'})
set(gca,'YLim',[30 100]); set(gca,'YTick',30:20:90); colormap('jet'); caxis([-2 4])
xlabel('Time (UT)'); ylabel('Altitude (km)')
yyaxis right;
semilogy(fig3c.x,fig3c.xflux,'-b','lineWidth',1);
ax=gca; ax.YAxis(2).Color='b';
ylabel('Fux (W/m^2)')

subplot(414)
imagesc(fig3d.t,fig3d.alt(1:151),log10(fig3d.ne)); axis xy;
set(gca,'XTick',[1:600:length(fig3c.t)]); set(gca,'XTickLabel',{'17:20' '17:30' '17:40' '17:50' '18:00' '18:10' '18:20' '18:30'})
set(gca,'YLim',[30 100]); set(gca,'YTick',30:20:90); colormap('jet'); caxis([-2 5])
xlabel('Time (UT)'); ylabel('Altitude (km)')
yyaxis right;
semilogy(fig3d.x,fig3d.xflux,'-b','lineWidth',1);
ax=gca; ax.YAxis(2).Color='b';
ylabel('Flux (W/m^2)')