clear 
clear path 
whos
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 w3d
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

figure;hold on

load iterations
load tval

tvalanom=tval-20;
sizeit=size(iterations);
time=zeros(sizeit(1));
for i=1:sizeit(1)
iteration=eval([iterations((i),1:10) ]);
%time(i)=(iteration-su_its)*delt;
time(i)=(iteration-su_its)*delt/60;
end
for k=1:10
%plot(time*.60,tvalanom(:,k)+(10-k)-1,'g-');hold on
plot(time*.60,tvalanom(:,k),'g-');hold on
end
%plot(time,tvalanom(:,2),'b-');hold on
%plot(time,tvalanom(:,3),'g-');hold on


title(descriptor)
%xlabel('time (seconds)')
xlabel('time (minutes)')
ylabel('central temperature anomaly (degrees)')
grid