clear
%close

fig=figure('Units','centimeters','Position',[0.0 0.0 21.0 26.0]...
,'PaperUnits','centimeters','PaperPosition',[0.0 0.0 21.0 26.0]);

load wn_dat;
size(wn_dat)
load qxrr_dat;
load qxri_dat;
load qxzr_dat;
load qxzi_dat;

subplot(2,2,1)
plot(wn_dat,qxrr_dat,'-')
xlabel('Wavenumber')
ylabel('Real(Qx,r)')
title('Real(Qx,r)')

subplot(2,2,2)
plot(wn_dat,qxri_dat,'-')
xlabel('Wavenumber')
ylabel('Imag(Qx,r)')
title('Imag(Qx,r)')

subplot(2,2,3)
plot(wn_dat,qxzr_dat,'-')
xlabel('Wavenumber')
ylabel('Real(Qx,z)')
title('Real(Qx,z)')

subplot(2,2,4)
plot(wn_dat,qxzi_dat,'-')
xlabel('Wavenumber')
ylabel('Imag(Qx,z)')
title('Imag(Qx,z)')

print integ1.ps




