% Lévêque correlations for near entrance fully developed flows
function [Nu, A] = LevNuL(Wall,fReD,Z)
  switch Wall
  case 'T'
    A=0.40377;           
  case {'H','H1','H2'}
    A=0.43399;           
  otherwise
    error('Wall= ''T'' or ''H''');
  end 
  Nu=A*(fReD./Z).^(1/3); 
end
