timour_matlab/mscripts/plottwfields.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(t(:,:,k)'));shading flat; colorbar; axis square
title(['t timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
drawnow
subplot(2,2,2)
pcolor(squeeze(w(:,:,k)'));shading flat; colorbar; axis square
title(['w timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
drawnow

pause

%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
end

return
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(t(:,:,k)'));shading flat; colorbar; axis square
21	title(['t timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
22	drawnow
23	subplot(2,2,2)
24	pcolor(squeeze(w(:,:,k)'));shading flat; colorbar; axis square
25	title(['w timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
26	drawnow
27
28	pause
29
30	%subplot(2,2,3)
31	%pcolor(squeeze((v(:,:,k)'.^2+u(:,:,k)'.^2).^0.5));shading flat; colorbar; axis square
32	%title(['speed timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
33	%drawnow
34	end
35
36	return