

pause off;
nx=129;
nz=129;
r=1;
ncol=64;
ns=6;
isx=2;
isz=3;

% number of snapshots

disp(' Loading snapshots ... ')
load -ascii 'snaps'; 

% resampling into arrays

for i=1:ns;
s=reshape(snaps((i-1)*nx*nz+1:i*nx*nz),nz,nx);s=0.5*(s/max(max(s))+r)*ncol;
subplot(isx,isz,i),image(s),colormap('jet'),axis('square'),pause
end
