1 |
clear |
2 |
clear path |
3 |
whos |
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 w3d |
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 |
|
25 |
sizeit=size(iterations); |
26 |
centralw=zeros(sizeit(1),Nz+1); |
27 |
|
28 |
for i=1:sizeit(1) |
29 |
ufilename=(['U.' iterations((i),1:10) ]); |
30 |
vfilename=(['V.' iterations((i),1:10) ]); |
31 |
u=rdmeta(ufilename,'b'); |
32 |
v=rdmeta(vfilename,'b'); |
33 |
w=zeros(Nx,Ny,Nz+1); |
34 |
for k=Nz:-1:1, |
35 |
w(1:Nx-1,1:Ny-1,k)= w(1:Nx-1,1:Ny-1,k+1)... |
36 |
-u(2:Nx,1:Ny-1,k)*dz(k)/dm +u(1:Nx-1,1:Ny-1,k)*dz(k)/dm ... |
37 |
-v(1:Nx-1,2:Ny,k)*dz(k)/dm +v(1:Nx-1,1:Ny-1,k)*dz(k)/dm ; |
38 |
end |
39 |
centralw(i,:)=w(Nx/2,Ny/2,:); |
40 |
plot(eval(iterations((i),1:10)),centralw(i,Nz/2),'o');drawnow |
41 |
end |
42 |
|
43 |
title(descriptor) |
44 |
xlabel('timestep') |
45 |
ylabel('central mid-depth w (m/s)') |
46 |
grid |
47 |
|