
sn=5;
hrs=4;
hrs=input('hrs');
%hrs=str2num(h);
stuff2=jd(hrs);
zxx=[.5 1 2 3 4 5 6 7 8 9 10 11 12 14 16 18 20];


[P,S]=polyfit(log(z(1:sn,hrs)),ws(1:sn,hrs),1);
zpl=[.001 20];	
uplx=polyval(P,log(zpl));
[Pn,Sn]=polyfit(log(z(1:sn,hrs)),wsn(1:sn,hrs),1);
uplxn=polyval(Pn,log(zpl));

upzz=polyval(Pn,log(zxx));
%figure;
%semilogy(ws(1:sn,hrs),z(1:sn,hrs),'x',uplx,zpl,wsn(1:sn,hrs),z(1:sn,hrs),'o',uplxn,zpl);
semilogy(wsn(1:sn,hrs),z(1:sn,hrs),'o',uplxn,zpl);

axis([0 15 .001 20]);
xlabel('U (m/s)');
ylabel('z (m)');
title([stuff2]);
%axis([0 median(ws(1:sn,hrs))+3 .001 20]);