function Bw_brutsaert = Bw_brutsaert(zu,L)
% correction function for stability, from Brutsaert and Parlange (1996)
% JGR (29,585-29,589)
% zu=height of wind measurement (m) ie. height of mixed layer 
% L=Obukhov length scale (m)
% 
%use surface layer values of a&b and the log profiles, eqn (7)
d0 = 0; 	%zero plane displacement
a = 0.374;
b = 2.408;

Bw_brutsaert = a*log(-(zu-d0)./L) + b;