timour_matlab/mscripts/plotw.m

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
subplot(2,1,1)
plot(squeeze(w(Nx/2,Ny/2,1:Nz)),'o-');
%set(gca,'YLim',[mdep(length(mdep)) 0]);
xlabel('depth , gridpoints')
ylabel('w m/s')
grid

figure

cmin=-1.e-4;
cmax=1.e-4,
V=[cmin cmax];

title='in-situ vertical velocity at the 4rd level from the top';
        imagesc(lat,long,squeeze(w(:,:,4))');shading flat;axis image;caxis(V);colorbar('vertical');
        set(gca,'ydir','norm')
        text(0,-30,descriptor);
        text(0,130,title);
        text(0,-20,'timestep');text(50,-20,num2str(it))


%anchor1=w(Nx/2,Ny/2,Nz/2);
%text((anchor1*0.9),mdep(Nz/2),num2str(anchor1));
%anchor2=eval(it);
%text((anchor1*0.9),mdep(round(3*Nz/4)),['timestep  ' num2str(anchor2)] );

%subplot(2,1,2)
%for i=1:length(mdep)
%hold on;plot(squeeze(w(:,Nx/2,i)),':');
%end
%plot(squeeze(w(:,Nx/2,2)),'--');
%plot(squeeze(w(:,Nx/2,Nz/2)),'-');
%set(gca,'XLim',[0 Nx]);
%grid on
%ylabel('w m/s')
%xlabel('distance across the domain (gridpoints)')
%
%subplot(2,1,1)
%title(descriptor)
1	edhill	1.1	clear
2			clear path
3			whos
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			figure
15			subplot(2,1,1)
16			plot(squeeze(w(Nx/2,Ny/2,1:Nz)),'o-');
17			%set(gca,'YLim',[mdep(length(mdep)) 0]);
18			xlabel('depth , gridpoints')
19			ylabel('w m/s')
20			grid
21
22			figure
23
24			cmin=-1.e-4;
25			cmax=1.e-4,
26			V=[cmin cmax];
27
28			title='in-situ vertical velocity at the 4rd level from the top';
29			imagesc(lat,long,squeeze(w(:,:,4))');shading flat;axis image;caxis(V);colorbar('vertical');
30			set(gca,'ydir','norm')
31			text(0,-30,descriptor);
32			text(0,130,title);
33			text(0,-20,'timestep');text(50,-20,num2str(it))
34
35
36			%anchor1=w(Nx/2,Ny/2,Nz/2);
37			%text((anchor1*0.9),mdep(Nz/2),num2str(anchor1));
38			%anchor2=eval(it);
39			%text((anchor10.9),mdep(round(3Nz/4)),['timestep ' num2str(anchor2)] );
40
41			%subplot(2,1,2)
42			%for i=1:length(mdep)
43			%hold on;plot(squeeze(w(:,Nx/2,i)),':');
44			%end
45			%plot(squeeze(w(:,Nx/2,2)),'--');
46			%plot(squeeze(w(:,Nx/2,Nz/2)),'-');
47			%set(gca,'XLim',[0 Nx]);
48			%grid on
49			%ylabel('w m/s')
50			%xlabel('distance across the domain (gridpoints)')
51			%
52			%subplot(2,1,1)
53			%title(descriptor)