hrs1=1;%find(jd==508);
hrs2=8064;%find(jd==514);
%hrs1=110;hrs2=111;
sn=5;
hrs=4;
wsy=[ws1 ws2 ws3 ws4 ws5]';
zt=[z1 z2 z3 z4 z5 ]';
tvy=[tv1 tv2 tv3 tv4 tv5]';
ty=[tcor1 tcor2 tcor3 tcor4 tcor5]';
qvy=[q1 q2 q3 q4 q5]';


%hrs=input('hrs');
clear wsx; clear zs;clear slopt;clear slopq;clear jdt;
jamx=0;
hold off;
zref=5;%(2*3.2*5.1*9*17)^.2;
wt=hsm./rhoa/cpa;
tsy=-wt./us;
Lt=-(tt+tdk)*2.5./wt.*us.^3/grav;%sensbile heat flux L
zpl=[1 2 4 6 8 10 15 20];
for hrs=hrs1:hrs2,%cnt_mx,
   ii=find(isfinite(tvy(:,hrs)));
   [sn a]=size(ii);
   if sn>3
      zs=zt(ii,hrs);
   	tvx=tvy(ii,hrs);
		[P,S]=polyfit(log(zs),tvx,2);
	%	zpl=[1 zs' 20];	
		tplx=polyval(P,log(zpl));
      slopt(hrs)=P(2)+2*P(1)*log(zref);
      
      qvx=qvy(ii,hrs);
		[P,S]=polyfit(log(zs),qvx,2);
		qplx=polyval(P,log(zpl));
      slopq(hrs)=P(2)+2*P(1)*log(zref);
      
     jdt(hrs)=jd(hrs);
      
%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);
%   figure;semilogy(tvx./tsy(hrs)*.4,zs,'o',tvx./tsy(hrs)*.4+psimx,zs,'x',uplx,zpl,uplxn,zpl);
	%semilogy(tvx,zs,'o',tplx,zpl);

  %semilogy(uplxn,zpl);
   	if jamx==1,
         axis([-15 5 1 20]);
         xlabel('Tv (C)');
			ylabel('z (m)');
         hold on;
         jamx=0;
      end;%axis([0 median(ws(1:sn,hrs))+3 .001 20]); 
   else;
      slopt(hrs)=NaN;
      slopq(hrs)=NaN;
      jdt(hrs)=jd(hrs);

	end;   
end;


figure;plot(jd,-(qss'-q4),jdt,slopq*5,'o',[508 515],[0 0 ]);axis([508 514 -2 2]);xlabel('Julian Day (1997)');ylabel('10*dq/dz at z=5 m');
figure;plot(jd,-(tscv-tv4),jdt,slopt*5,'o',[508 515],[0 0 ]);axis([508 514 -5 5]);xlabel('Julian Day (1997)');ylabel('10*dThv/dz at z=5 m');

figure;xm=-1;xx=2;plot(slopt*5,hsbm./median(wsy),'.',[xm xx],[0 0]);axis([xm xx -6 5]);
ylabel('Hsbm/wspd');xlabel('dThv/dz*5');

figure;xm=-1;xx=2;plot(slopt*5,hs./median(wsy),'.',[xm xx],[0 0]);axis([xm xx -6 5]);
ylabel('Hs/wspd');xlabel('dThv/dz*5');


figure;xm=-1;xx=2;plot(slopt*5,hsm./median(wsy),'.',[xm xx],[0 0]);axis([xm xx -6 5]);
ylabel('Hsm/wspd');xlabel('dThv/dz*5');

figure;xm=-2;xx=4;plot(-dtt,hsm./median(wsy),'.',[xm xx],[0 0]);axis([xm xx -6 5]);
ylabel('Hsm/wspd');xlabel('T-Ts');


figure;xm=-.5;xx=.5;plot(slopq*5,hlbm./median(wsy),'.',[xm xx],[0 0]);axis([xm xx -2 2]);
ylabel('Hlbm/wspd');xlabel('dq/dz*5');

figure;xm=-.5;xx=.5;plot(slopq*5,hl./median(wsy),'.',[xm xx],[0 0]);axis([xm xx -2 2]);
ylabel('Hl/wspd');xlabel('dq/dz*5');