timour_matlab/mscripts/plotuvxz.m

function plotfields( iter )
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


for k=1:Nz
k
figure(1)
subplot(2,2,1)
pcolor(squeeze(u(:,:,k)'));shading flat; colorbar; axis square
title(['u timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
drawnow
subplot(2,2,2)
pcolor(squeeze(v(:,:,k)'));shading flat; colorbar; axis square
title(['v timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
drawnow
subplot(2,2,3)
pcolor(squeeze((v(:,:,k)'.^2+u(:,:,k)'.^2).^0.5));shading flat; colorbar; axis square
title(['speed timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
drawnow

pause

end

return
1	edhill	1.1	function plotfields( iter )
2			clear
3			clear path
4
5			global Nx Ny Nz
6			global lat long dz dm mdep
7			global delt_su su_its t_su delt
8			global descriptor this_path
9			global f deltaf Q beta r_expt r_heat H
10			global time rots it
11			global g Cp rho_bar alpha
12			global u v t w
13			global iterations
14
15
16			for k=1:Nz
17			k
18			figure(1)
19			subplot(2,2,1)
20			pcolor(squeeze(u(:,:,k)'));shading flat; colorbar; axis square
21			title(['u timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
22			drawnow
23			subplot(2,2,2)
24			pcolor(squeeze(v(:,:,k)'));shading flat; colorbar; axis square
25			title(['v timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
26			drawnow
27			subplot(2,2,3)
28			pcolor(squeeze((v(:,:,k)'.^2+u(:,:,k)'.^2).^0.5));shading flat; colorbar; axis square
29			title(['speed timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
30			drawnow
31
32			pause
33
34			end
35
36			return