% Correlation for entry convection average Nusselt number
function [NuL, p, S] = LamEntryNuL(Wall,Pr,fReD,NuInf,Z)
  GzNu=GrtzNuL(Wall,fReD,NuInf,Z);                    
  LevNu=LevNuL(Wall,fReD,Z);                          
  mLevNu=mLevNuL(Wall,Pr,Z);                          

  S=exp(NuInf*1.78-7.42);                             
  p=0.511*NuInf.^2-4.06*NuInf+9.75;                   
  F=(LevNu./mLevNu).*(1+(S*Pr*Z).^(p/6)).^(1/p);      
  NuL=GzNu./tanh(F);                                  
end
