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 qxr_dat;
load vxr_dat;
load wxr_dat;
load qzr_dat;
load wzr_dat;
load qxi_dat;
load vxi_dat;
load wxi_dat;
load qzi_dat;
load wzi_dat;

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

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

subplot(3,2,3)
plot(wn_dat,vxr_dat,'-',wn_dat,vxi_dat,'--')
xlabel('Wavenumber')
ylabel('Vx')

subplot(3,2,5)
plot(wn_dat,qzr_dat,'-',wn_dat,qzi_dat,'--')
xlabel('Wavenumber')
ylabel('Qz')

subplot(3,2,6)
plot(wn_dat,wzr_dat,'-',wn_dat,wzi_dat,'--')
xlabel('Wavenumber')
ylabel('Wz')

print integ.ps




