wrong normalization of spinor states in states_elec_generate_random

Wave functions are normalized wrongly by states_elec_generate_random due to:

        alpha = TOCMPLX(sqrt(M_HALF + st%spin(3, ist, ik)), M_ZERO)
        beta  = TOCMPLX(sqrt(M_ONE - abs(alpha)**2), M_ZERO)
        if(abs(alpha) > M_ZERO) then
          beta = TOCMPLX(st%spin(1, ist, ik) / abs(alpha), st%spin(2, ist, ik) / abs(alpha))
        end if
        zpsi(1:mesh%np, 1) = alpha*zpsi(1:mesh%np, 1)
        zpsi(1:mesh%np, 2) = beta*zpsi(1:mesh%np, 2)