We present a nonperturbative calculation of all multifractal scaling e
xponents at strong disorder for critical wave functions of Dirac fermi
ons interacting with a non-Abelian random vector potential in two dime
nsions. The results, valid for an arbitrary number of fermionic flavor
s, are obtained by deriving from conformal field theory an effective G
aussian model for the wave function amplitudes and mapping it to the t
hermodynamics of a single particle in a random potential. Our spectrum
confirms that the wave functions remain delocalized in the presence o
f strong disorder.