%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) ;
iterations
iterations1=iterations(1:88,1:10);
iterations2=iterations(89:length(iterations),1:10);
delt1=0.1;
delt2=0.05;
su_its=10000;
for i=1:length(iterations1)
time1(i)=(eval(iterations1(i,1:10))-su_its)*delt1;
end
for i=1:length(iterations2)
time2(i)=(eval(iterations2(i,1:10))-eval(iterations1(88,1:10)))*delt2+time1(88);
end
time=[time1 time2];