% "short" plot of an "inner" packet...

function [] = plot_inner_dataset(num)
global idata data

odata=data;

cd /var/tmp

eval (sprintf('load packet%d_r', num));
eval (sprintf('load packet%d_i', num));
data = fftshift( eval( sprintf('(packet%d_r + 1i*packet%d_i)', num, num))' );
% Q=fftshift(fft2(data));
eval (sprintf('clear packet%d_r packet%d_i', num, num));

sz = size(data);
ii = (sz(1)/2);
jj = (sz(2)/2);

% figure(1);
mesh([-ii:(ii-1)],[-jj:(jj-1)],abs(data),angle(data));
set(gca,'XLim',[-ii ii],'YLim',[-jj jj],  ...
        'XTick',[-ii -ii/2 0 ii/2 ii],'YTick',[-jj -jj/2 0 jj/2 jj]);
%set(gca,'Visible','off');
title(sprintf('wave function after %d segments (c-sp)', num));

% figure(2);
% mesh(abs(Q),angle(Q));
% title(sprintf('momentum-space packet after %d segments', num));
