timour_matlab/mscripts/plotuv1.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


tt(1:Nx/2,:,:)=t(Nx/2+1:Nx,:,:);
uu(1:Nx/2,:,:)=u(Nx/2+1:Nx,:,:);
vv(1:Nx/2,:,:)=v(Nx/2+1:Nx,:,:);

%tt(1:Nx/2,:,:)=t(Nx/2+1:Nx,:,:);
%uu(1:Nx/2,:,:)=u(Nx/2+1:Nx,:,:);
%vv(1:Nx/2,:,:)=v(Nx/2+1:Nx,:,:);

cmin=20;
cmax=32;
V=[cmin cmax];

counter=0;
for i=2:4:Nx,
counter=counter+1;
x(counter)=i;
end

counter=0;
for i=2:4:Ny,
counter=counter+1;
y(counter)=i;
end

for k=1:Nz,
figure(1)

caxis('manual')
        imagesc(x,y,squeeze(t(1:Nx,:,k))');shading flat;axis image;caxis(V);colorbar('horizontal');
        set(gca,'ydir','norm')
hold on

quiver(x,y,u(2:4:Nx,2:4:Ny,k)',v(2:4:Nx,2:4:Ny,k)')
title(['velocity and temperature at the timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
hold off

drawnow
pause
end

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			tt(1:Nx/2,:,:)=t(Nx/2+1:Nx,:,:);
16			uu(1:Nx/2,:,:)=u(Nx/2+1:Nx,:,:);
17			vv(1:Nx/2,:,:)=v(Nx/2+1:Nx,:,:);
18
19			%tt(1:Nx/2,:,:)=t(Nx/2+1:Nx,:,:);
20			%uu(1:Nx/2,:,:)=u(Nx/2+1:Nx,:,:);
21			%vv(1:Nx/2,:,:)=v(Nx/2+1:Nx,:,:);
22
23			cmin=20;
24			cmax=32;
25			V=[cmin cmax];
26
27			counter=0;
28			for i=2:4:Nx,
29			counter=counter+1;
30			x(counter)=i;
31			end
32
33			counter=0;
34			for i=2:4:Ny,
35			counter=counter+1;
36			y(counter)=i;
37			end
38
39			for k=1:Nz,
40			figure(1)
41
42			caxis('manual')
43			imagesc(x,y,squeeze(t(1:Nx,:,k))');shading flat;axis image;caxis(V);colorbar('horizontal');
44			set(gca,'ydir','norm')
45			hold on
46
47			quiver(x,y,u(2:4:Nx,2:4:Ny,k)',v(2:4:Nx,2:4:Ny,k)')
48			title(['velocity and temperature at the timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
49			hold off
50
51			drawnow
52			pause
53			end
54