% Correlation for entry convection local Nusselt number
function [Nux, NuL] = LamEntryNux(Wall,Pr,fReD,NuInf,Z)
  [NuL, p, S]=LamEntryNuL(Wall,Pr,fReD,NuInf,Z);            
  [GzNux, GzNuL]=GrtzNux(Wall,fReD,NuInf,Z);                
  LevNu=LevNuL(Wall,fReD,Z);                                
  mLevNu=mLevNuL(Wall,Pr,Z);                                
  F=(LevNu./mLevNu).*(1+(S*Pr*Z).^(p/6)).^(1/p);            
  Add=F.*(2-(1+(S*Pr*Z).^(p/6)).^-1)/3./sinh(2*F);
  Add(isnan(Add))=0;                                        % if Pr=Inf

  switch Wall
  case 'T'
    Nux=NuL.*(GzNux./GzNuL-Add);                            
  case {'H','H1','H2'}
    Nux=NuL./(GzNuL./GzNux+Add);                            
  otherwise
    error('Wall= ''T'' or ''H''');
  end
end
