0001 function [znm,vm] = VMT_ComputeNormProf(zn,v,n)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013 zn_vec = reshape(zn,size(zn,1)*size(zn,2),1);
0014 v_vec = reshape(v,size(v,1)*size(v,2),1);
0015
0016
0017 cell_breaks = 0:1/n:1;
0018 dum = diff(cell_breaks);
0019 znm = cumsum([(cell_breaks(2)/2) dum(2:end)]);
0020
0021 vm = nan*ones(1,n);
0022 obs = nan*ones(1,n);
0023 for i = 1:n
0024 indx = find(zn_vec >= cell_breaks(i) & zn_vec < cell_breaks(i+1));
0025 vm(i) = nanmedian(v_vec(indx));
0026 obs(i) = length(indx);
0027 end
0028
0029
0030
0031 indx = find(obs < 0.2*nanmedian(obs));
0032
0033
0034
0035 if 0
0036 figure(1); clf
0037 plot(v_vec,zn_vec,'k.'); hold on
0038 plot(vm,znm,'bs-'); hold on
0039 plot(vm(indx),znm(indx),'ro')
0040 xlabel('velocity')
0041 ylabel('normalized height above bottom')
0042 xlim([0 max(v_vec)])
0043 ylim([0 1])
0044 end
0045
0046
0047 vm(indx) = nan;
0048
0049
0050
0051
0052
0053