  Results of large-scale numerical simulations are reported on the Anderson
localization in a two-dimensional square lattice tight-binding model with
random flux. Localization lengths, fluctuations of the conductance, and the
density of states are computed for quasi-one-dimensional geometry. Numerical
results indicate that the model exhibits the same critical behavior as the one
studied by Gade and Wegner. It is argued that all the states except a
zero-energy state are localized and the density of states has a singularity in
the center of the band. The energy scale below which the density of states
increases is found to be extremely small.
