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