%%%%%%%%%%%%%%%%%%%%%
%  Function slicer  %
%%%%%%%%%%%%%%%%%%%%%
function sli(z1,z2)
global posi T D;

clf;
if exist('posi')==0
 load posi.out;
end

if z1<z2
 z3 = z1;
 z1 = z2;
 z2 = z3;
end;

for i=1:216
 if (posi(i,3)<=z1)&(posi(i,3)>=z2)
  bx = [bx posi(i,1)];
  by = [by posi(i,2)];
  bz = [bz posi(i,3)];
end;
end;

plot(bx,by,'o');
xlabel('X');
ylabel('Y');

ss=input('T = ','s');
if isstr(ss) T = ss; end;
ss=input('Density = ','s');
if isstr(ss) D=ss; end;

z1 = num2str(z1);
z2 = num2str(z2);
title(['A slice at T = ' T 'K Density = ' D ' between ' z2 '<Z<' z1]);

PRINT;

end;




