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	edhill	1.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