timour_matlab/mscripts/plotuv.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


cmin=20;
cmax=28;
V=[cmin cmax];

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

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

nz=5
caxis('manual')
%       colormap gray
%        cmap=flipud(colormap);
        cmap=zeros(64,3);
        colormap(cmap);

%        imagesc(lat,long,squeeze(t(:,:,nz))');shading flat;axis image;caxis(V);
        contour(squeeze(t(:,:,nz))',21:1:30);shading flat;axis image;caxis(V);
%       colorbar('horizontal');
        set(gca,'ydir','norm','FontSize',15,'XtickLabel','|','YTickLabel','|')
        xlabel('X','FontSize',20)
        ylabel('Y','FontSize',20,'Rotation',[0])
hold on

quiver(x,y,u(2:8:Nx,2:4:Ny-2,nz)',v(2:8:Nx,2:4:Ny-2,nz)')
% title(['velocity and temperature at the timestep ' num2str(eval(it)) ' level ' int2str(5) ]);
hold off
1	edhill	1.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			cmin=20;
16			cmax=28;
17			V=[cmin cmax];
18
19			counter=0;
20			for i=2:8:Nx,
21			counter=counter+1;
22			x(counter)=i;
23			end
24
25			counter=0;
26			for i=2:4:Ny,
27			counter=counter+1;
28			y(counter)=i;
29			end
30
31			nz=5
32			caxis('manual')
33			% colormap gray
34			% cmap=flipud(colormap);
35			cmap=zeros(64,3);
36			colormap(cmap);
37
38			% imagesc(lat,long,squeeze(t(:,:,nz))');shading flat;axis image;caxis(V);
39			contour(squeeze(t(:,:,nz))',21:1:30);shading flat;axis image;caxis(V);
40			% colorbar('horizontal');
41			set(gca,'ydir','norm','FontSize',15,'XtickLabel','\|','YTickLabel','\|')
42			xlabel('X','FontSize',20)
43			ylabel('Y','FontSize',20,'Rotation',[0])
44			hold on
45
46			quiver(x,y,u(2:8:Nx,2:4:Ny-2,nz)',v(2:8:Nx,2:4:Ny-2,nz)')
47			% title(['velocity and temperature at the timestep ' num2str(eval(it)) ' level ' int2str(5) ]);
48			hold off