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