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 |
|
24 |
sizeit=size(iterations); |
25 |
totalt=zeros(sizeit(1),1); |
26 |
figure;hold on |
27 |
for i=1:sizeit(1) |
28 |
tfilename=(['T.' iterations((i),1:10) ]); |
29 |
t=rdmeta(tfilename,'b'); |
30 |
totalt(i)=sum(sum(sum(t(:,:,:)))/(Nx*Ny*Nz*20); |
31 |
plot(eval(iterations((i),1:10)),totalt(i,:),'.');drawnow |
32 |
end |
33 |
|
34 |
|
35 |
|