% Fully developed laminar pipe flow with constant wall temperature. 
function [h, NuL, f, ReD] = LamPipeNuT(um,D,L,mat,Tm)
  [Z, ReD, Pr, km] = InvGz(um,D,L,mat,Tm);   
  if max(ReD) > 2300 warning = 'LamPipeNuT: ReD > 2300'
  end
  f   = 64./ReD;                            % Table 11-1
  NuL = GrtzNuL('T',64,3.66,Z);             
  h   = NuL.*km./D;
end
