1 |
edhill |
1.1 |
clear |
2 |
|
|
clear path |
3 |
|
|
clear iterations |
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 |
|
|
figure;hold on |
24 |
|
|
%set(gca,'XLim',[10000 35000]) |
25 |
|
|
%set(gca,'YLim',[20 35]) |
26 |
|
|
sizeit=size(iterations); |
27 |
|
|
centraltheta=zeros(sizeit(1),Nz); |
28 |
|
|
|
29 |
|
|
for i=1:sizeit(1) |
30 |
|
|
tfilename=(['T.' iterations((i),1:10) ]); |
31 |
|
|
t=rdmeta(tfilename,'b'); |
32 |
|
|
|
33 |
|
|
|
34 |
|
|
centralt(i,:)=t(Nx/2,Ny/2,:)-20; |
35 |
|
|
plot(eval(iterations((i),1:10)),centralt(i,:),'.');drawnow |
36 |
|
|
tincol(i)=sum((centralt(i,:)).*dz); |
37 |
|
|
time(i)=(eval(iterations(i,1:10))-su_its)*delt; |
38 |
|
|
|
39 |
|
|
|
40 |
|
|
end |
41 |
|
|
tcol=centralt; |
42 |
|
|
title(descriptor) |
43 |
|
|
xlabel('timestep') |
44 |
|
|
ylabel('central temperature') |
45 |
|
|
grid |
46 |
|
|
|
47 |
|
|
computertime=time; |
48 |
|
|
|
49 |
|
|
%save centralt centralt |
50 |
|
|
%save iterations iterations |
51 |
|
|
%save tcol tcol |
52 |
|
|
%save tincol tincol |
53 |
|
|
%save mdep mdep |
54 |
|
|
save computer tcol tincol mdep computertime |