% Modified Lévêque correlation for near entrance of developing flows
function Nu = mLevNuL(Wall,Pr,Z)
  Pr=iff(Pr==0,1e-8,Pr);
  switch Wall
  case 'T'
    B=1.122; C=0.0472;              
  case {'H','H1','H2'}
    B=1.288; C=0.0206;              
  otherwise
    error('Wall= ''T'' or ''H''');
  end
  g=B*(1+(C./Pr(:)).^(3/4)).^(2/3); 
  fLReD=5.312./g./sqrt(Pr(:).*Z);   
  Nu=LevNuL(Wall,fLReD,Z);          
end
