timour_matlab/mscripts/working1.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')
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			sizeit=size(iterations);
27			for i=1:sizeit(1)
28			iter(i)=eval(iterations(i,1:10));
29			end
30			nitstart=find(iter==eval(itstart))
31
32
33			figure;hold on
34
35			for i=nitstart:sizeit(1)-2
36			tfilename=(['T.' iterations((i),1:10) ]);
37			iterations((i),1:10)
38			i
39			t=rdmds(tfilename,'b');
40
41			meantee=mean(mean(mean(t(1:Nx/2,1:Ny-1,1))));
42			meantee1=mean(mean(mean(t(1:Nx-1,1:Ny-1,1:Nz-1))));
43
44			% tval(i)=meantee-meantee1;
45			tval(i)=meantee1-20;
46
47			% tval(i)=max(max(max(t)))-meantee1;
48			% plot(eval(iterations((i),1:10)),tval(i),'.');drawnow
49
50			% time(i)=eval(iterations(i,1:10));
51
52			if i==1
53			time(i)=100;
54			elseif i<=30
55			time(i)=time(i-1)+100;
56			else
57			time(i)=time(i-1)+50;
58			end
59
60			% time(i)=i*50;
61			plot(time(i),tval(i),'.');drawnow
62
63
64			end
65
66
67
68			tcol=tval;
69			title(descriptor)
70			xlabel('timestep')
71			ylabel('central temperature')
72			grid
73
74			computertime=time;
75
76			%save centralt centralt
77			%save iterations iterations
78			%save tcol tcol
79			%save tincol tincol
80			%save mdep mdep
81			save working tcol tincol mdep computertime