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 vxrr_dat;
load wxrr_dat;
load qzrr_dat;
load wzrr_dat;
load qxri_dat;
load vxri_dat;
load wxri_dat;
load qzri_dat;
load wzri_dat;

subplot(3,2,1)
plot(wn_dat,qxrr_dat,'-',wn_dat,qxri_dat,'--')
xlabel('Wavenumber')
ylabel('Qx,r')
title('Solid Lines: Real Parts')

subplot(3,2,2)
plot(wn_dat,wxrr_dat,'-',wn_dat,wxri_dat,'--')
xlabel('Wavenumber')
ylabel('Wx,r')
title('Dashed Lines: Imaginary Parts')

subplot(3,2,3)
plot(wn_dat,vxrr_dat,'-',wn_dat,vxri_dat,'--')
xlabel('Wavenumber')
ylabel('Vx,r')

subplot(3,2,5)
plot(wn_dat,qzrr_dat,'-',wn_dat,qzri_dat,'--')
xlabel('Wavenumber')
ylabel('Qz,r')

subplot(3,2,6)
plot(wn_dat,wzrr_dat,'-',wn_dat,wzri_dat,'--')
xlabel('Wavenumber')
ylabel('Wz,r')

print integr.ps




