% Convection from rotating disk
function [hT, NuT, ReR] = RotDiskNuT(Omega,R,ReCrit,mat,Tf)
  [ReR, Pr, kf] = Re(Omega*R,R,mat,Tf);                             

  NuLam=@(ReR,Pr)0.585*ReR.^0.5./(0.6./Pr+0.95./Pr.^(1/3));         
  NuTurb=@(ReR,Pr)0.016154*ReR.^0.8.*Pr.^(1/3);                     

  NuT = MixNu('T',NuLam,NuTurb,ReCrit,ReR,Pr,1/2);                  
  hT = NuT.*kf./R;
end
