% Natural convection in which flow is returned by horizontal plate
function [h, NuL, RaL] = HorzPlate1NatNuT(Ts,Tinf,L,mat)
  [RaL Pr kf]  = Ra(Ts,Tinf,L,mat); 
  if min(RaL) < 1e4  warning = 'HorzPlate1NatNuT: RaL < 1e4'  
  end
  if max(RaL) > 1e11 warning = 'HorzPlate1NatNuT: RaL > 1e11' 
  end
  RaCrit=1e7;
  NuL=pick_fun(RaL,RaL<RaCrit,@(RaL)0.54*RaL.^(1/4),... 
                             @(RaL)0.15*RaL.^(1/3));    
  h    = NuL.*kf./L;                                    
end
