disp('bulk')
fclose('all');
flist=fopen('\matlab5\cwf\progs\wavhr2_5.txt','r');
stuff=fgetl(flist);
stuff=fgetl(flist);
nx=200;
zu=15;
zt=15;
zq=15;
	
for i = 1:nx
	P=1010;                     %air pressure
	us=0;                       %surface current
	zi=600;                    %inversion ht
	stuff=fgetl(flist);
	u=eval(stuff(16:21));
	ts=eval(stuff(38:42));
	t=eval(stuff(43:48));
	qs=eval(stuff(49:54));
	q=eval(stuff(55:60));
	Rs=eval(stuff(112:117));
	Rl=eval(stuff(118:123));
	rain=eval(stuff(124:129));
	x=[u us ts t qs q Rs Rl rain zi  P zu zt zq] ;
	y=cool_lkb(x);
	hsb=y(1);
	hlb=y(2);
	tau=y(3);
	zo=y(4);
	zot=y(5);
	zoq=y(6);
	L=y(7);
	usr=y(8);
	tsr=y(9);
	qsr=y(10);
	dter=y(11);
	dqer=y(12);
	tkt=y(13);
	RF=y(14);
	Cd=y(15);
	Ch=y(16);
	Ce=y(17);
	Cdn_10=y(18);
	Chn_10=y(19);
	Cen_10=y(20);
	e=[stuff(16:21) stuff(38:60) stuff(112:123) stuff(71:75) stuff(86:90) stuff(105:111) '  ' num2str(hsb) '  ' num2str(hlb) '  ' num2str(tau) '  ' num2str(dter) '  ' num2str(tkt) '  ' num2str(Cen_10)];
	disp(e)
end; %  data line loop
	
	
	
	
	
	
	
	