No veo la contradicción que indicás. Para un número complejo z=a+jb, la fase se calcula como decís: phi=arctg(im(z)/re(z)). Sin excepciones.
Calculando la impedancia equivalente del paralelo de R y L, y aplicando lo anterior se llega al resultado que figura en la solución.
Calculando la impedancia equivalente del paralelo de R y L, y aplicando lo anterior se llega al resultado que figura en la solución.
Saludos,
Pablo