timour_matlab/mscripts/plotuv2.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,:,:);
cmin=20;
cmax=32;
V=[cmin cmax];

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

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

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

quiver(x,y,uu(2:4:Nx/2,2:4:Ny,5)',vv(2:4:Nx/2,2:4:Ny,5)')
title(['velocity and temperature at the timestep ' num2str(eval(it)) ' level ' int2str(5) ]);
hold off
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	cmin=20;
19	cmax=32;
20	V=[cmin cmax];
21
22	counter=0;
23	for i=Nx/2+2:4:Nx,
24	counter=counter+1;
25	x(counter)=i;
26	end
27
28	counter=0;
29	for i=2:4:Ny,
30	counter=counter+1;
31	y(counter)=i;
32	end
33
34	caxis('manual')
35	imagesc(x,y,squeeze(tt(:,:,5))');shading flat;axis image;caxis(V);colorbar('horizontal');
36	set(gca,'ydir','norm')
37	hold on
38
39	quiver(x,y,uu(2:4:Nx/2,2:4:Ny,5)',vv(2:4:Nx/2,2:4:Ny,5)')
40	title(['velocity and temperature at the timestep ' num2str(eval(it)) ' level ' int2str(5) ]);
41	hold off