clear 
clear path 

global Nx Ny Nz
global lat long dz dm mdep
global delt_su su_its t_su delt
global descriptor this_path
global f deltaf Q beta r_expt r_heat H
global time rots it
global g Cp rho_bar alpha
global u v t w
global iterations


param_file_name =   ...
   input(' Please enter the name of the m-file with the parameters for this run : ','s') ;
feval(param_file_name) 

sizeit=size(iterations);
  for i=1:round(sizeit(1)/16)
figurenumber=i;
figurenumber
      figure(i)
title(descriptor)
for j=1:16;
figurenumber
subplotnumber=j;
subplotnumber
((i-1)*16)+j
   it = iterations(((figurenumber-1)*16)+j,1:10);
time   = (it-su_its)*delt ;
rots = time*f/(4*pi);
   path   = this_path 
   cmdstr=['cd ' path ];
   eval(cmdstr);
   path=pwd
   % T (oC)
   tfilename=(['T.' it ]);
   t=rdmeta(tfilename,'b');

       subplot(4,4,j)
      polyplot
      drawnow
set(gca,'XTick',[])           
set(gca,'YTick',[])
title(it)

   end
end