function Re_theta = Hall_Gibbings(lambda,Tu)
    % Hall and Gibbings 1972

    Re_theta = 190+exp(6.88-1.03*Tu);

end

