function Re_theta = Dunham(lambda,Tu)
    % Dunham 1972

    Re_theta = (0.27+0.73*exp(-0.8*Tu)).*(550+680./(1+Tu-21*lambda));

end

