We analyze the behavior of rational inner functions on the unit bidisk near singularities on the distinguished boundary T-2 using level sets. We show that the unimodular level sets of a rational inner function phi can be parametrized with analytic curves and connect the behavior of these analytic curves to that of the zero set of phi. We apply these results to obtain a detailed description of the fine numerical stability of phi: for instance, we show that partial derivative phi/partial derivative z(1) and partial derivative phi/partial derivative z(2) always possess the same L-p-integrability on T-2, and we obtain combinatorial relations between intersection multiplicities at singularities and vanishing orders for branches of level sets. We also present several new methods of constructing rational inner functions that allow us to prescribe properties of their zero sets, unimodular level sets, and singularities.