% Laminar pipe flow with combined thermal and hydrodynamic entry region. 
function [h, NuL, f, ReD] = LamPipeEntryNuH(um,D,L,mat,Tm)
  [Z, ReD, Pr, km] = InvGz(um,D,L,mat,Tm);           
  if max(ReD) > 2300 warning = 'LamPipeEntryNuH: ReD > 2300'
  end
  if Pr < 0.6        warning = 'LamPipeEntryNuH: Pr < 0.6'
  end
  f   = (64./ReD).*(1+(D.*ReD./21.7./L).^1.1).^0.46; 
  NuL = LamEntryNuL('T',Pr,64,4.36,Z);               
  h   = NuL*km./D;
end
