timour_matlab/mscripts/plotmeanw.m


figure
subplot(2,1,1)
plot(squeeze(meanw(Ny/2,Nx/2,1:length(mdep))),(mdep),'o-');
set(gca,'YLim',[mdep(length(mdep)) 0]);
xlabel('time mean w m/s')
ylabel('depth m')
grid
anchor1=meanw(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(meanw(:,Nx/2,i)),':');
end
plot(squeeze(meanw(:,Nx/2,2)),'--');
plot(squeeze(mean =w(:,Nx/2,Nz/2)),'-');
set(gca,'XLim',[0 Nx]);
grid on
ylabel('time mean w m/s')
xlabel('distance across the domain (gridpoints)')

subplot(2,1,1)
title(descriptor)
1
2	figure
3	subplot(2,1,1)
4	plot(squeeze(meanw(Ny/2,Nx/2,1:length(mdep))),(mdep),'o-');
5	set(gca,'YLim',[mdep(length(mdep)) 0]);
6	xlabel('time mean w m/s')
7	ylabel('depth m')
8	grid
9	anchor1=meanw(Nx/2,Ny/2,Nz/2);
10	text((anchor1*0.9),mdep(Nz/2),num2str(anchor1));
11	%anchor2=eval(it);
12	%text((anchor10.9),mdep(round(3Nz/4)),['timestep ' num2str(anchor2)] );
13
14	subplot(2,1,2)
15	for i=1:length(mdep)
16	hold on;plot(squeeze(meanw(:,Nx/2,i)),':');
17	end
18	plot(squeeze(meanw(:,Nx/2,2)),'--');
19	plot(squeeze(mean =w(:,Nx/2,Nz/2)),'-');
20	set(gca,'XLim',[0 Nx]);
21	grid on
22	ylabel('time mean w m/s')
23	xlabel('distance across the domain (gridpoints)')
24
25	subplot(2,1,1)
26	title(descriptor)