timour_matlab/mscripts/plottopflow.m

function flow
clear 
clear path 
whos
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

figure;k=1;
        uxy=u(:,:,k);
        ubx=(uxy(1:(Nx-1),:)+uxy(2:Nx,:))*0.5;
        ubx(Nx,:)=(uxy(Nx,:)+uxy(1,:))*0.5;
        uxy(find(uxy==0))=uxy(find(uxy==0))*NaN;
        ubx(find(ubx==0))=ubx(find(ubx==0))*NaN;
        % V
        vxy=v(:,:,k);
        vby=(vxy(:,1:(Ny-1))+vxy(:,2:Ny))*0.5;
        vby(:,Ny)=(vxy(:,1)+vxy(:,Ny))*0.5;
        vxy(find(vxy==0))=vxy(find(vxy==0))*NaN;
        vby(find(vby==0))=vby(find(vby==0))*NaN;

        % Scale arrow
        ubx(6,6) = 0.01;
        vby(6,6) = 0.0;
%        [x,y]=meshgrid(long(2:4:Ny-2),lat(2:4:Nx-2));
        [x,y]=meshgrid(long,lat);
%        quiver(x,y,ubx((2:4:Nx-2),(2:4:Ny-2))',vby((2:4:Nx-2),(2:4:Ny-2))',2); axis square
        quiver(x,y,ubx,vby,2); axis square
        title(['speed timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
        text(6,8,'1 cm/s')
        set(gca,'Xlim',[0,max(long)]);
        set(gca,'Ylim',[0,max(lat)]);
        drawnow
return
1	function flow
2	clear
3	clear path
4	whos
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	figure;k=1;
16	uxy=u(:,:,k);
17	ubx=(uxy(1:(Nx-1),:)+uxy(2:Nx,:))*0.5;
18	ubx(Nx,:)=(uxy(Nx,:)+uxy(1,:))*0.5;
19	uxy(find(uxy==0))=uxy(find(uxy==0))*NaN;
20	ubx(find(ubx==0))=ubx(find(ubx==0))*NaN;
21	% V
22	vxy=v(:,:,k);
23	vby=(vxy(:,1:(Ny-1))+vxy(:,2:Ny))*0.5;
24	vby(:,Ny)=(vxy(:,1)+vxy(:,Ny))*0.5;
25	vxy(find(vxy==0))=vxy(find(vxy==0))*NaN;
26	vby(find(vby==0))=vby(find(vby==0))*NaN;
27
28	% Scale arrow
29	ubx(6,6) = 0.01;
30	vby(6,6) = 0.0;
31	% [x,y]=meshgrid(long(2:4:Ny-2),lat(2:4:Nx-2));
32	[x,y]=meshgrid(long,lat);
33	% quiver(x,y,ubx((2:4:Nx-2),(2:4:Ny-2))',vby((2:4:Nx-2),(2:4:Ny-2))',2); axis square
34	quiver(x,y,ubx,vby,2); axis square
35	title(['speed timestep ' num2str(eval(it)) ' level ' int2str(k) ]);
36	text(6,8,'1 cm/s')
37	set(gca,'Xlim',[0,max(long)]);
38	set(gca,'Ylim',[0,max(lat)]);
39	drawnow
40	return