% Natural convection in which flow passes horizontal plate
function [h, NuL, RaL] = HorzPlate2NatNuT(Ts,Tinf,L,mat)
  [RaL Pr kf]  = Ra(Ts,Tinf,L,mat); 
  if min(RaL) < 1e5  warning = 'HorzPlate2NatNuT: RaL < 1e5'  
  end
  if max(RaL) > 1e10 warning = 'HorzPlate2NatNuT: RaL > 1e10' 
  end
  NuL  = 0.52*RaL.^.2;                                  
  h    = NuL.*kf./L;                                    
end
