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) ;
path   = this_path 
cmdstr=['cd ' path ];
eval(cmdstr);
path=pwd

iterations

% itstart = input(' Please enter start iteration : ','s')
itstart='0000002000'
sizeit=size(iterations);
for i=1:sizeit(1)
iter(i)=eval(iterations(i,1:10));
end
nitstart=find(iter==eval(itstart))


figure;hold on

for i=nitstart:sizeit(1)-2
tfilename=(['T.' iterations((i),1:10) ]);
 iterations((i),1:10)
 i
t=rdmds(tfilename,'b');

meantee=mean(mean(mean(t(1:Nx/2,1:Ny-1,1))));
meantee1=mean(mean(mean(t(1:Nx-1,1:Ny-1,1:Nz-1))));

% tval(i)=meantee-meantee1;
%  tval(i)=meantee1-20;

 tval(i)=max(max(max(t)))-meantee1;
% plot(eval(iterations((i),1:10)),tval(i),'.');drawnow

% time(i)=eval(iterations(i,1:10));

if i==1
time(i)=100;
elseif i<=30
time(i)=time(i-1)+100;
else
time(i)=time(i-1)+50;
end

% time(i)=i*50;
plot(time(i),tval(i),'.');drawnow


end



tcol=tval;
title(descriptor)
xlabel('timestep')
ylabel('central temperature')
grid

computertime=time;

%save centralt centralt
%save iterations iterations
%save tcol tcol
%save tincol tincol
%save mdep mdep 
save working tcol tincol mdep computertime