%
% Generates jpeg images from output_movie file
%
fil=input('Enter input file:')
fid = fopen(fil);
t=[];
a=fscanf(fid, '%g', [1,1]);
b=fscanf(fid, '%g', [1,1]);
c=fscanf(fid, '%g', [1,1]);
d=fscanf(fid, '%g', [1,1]);
m=fscanf(fid, '%g', [1,1]);
n=fscanf(fid, '%g', [1,1]);
nstop=fscanf(fid, '%g', [1,1]);
x=fscanf(fid, '%g', [1,m]);
y=fscanf(fid, '%g', [1,n]);

f=[];
for j = 1:1:nstop
j
  t=fscanf(fid, '%g', [1,m*n]);
l = 0;
  for i = 1:1:m
  for k = 1:1:n
  l = l+1;
       f(i,k) = t(1,l);
  end
  end
  surf(x,y,f);
  axis([a b c d -1 1]);
%axis([0 32 0 15 -0.1e-6 0.1e-6]);
  ofil = sprintf('%s.%04d.jpg',fil,j);
  print('-djpeg','-r70',ofil);
  sprintf('%s.%04d.jpg written.',fil,j)
  pause(0.01);
end
fclose(fid);