%% This script and data support Figure 10 of the manuscript %% Figure 10, v1 close all clear all clc load('v1.mat') FS=10; h=imagesc(xp,zp,parmp); set(h,'alphadata',~isnan(parmp)); set(gca,'FontSize',FS); axis xy; colormap(parula(256)); c=colorbar; parmlbl = 'Field-aligned drift velocity (m/s)'; xlabel(c,parmlbl); xlabel('magnetic latitude (deg.)'); ylabel('altitude (km)'); xlim([14 22]) ylim([130 600]) caxis([-145 145]) UThrs=floor(t/3600); UTmin=floor((t/3600-UThrs)*60); UTsec=floor((t/3600-UThrs-UTmin/60)*3600); UThrsstr=num2str(UThrs); UTminstr=num2str(UTmin); if (numel(UTminstr)==1) UTminstr=['0',UTminstr]; end UTsecstr=num2str(UTsec); if (numel(UTsecstr)==1) UTsecstr=['0',UTsecstr]; end timestr=[UThrsstr,':',UTminstr,':',UTsecstr]; strval=sprintf('%s \n %s',[num2str(ymd(2)),'/',num2str(ymd(3)),'/',num2str(ymd(1))], ... [timestr,' UT']); title(strval); %% Figure 10, ne close all clear all clc load('ne1.mat') FS=10; h=imagesc(xp,zp,parmp); set(h,'alphadata',~isnan(parmp)); set(gca,'FontSize',FS); axis xy; colormap(parula(256)); c=colorbar; parmlbl = 'Electron density perturbations (%)'; xlabel(c,parmlbl); xlabel('magnetic latitude (deg.)'); ylabel('altitude (km)'); xlim([14 22]) ylim([130 600]) caxis([-30 30]) UThrs=floor(t/3600); UTmin=floor((t/3600-UThrs)*60); UTsec=floor((t/3600-UThrs-UTmin/60)*3600); UThrsstr=num2str(UThrs); UTminstr=num2str(UTmin); if (numel(UTminstr)==1) UTminstr=['0',UTminstr]; end UTsecstr=num2str(UTsec); if (numel(UTsecstr)==1) UTsecstr=['0',UTsecstr]; end timestr=[UThrsstr,':',UTminstr,':',UTsecstr]; strval=sprintf('%s \n %s',[num2str(ymd(2)),'/',num2str(ymd(3)),'/',num2str(ymd(1))], ... [timestr,' UT']); title(strval); %% Figure 10, temperature close all clear all clc load('t1.mat') FS=10; h=imagesc(xp,zp,parmp); set(h,'alphadata',~isnan(parmp)); set(gca,'FontSize',FS); axis xy; colormap(parula(256)); c=colorbar; parmlbl = 'Ion temperature perturbations (%)'; xlabel(c,parmlbl); xlabel('magnetic latitude (deg.)'); ylabel('altitude (km)'); xlim([14 22]) ylim([130 600]) caxis([-30 30]) UThrs=floor(t/3600); UTmin=floor((t/3600-UThrs)*60); UTsec=floor((t/3600-UThrs-UTmin/60)*3600); UThrsstr=num2str(UThrs); UTminstr=num2str(UTmin); if (numel(UTminstr)==1) UTminstr=['0',UTminstr]; end UTsecstr=num2str(UTsec); if (numel(UTsecstr)==1) UTsecstr=['0',UTsecstr]; end timestr=[UThrsstr,':',UTminstr,':',UTsecstr]; strval=sprintf('%s \n %s',[num2str(ymd(2)),'/',num2str(ymd(3)),'/',num2str(ymd(1))], ... [timestr,' UT']); title(strval); %% Figure 10, ,v1 close all clear all clc load('v2.mat') FS=10; h=imagesc(xp,zp,parmp); set(h,'alphadata',~isnan(parmp)); set(gca,'FontSize',FS); axis xy; colormap(parula(256)); c=colorbar; parmlbl = 'Field-aligned drift velocity (m/s)'; xlabel(c,parmlbl); xlabel('magnetic latitude (deg.)'); ylabel('altitude (km)'); xlim([14 22]) ylim([130 600]) caxis([-200 200]) UThrs=floor(t/3600); UTmin=floor((t/3600-UThrs)*60); UTsec=floor((t/3600-UThrs-UTmin/60)*3600); UThrsstr=num2str(UThrs); UTminstr=num2str(UTmin); if (numel(UTminstr)==1) UTminstr=['0',UTminstr]; end UTsecstr=num2str(UTsec); if (numel(UTsecstr)==1) UTsecstr=['0',UTsecstr]; end timestr=[UThrsstr,':',UTminstr,':',UTsecstr]; strval=sprintf('%s \n %s',[num2str(ymd(2)),'/',num2str(ymd(3)),'/',num2str(ymd(1))], ... [timestr,' UT']); title(strval); %% Figure 10, ne close all clear all clc load('ne2.mat') FS=10; h=imagesc(xp,zp,parmp); set(h,'alphadata',~isnan(parmp)); set(gca,'FontSize',FS); axis xy; colormap(parula(256)); c=colorbar; parmlbl = 'Electron density perturbations (%)'; xlabel(c,parmlbl); xlabel('magnetic latitude (deg.)'); ylabel('altitude (km)'); xlim([14 22]) ylim([130 600]) caxis([-30 30]) UThrs=floor(t/3600); UTmin=floor((t/3600-UThrs)*60); UTsec=floor((t/3600-UThrs-UTmin/60)*3600); UThrsstr=num2str(UThrs); UTminstr=num2str(UTmin); if (numel(UTminstr)==1) UTminstr=['0',UTminstr]; end UTsecstr=num2str(UTsec); if (numel(UTsecstr)==1) UTsecstr=['0',UTsecstr]; end timestr=[UThrsstr,':',UTminstr,':',UTsecstr]; strval=sprintf('%s \n %s',[num2str(ymd(2)),'/',num2str(ymd(3)),'/',num2str(ymd(1))], ... [timestr,' UT']); title(strval); %% Figure 10, temperature close all clear all clc load('t2.mat') FS=10; h=imagesc(xp,zp,parmp); set(h,'alphadata',~isnan(parmp)); set(gca,'FontSize',FS); axis xy; colormap(parula(256)); c=colorbar; parmlbl = 'Ion temperature perturbations (%)'; xlabel(c,parmlbl); xlabel('magnetic latitude (deg.)'); ylabel('altitude (km)'); xlim([14 22]) ylim([130 600]) caxis([-30 30]) UThrs=floor(t/3600); UTmin=floor((t/3600-UThrs)*60); UTsec=floor((t/3600-UThrs-UTmin/60)*3600); UThrsstr=num2str(UThrs); UTminstr=num2str(UTmin); if (numel(UTminstr)==1) UTminstr=['0',UTminstr]; end UTsecstr=num2str(UTsec); if (numel(UTsecstr)==1) UTsecstr=['0',UTsecstr]; end timestr=[UThrsstr,':',UTminstr,':',UTsecstr]; strval=sprintf('%s \n %s',[num2str(ymd(2)),'/',num2str(ymd(3)),'/',num2str(ymd(1))], ... [timestr,' UT']); title(strval); %% Figure 10, v1 close all clear all clc load('v3.mat') FS=10; h=imagesc(xp,zp,parmp); set(h,'alphadata',~isnan(parmp)); set(gca,'FontSize',FS); axis xy; colormap(parula(256)); c=colorbar; parmlbl = 'Field-aligned drift velocity (m/s)'; xlabel(c,parmlbl); xlabel('magnetic latitude (deg.)'); ylabel('altitude (km)'); xlim([14 22]) ylim([130 600]) caxis([-50 50]) UThrs=floor(t/3600); UTmin=floor((t/3600-UThrs)*60); UTsec=floor((t/3600-UThrs-UTmin/60)*3600); UThrsstr=num2str(UThrs); UTminstr=num2str(UTmin); if (numel(UTminstr)==1) UTminstr=['0',UTminstr]; end UTsecstr=num2str(UTsec); if (numel(UTsecstr)==1) UTsecstr=['0',UTsecstr]; end timestr=[UThrsstr,':',UTminstr,':',UTsecstr]; strval=sprintf('%s \n %s',[num2str(ymd(2)),'/',num2str(ymd(3)),'/',num2str(ymd(1))], ... [timestr,' UT']); title(strval); %% Figure 10, ne close all clear all clc load('ne3.mat') FS=10; h=imagesc(xp,zp,parmp); set(h,'alphadata',~isnan(parmp)); set(gca,'FontSize',FS); axis xy; colormap(parula(256)); c=colorbar; parmlbl = 'Electron density perturbations (%)'; xlabel(c,parmlbl); xlabel('magnetic latitude (deg.)'); ylabel('altitude (km)'); xlim([14 22]) ylim([130 600]) caxis([-6 6]) UThrs=floor(t/3600); UTmin=floor((t/3600-UThrs)*60); UTsec=floor((t/3600-UThrs-UTmin/60)*3600); UThrsstr=num2str(UThrs); UTminstr=num2str(UTmin); if (numel(UTminstr)==1) UTminstr=['0',UTminstr]; end UTsecstr=num2str(UTsec); if (numel(UTsecstr)==1) UTsecstr=['0',UTsecstr]; end timestr=[UThrsstr,':',UTminstr,':',UTsecstr]; strval=sprintf('%s \n %s',[num2str(ymd(2)),'/',num2str(ymd(3)),'/',num2str(ymd(1))], ... [timestr,' UT']); title(strval); %% Figure 10, temperature close all clear all clc load('t3.mat') FS=10; h=imagesc(xp,zp,parmp); set(h,'alphadata',~isnan(parmp)); set(gca,'FontSize',FS); axis xy; colormap(parula(256)); c=colorbar; parmlbl = 'Ion temperature perturbations (%)'; xlabel(c,parmlbl); xlabel('magnetic latitude (deg.)'); ylabel('altitude (km)'); xlim([14 22]) ylim([130 600]) caxis([-5 5]) UThrs=floor(t/3600); UTmin=floor((t/3600-UThrs)*60); UTsec=floor((t/3600-UThrs-UTmin/60)*3600); UThrsstr=num2str(UThrs); UTminstr=num2str(UTmin); if (numel(UTminstr)==1) UTminstr=['0',UTminstr]; end UTsecstr=num2str(UTsec); if (numel(UTsecstr)==1) UTsecstr=['0',UTsecstr]; end timestr=[UThrsstr,':',UTminstr,':',UTsecstr]; strval=sprintf('%s \n %s',[num2str(ymd(2)),'/',num2str(ymd(3)),'/',num2str(ymd(1))], ... [timestr,' UT']); title(strval);