

disp(' ')
disp('  ----------------------- ')
disp(' ')
disp('  Basic snaphot plotting ')
disp(' ')
disp('  ----------------------- ')
disp(' ')
disp(' ')

nx=100;
nz=100;

load rot_xs
load rot_zs
load rot_rys

snapx=rot_xs;
snapz=rot_zs;
roty=rot_rys;

ns=floor(length(snapx)/(nx*nz));

disp(' ')
disp(sprintf(' %g snapshots will be displayed ... ',ns))
disp(' ')

j0=1;
jinc=1;

isnap=0;
for j=j0:jinc:ns;
isnap=isnap+1;
sx=reshape(snapx((j-1)*nz*nx+1:j*nz*nx),nz,nx);
sz=reshape(snapz((j-1)*nz*nx+1:j*nz*nx),nz,nx);
ry=reshape(roty((j-1)*nz*nx+1:j*nz*nx),nz,nx);

subplot(131),imagesc(sx),axis('equal'),axis('off'),
subplot(132),imagesc(sz),axis('equal'),axis('off'),
subplot(133),imagesc(ry),axis('equal'),axis('off'),
colormap gray
disp(' Press any key to continue ')
pause
end

