% Forced convection with jet impingement
function [h, NuD, ReD] = JetNuT(U,D,H,R,mat,Tf)
  [ReD,Pr,kf] = Re(U,D,mat,Tf);                       
  H_D  = H./D; R_D  = R./D;
  if     min(H_D) < 2   warning = 'JetNuT: H/D < 2'   
  elseif max(H_D) > 12  warning = 'JetNuT: H/D > 12'  
  end
  if     min(R_D) < 2.5 warning = 'JetNuT: R/D < 2.5' 
  elseif max(R_D) > 7.5 warning = 'JetNuT: R/D > 7.5' 
  end
  if     min(ReD) < 2e3 warning = 'JetNuT: ReD < 2e3' 
  elseif max(ReD) > 4e5 warning = 'JetNuT: ReD > 4e5' 
  end
  G    = (1-1.1./R_D)./R_D./(1+0.1*(H_D-6)./R_D) ;    
  F1   = 2*(ReD.*(1+0.005*ReD.^0.55)).^0.5;           
  NuD  = G.*F1*Pr^0.42;                               
  h    = NuD*kf./D;                                   
end
