function p=plotraw(d,p)
% function p=plotraw(d,p)
% plot some results
%
%======================================================================
% P L O T R A W . M
% doc: Fri Jan 5 15:38:43 2007
% dlm: Wed Apr 15 11:30:50 2015
% (c) 2007 M. Visbeck with contribs from A. Thurnherr
% uE-Info: 19 76 NIL 0 0 72 0 2 4 NIL ofnI
%======================================================================
% MODIFICATIONS BY ANT:
% Jan 5, 2007: - fixed checkbeam() as suggested by B. Huber
% Jun 3, 2013: - BUG: top panel of Fig. 2 was wrong for dual-headed
% LADCPs with different UL/DL bin sizes
% Apr 15, 2015: - BUG: ax() is always used but only sometimes defined()
% reported by Achim Randelhoff (yahoo email 04/23/14)
pmax=200;
orient tall
if length(d.time_jul)>pmax
ii=fix(linspace(2,length(d.time_jul)-2,pmax));
else
ii=1:length(d.time_jul);
end
clf
subplot(411)
z=-d.izm(:,1)+d.z(1);
zz=[];
rw=[];
% find three beam solutions
% a profile is considered 3-beam is half of the good data are 3-beam
iz=fliplr(d.izu);
n3bu=0;
n3bd=0;
if length(iz)>1
zz=[zz;-z(iz)];
rw=[rw;d.rw(iz,ii)+d.weight(iz,ii)*0];
% check for 3-beam solutions
iw=sum(~isnan(d.rw(iz,ii)));
ie=sum(~isnan(d.re(iz,ii)));
i3bu=find(iw>(2*ie));
n3bu=length(i3bu)/length(ii)*100;
rw=[rw;rw(1,:)*NaN];
zz=[zz;0];
end
iz=d.izd;
if length(iz)>1
zz=[zz;-z(iz)];
rw=[rw;d.rw(iz,ii)+d.weight(iz,ii)*0];
% check for 3-beam solutions
iw=sum(~isnan(d.rw(iz,ii)));
ie=sum(~isnan(d.re(iz,ii)));
i3bd=find(iw>(2*ie));
n3bd=length(i3bd)/length(ii)*100;
end
% plot vertical velocities
col=jet(128);
col=([[1 1 1]; col]);
colormap(col)
pcolorn(ii,zz,rw) % this is neccessary even for dual-headed data sets
if (length(d.izu) > 0)
pcolorn(ii,zz(1:length(d.izu)),rw(1:length(d.izu),:))
hold on, shading flat
ax = axis;
ax(1) = ii(1);
ax(2) = ii(end);
ax(4) = -zz(end);
axis(ax);
pcolorn(ii,zz(length(d.izu)+1:end),rw(length(d.izu)+1:end,:))
else
ax = axis; % alternative solution to bug reported by Achim Randelhoff yahoo email 04/23/2014
end
hold on
plot([ii(1),ii(end)],[0 0],'-k')
colorbar('horiz')
% mark 3-beam solutions
if n3bu>10
l3b=NaN*ii;
l3b(i3bu)=ax(3);
plot(ii,l3b,'-r','linewidth',8)
text(mean(ax(1:2)),ax(4),['found ',int2str(n3bu),'% profiles 3 beam solutions'],...
'VerticalAlignment','top','HorizontalAlignment','center','fontsize',14)
end
if n3bd>10
l3b=NaN*ii;
l3b(i3bd)=ax(4);
plot(ii,l3b,'-r','linewidth',8)
text(mean(ax(1:2)),ax(3),['found ',int2str(n3bd),'% profiles 3 beam solutions'],...
'VerticalAlignment','bottom','HorizontalAlignment','center','fontsize',14)
end
ylabel('range [m]')
xlabel('ensemble')
if existf(p,'name')
streamer([p.name,' Figure 2']);
end
title(' W as function of bindepth and time')
if existf(d,'tsd_m')
subplot(427)
plot(d.tsd_m(1:length(d.izd),:),-d.zd)
axis tight
ax=axis;
if existf(d,'tsu_m')
hold on
plot(d.tsu_m(1:length(d.izu),:),d.zu)
axis tight
ax=axis;
plot(ax(1:2),ax(1:2)*0,'-k')
end
t=d.tsd_m;
checkbeam(t,ax,1)
if existf(d,'tsu_m')
t=d.tsu_m;
checkbeam(t,ax,0)
end
axis(ax)
ylabel('distance [m]')
xlabel('echo amplidute [dB]')
title('Beam Performance')
end
if existf(d,'cmd_m')
subplot(428)
plot(d.cmd_m(1:length(d.izd),:),-d.zd)
axis tight
ax=axis;
[dum,dum,dum,x,y]=makebars(-d.zd,sum(isfinite(d.weight(d.izd,:))'));
hold on
fill(-y/max(y)*10,x,'r')
if existf(d,'cmu_m')
hold on
plot(d.cmu_m(1:length(d.izu),:),d.zu)
plot(ax(1:2),ax(1:2)*0,'-k')
[dum,dum,dum,x,y]=makebars(d.zu,sum(isfinite(d.weight(d.izu,:))'));
fill(-y/max(y)*10,x,'g')
axis tight
ax=axis;
for i=1:4
[dum,ir]=min(abs(d.cmu_m(1:length(d.izu),i)-max(d.cmu_m(1,:))*0.3));
text((0.12*i+0.27)*ax(2),ax(4),int2str(d.zu(ir)),'VerticalAlignment','top')
p.up_range(i)=d.zu(ir);
end
text(0.01*ax(2),ax(4),['#',int2str(p.up_sn),...
' range:'],'VerticalAlignment','top')
end
ax(1)=-13;
for i=1:4
[dum,ir]=min(abs(d.cmd_m(1:length(d.izd),i)-max(d.cmd_m(1,:))*0.3));
p.dn_range(i)=d.zd(ir);
text((0.12*i+0.27)*ax(2),ax(3),int2str(d.zd(ir)),'VerticalAlignment','bottom')
end
text(0.01*ax(2),ax(3),['#',int2str(p.down_sn),...
' range:'],'VerticalAlignment','bottom')
axis(ax)
ylabel('distance [m]')
xlabel('correlation ')
title('Range of good data')
end
if existf(d,'z')
subplot(813)
plot(d.z)
ylabel('depth')
ax=axis;
ax(3)=min(d.z*1.05);
ax(4)=0;
ax(2)=length(d.z);
axis(ax)
end
if existf(d,'tilt')
subplot(814)
plot(d.tilt(1,:))
ylabel('tilt [deg]')
ax=axis;
ax(2)=length(d.z);
ax(4)=30;
axis(ax)
end
if existf(d,'hdg')
subplot(815)
plot(d.hdg(1,:))
ylabel('heading [deg]')
ax=axis;
ax(4)=360;
ax(2)=length(d.z);
axis(ax)
set(gca,'YTick',[0 90 180 270 360])
end
if existf(d,'xmv')
subplot(816)
plot(d.xmv(1,:))
text(length(d.xmv)/2,mean(d.xmv(1,:)),[' mean: ',num2str(fix(p.xmv(1)*10)/10)])
ylabel('X-mit volt down')
xlabel('ensemble')
ax=axis;
ax(2)=length(d.z);
axis(ax)
end
%=============================================
function checkbeam(t,ax,do)
% check beam performance
%%%bl=size(t,2);
bl = length(t); %%% BH fix
iend=fix(bl/2):bl;
tax=mean(ax(1:2));
if do
tay=ax(3);
tflag='bottom';
else
tay=ax(4);
tflag='top';
end
for i=1:4;
% first correct for source level
t(:,i)=t(:,i)-mean(t(:,i));
% find noise level
tn(i)=medianan(t(iend,i),2);
s2n(i)=mean(t(1:2,i))-tn(i);
end
ifail=s2n<max(s2n)*0.5;
ibad=~ifail & s2n<max(s2n)*0.65;
iweak=~ifail & ~ibad & s2n<max(s2n)*0.8;
for i=1:4;
text(ax(1)+0.2*i*diff(ax(1:2)),tay,[int2str(s2n(i)./max(s2n)*100),'%'],...
'VerticalAlignment',tflag)
end
if sum(ifail)>0
it=find(ifail==1);
text(tax,tay*0.5,[' beam ',int2str(it),' broken'],...
'VerticalAlignment',tflag,'HorizontalAlignment','center','fontsize',18)
end
if sum(ibad)>0
it=find(ibad==1);
text(tax,tay*0.65,[' beam ',int2str(it),' bad'],...
'VerticalAlignment',tflag,'HorizontalAlignment','center','fontsize',15)
end
if sum(iweak)>0
it=find(iweak==1);
text(tax,tay*0.8,[' beam ',int2str(it),' weak'],...
'VerticalAlignment',tflag,'HorizontalAlignment','center','fontsize',12)
end
return