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	edhill	1.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)