
sn=5;
hrs=4;
%hrs=input('hrs');
for hrs=1:n,%cnt_mx,
   
	%hrs=str2num(h);
	stuff2=jd(hrs);
	[P,S]=polyfit(log(zs(1:sn,hrs)),wsx(1:sn,hrs)/usy(hrs)*.4,1);
	zpl=[.001 20];	
	uplx=polyval(P,log(zpl));
	[Pn,Sn]=polyfit(log(zs(1:sn,hrs))-psimx(1:sn,hrs),wsx(1:sn,hrs)/usy(hrs)*.4,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)*.4,zs(1:sn,hrs)/exp(psimx(1:sn,hrs)),'x');
  %semilogy(uplxn,zpl);
   if hrs==1,
		axis([0 30 .001 50]);
		xlabel('kU/(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=[2. 3.2 5.1 8.9 15];
%semilogy(umm,log(zz),'x')

