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	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