%% 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);