disp('med_hl_plt_AMMA');
%ysigoph=y1(:,133)';
clear hlmm hlbm nn dqm hlmmx hlbmx dum
uu=0:15;
for i=1:16
   ii=find(floor(v1)==uu(i));
   jj=find(isfinite(hll(ii)) & isfinite(yhlb(ii)) );
   if isempty(jj)
      hlmm(i)=NaN;
      hlbm(i)=NaN;
      hlmmx(i)=NaN;
      hlbmx(i)=NaN;
      nn(i)=0;
      dqm(i)=NaN;
      dum(i)=NaN;
   else
   	hlmm(i)=median(hll(ii(jj))+hl_webar(ii(jj))');
   	hlbm(i)=median(yhlb(ii(jj)));
   	nn(i)=length(jj);
      dqm(i)=median(yqs(ii(jj))-yq(ii(jj)));%
     	hlmmx(i)=mean(hll(ii(jj))+hl_webar(ii(jj))');
   	hlbmx(i)=mean(yhlb(ii(jj)));
		dum(i)=mean(v1(ii(jj)));
   end;
end;

figure;plot(dum,hlbm,'-x',dum,hlmm,'-o',dum,hlbmx,'x',dum,hlmmx,'o');xlabel('U (m/s)');ylabel('Hl (W/m^2)');title([jdollars  '  line-median  x-bulk  o-turb']);
figure;plot(hlbm,hlmm,'o',hlbmx,hlmmx,'x',[0 200],[0 200]);axis([0 200 0 200]);xlabel('Bulk Hl (W/m^2)');ylabel('Turb Hl (W/m^2)');title([jdollars  '  x-mean  o-median'  ]);



