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	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)