% Snapshots of vertical fluid velocities for 2 positions 
% over suthern and northern faults (Figure 3 of the manuscript)
set(gcf,'Position',[50,50,1000,800]);
load('w-3d-300sec.mat')
subplot(2,2,1)
plot(squeeze(www(840,840,:)),70.25:0.25:70.25+200*0.25-0.25)
xlabel('m/s')
ylabel('Altitude (km)')
title('42.4°S/173.28°E at T_0 + 300 sec')
load('w-3d-325sec.mat')
subplot(2,2,2)
plot(squeeze(www(840,840,:)),70.25:0.25:70.25+200*0.25-0.25)
xlabel('m/s')
ylabel('Altitude (km)')
title('42.4°S/173.28°E at T_0 + 325 sec')
load('w-3d-355sec.mat')
subplot(2,2,3)
plot(squeeze(www(954,918,:)),70.25:0.25:70.25+200*0.25-0.25)
xlabel('m/s')
ylabel('Altitude (km)')
title('41.95°S/173.75°E at T_0 + 355 sec')
load('w-3d-380sec.mat')
subplot(2,2,4)
plot(squeeze(www(954,918,:)),70.25:0.25:70.25+200*0.25-0.25)
xlabel('m/s')
ylabel('Altitude (km)')
title('41.95°S/173.75°E at T_0 + 380 sec')