%
%
freq = 100.0
range = [100 500 1000 1500 5000 10000 15000 20000];
figctr = 1;

for rn = range,
 eval(['load /tmp/k.', num2str(freq), '.0.', num2str(rn), '.out -ascii']);
 eval(['load /tmp/r.', num2str(freq), '.0.', num2str(rn), '.out -ascii']);
 ki = k(:,2) + i .* k(:,3);
 ri = r(:,2) + i .* r(:,3);

 figure(figctr);
 plot(k(:,1), abs(ki), r(:,1), abs(ri));
 title (['kraken/ram comparison, ', num2str(rn), ' meters']);
 legend ('kraken', 'ram');
 figctr = figctr + 1;
end