% This script loads and plots a meridional slice (with altitude) of absolute temperature
% perturbations at specified epoch. Data is structured as [time,latitude,altitude], where
% t1=05:46:29 UT, dt = 5 sec and altitude grid size is 1 km.
% To avoid problems with loading big files, the output is split into 3
% files

close all;clear all; clc;

load('temps1.mat')
load('map.mat')
slice =1;
t = 1100;
pcolor(linspace(cell2mat(map{1,slice}(1)),cell2mat(map{1098,slice}(1)),1098),1:1:500,(squeeze(temps(t,1:1:1098,:))'));
shading interp
axis xy
ylabel('Altitude')
xlabel('Geographic longitude (°)')
colorbar