timour_matlab/mscripts/calctime.m

%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];


1	%clear
2	clear path
3
4	global Nx Ny Nz
5	global lat long dz dm mdep
6	global delt_su su_its t_su delt
7	global descriptor this_path
8	global f deltaf Q beta r_expt r_heat H
9	global time rots it
10	global g Cp rho_bar alpha
11	global u v t w
12	global iterations
13	param_file_name = ...
14	input(' Please enter the name of the m-file with the parameters for this run : ','s') ;
15	feval(param_file_name) ;
16	iterations
17	iterations1=iterations(1:88,1:10);
18	iterations2=iterations(89:length(iterations),1:10);
19	delt1=0.1;
20	delt2=0.05;
21	su_its=10000;
22	for i=1:length(iterations1)
23	time1(i)=(eval(iterations1(i,1:10))-su_its)*delt1;
24	end
25	for i=1:length(iterations2)
26	time2(i)=(eval(iterations2(i,1:10))-eval(iterations1(88,1:10)))*delt2+time1(88);
27	end
28	time=[time1 time2];
29
30
31
32
33
34
35