function Re_theta = Davis_Atkins(lambda,Tu)
    % Davis and Atkins 2024

    Re_theta = 440*Tu.^(-0.56);

end

