0;


set(0,'DefaultAxesFontName','Times')
set(0,'DefaultAxesFontSize',10)
widthCM=13.5; % width of plots in cm.  JFM width is 135mm.


fig = figure('units','centimeters','position',...
    [5 -5 widthCM 0.7*widthCM]);

tiledlayout(2,2)
nexttile

load('Zvalues04aspect.mat','omega_grid','Z_grid')
plot(omega_grid*4,log(Z_grid),'Color','k')
xlabel('$$\omega R$$')
ylabel('$$\log|Z|$$')
subtitle('(a)')

nexttile

load('Zvalues08aspect.mat','omega_grid','Z_grid')
plot(omega_grid*8,log(Z_grid),'Color','k')
xlabel('$$\omega R$$')
ylabel('$$\log|Z|$$')
subtitle('(b)')

nexttile

load('Zvalues16aspect.mat','omega_grid','Z_grid')
plot(omega_grid*16,log(Z_grid),'Color','k')
xlabel('$$\omega R$$')
ylabel('$$\log|Z|$$')
subtitle('(c)')


nexttile

load('Zvalues16aspect_bent.mat','omega_grid','Z_grid')
plot(omega_grid*16,real(log(Z_grid)),'Color','k')
xlabel('$$\omega R$$')
ylabel('$$\log|Z|$$')
subtitle('(d)')

exportgraphics(gcf,'inletimpedancetiles.pdf','ContentType','vector')

