nx=10;ny=5;nz=5;
f=load('fields.dat');
size(f)
% in to out, bottom to top
for ijim=1:length(f(:,1))
 it=0;
 for i=1:nz
  for j=1:ny
   for k=1:nx
    it=it+1;   
    state(k,j,i)=f(ijim,it);
   end
  end
 end 
 toplot(:,:)=state(1,:,:);
 figure(1);contour(toplot');drawnow
end