
sn=5;
hrs=4;
hrs=input('hrs');
%hrs=str2num(h);
stuff2=jd(hrs);
[P,S]=polyfit(log(z(1:sn,hrs)),ws(1:sn,hrs)/usx(hrs),1);
zpl=[.001 20];	
uplx=polyval(P,log(zpl));
[Pn,Sn]=polyfit(log(z(1:sn,hrs)),wsn(1:sn,hrs)/usx(hrs),1);
uplxn=polyval(Pn,log(zpl));
      
%figure;
semilogy(ws(1:sn,hrs)/usx(hrs),z(1:sn,hrs),'x',uplx,zpl,wsn(1:sn,hrs)/usx(hrs),z(1:sn,hrs),'o',uplxn,zpl);
axis([0 50 .001 20]);
xlabel('U/u_*');
ylabel('z (m)');
title([stuff2]);
%axis([0 median(ws(1:sn,hrs))+3 .001 20]);