
sn=5;
hrs=4;
%hrs=input('hrs');
for hrs=1:20,%cnt_mx,
   
	%hrs=str2num(h);
	stuff2=jd(hrs);
	[P,S]=polyfit(log(zs(1:sn,hrs)),wsx(1:sn,hrs)/usy(hrs),1);
	zpl=[.001 20];	
	uplx=polyval(P,log(zpl));
	[Pn,Sn]=polyfit(log(zs(1:sn,hrs)),wsn(1:sn,hrs)/usy(hrs),1);
	uplxn=polyval(Pn,log(zpl));
      
%figure;
	%semilogy(wsx(1:sn,hrs)/usy(hrs),zs(1:sn,hrs),'x',uplx,zpl,wsn(1:sn,hrs)/usy(hrs),z(1:sn,hrs),'o',uplxn,zpl);
	%semilogy(wsx(1:sn,hrs)/usy(hrs),zs(1:sn,hrs),'o');
  semilogy(uplxn,zpl);
   if hrs==1,
		axis([0 75 .001 20]);
		xlabel('U/u_*');
		ylabel('z (m)');
   	hold on;
   end;%axis([0 median(ws(1:sn,hrs))+3 .001 20]); 
end;
for j=1:5;umm(j)=median(wsn(j,:)./usy);end;
zz=mean(zs')%[2.1 3.2 5.1 8.9 14];
semilogy(umm,log(zz),'x')

