timour_matlab/mscripts/plottevolution2.m

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

1	edhill	1.1	clear
2			clear path
3			whos
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 w3d
12			global iterations
13
14			param_file_name = ...
15			input(' Please enter the name of the m-file with the parameters for this run : ','s') ;
16			feval(param_file_name) ;
17			path = this_path
18			cmdstr=['cd ' path ];
19			eval(cmdstr);
20			path=pwd
21
22			figure;hold on
23
24			load iterations
25			load tval
26
27			tvalanom=tval-20;
28			sizeit=size(iterations);
29			time=zeros(sizeit(1));
30			for i=1:sizeit(1)
31			iteration=eval([iterations((i),1:10) ]);
32			%time(i)=(iteration-su_its)*delt;
33			time(i)=(iteration-su_its)*delt/60;
34			end
35			for k=1:10
36			%plot(time*.60,tvalanom(:,k)+(10-k)-1,'g-');hold on
37			plot(time*.60,tvalanom(:,k),'g-');hold on
38			end
39			%plot(time,tvalanom(:,2),'b-');hold on
40			%plot(time,tvalanom(:,3),'g-');hold on
41
42
43			title(descriptor)
44			%xlabel('time (seconds)')
45			xlabel('time (minutes)')
46			ylabel('central temperature anomaly (degrees)')
47			grid
48