% Natural convection from a vertical plate of constant temperature 
function [h, NuL, RaL] = VertPlateNatNuT(Ts,Tinf,L,RaCrit,mat)
  [RaL, Pr, kf] = Ra(Ts,Tinf,L,mat);
  Cl=@(Pr)(Pr./4./(0.609+1.221*Pr.^0.5+1.238*Pr)).^.25;
  NuLamBL=@(RaL,Pr)Cl(Pr).*RaL.^0.25;                    
  NuLam=@(RaL,Pr)2.0./log(1+2.0./NuLamBL(RaL,Pr));       
  Ct=@(Pr)0.13*Pr.^0.22./(1+0.61*Pr.^0.81).^0.42;
  NuTurb=@(RaL,Pr)Ct(Pr).*RaL.^(1/3);                    
  NuL = MixNu('T',NuLam,NuTurb,RaCrit,RaL,Pr,0);         
  h = NuL.*kf./L;
end
