timour_matlab/mscripts/working2.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) ;
path   = this_path 
cmdstr=['cd ' path ];
eval(cmdstr);
path=pwd

iterations

% itstart = input(' Please enter start iteration : ','s')
itstart='0000002000'
sizeit=size(iterations);
for i=1:sizeit(1)
iter(i)=eval(iterations(i,1:10));
end
nitstart=find(iter==eval(itstart))


figure;hold on

for i=nitstart:sizeit(1)-2
tfilename=(['T.' iterations((i),1:10) ]);
 iterations((i),1:10)
 i
t=rdmds(tfilename,'b');

meantee=mean(mean(mean(t(1:Nx/2,1:Ny-1,1))));
meantee1=mean(mean(mean(t(1:Nx-1,1:Ny-1,1:Nz-1))));

% tval(i)=meantee-meantee1;
%  tval(i)=meantee1-20;

 tval(i)=max(max(max(t)))-meantee1;
% plot(eval(iterations((i),1:10)),tval(i),'.');drawnow

% time(i)=eval(iterations(i,1:10));

if i==1
time(i)=100;
elseif i<=30
time(i)=time(i-1)+100;
else
time(i)=time(i-1)+50;
end

% time(i)=i*50;
plot(time(i),tval(i),'.');drawnow


end


tcol=tval;
title(descriptor)
xlabel('timestep')
ylabel('central temperature')
grid

computertime=time;

%save centralt centralt
%save iterations iterations
%save tcol tcol
%save tincol tincol
%save mdep mdep 
save working tcol tincol mdep computertime
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
14
15	param_file_name = ...
16	input(' Please enter the name of the m-file with the parameters for this run : ','s') ;
17	feval(param_file_name) ;
18	path = this_path
19	cmdstr=['cd ' path ];
20	eval(cmdstr);
21	path=pwd
22
23	iterations
24
25	% itstart = input(' Please enter start iteration : ','s')
26	itstart='0000002000'
27	sizeit=size(iterations);
28	for i=1:sizeit(1)
29	iter(i)=eval(iterations(i,1:10));
30	end
31	nitstart=find(iter==eval(itstart))
32
33
34	figure;hold on
35
36	for i=nitstart:sizeit(1)-2
37	tfilename=(['T.' iterations((i),1:10) ]);
38	iterations((i),1:10)
39	i
40	t=rdmds(tfilename,'b');
41
42	meantee=mean(mean(mean(t(1:Nx/2,1:Ny-1,1))));
43	meantee1=mean(mean(mean(t(1:Nx-1,1:Ny-1,1:Nz-1))));
44
45	% tval(i)=meantee-meantee1;
46	% tval(i)=meantee1-20;
47
48	tval(i)=max(max(max(t)))-meantee1;
49	% plot(eval(iterations((i),1:10)),tval(i),'.');drawnow
50
51	% time(i)=eval(iterations(i,1:10));
52
53	if i==1
54	time(i)=100;
55	elseif i<=30
56	time(i)=time(i-1)+100;
57	else
58	time(i)=time(i-1)+50;
59	end
60
61	% time(i)=i*50;
62	plot(time(i),tval(i),'.');drawnow
63
64
65	end
66
67
68
69	tcol=tval;
70	title(descriptor)
71	xlabel('timestep')
72	ylabel('central temperature')
73	grid
74
75	computertime=time;
76
77	%save centralt centralt
78	%save iterations iterations
79	%save tcol tcol
80	%save tincol tincol
81	%save mdep mdep
82	save working tcol tincol mdep computertime