I bosonize the position-space correlators of flavor-diagonal scalar fermion bilinears in arbitrary generalizations of the Schwinger model with n F massless fermions coupled to n A gauge bosons for n F ≥ n A . For n A = n F , the fermion bilinears can be bosonized in terms of n F scalars with masses proportional to the gauge couplings. As in the Schwinger model, bosonization can be used to find all correlators, including those that are forbidden in perturbation theory by anomalous chiral symmetries, but there are subtleties when there is more than one gauge boson. The new result here is the general treatment of the dependence on gauge boson masses in models with more than one gauge symmetry. For n A < n F , there are fermion bilinears with nontrivial anomalous dimensions and there are unbroken chiral symmetries so some correlators vanish while others are non-zero due to chiral anomlies. Taking careful account of the dependence on the masses, I show how the n A < n F models emerge from n A = n F as gauge couplings (and thus gauge boson masses) go to zero. When this is done properly, the limit of zero gauge coupling is smooth. Our consistent treatment of gauge boson masses guarantees that anomalous symmetries are broken while unbroken chiral symmetries are preserved because correlators that break the non-anomalous symmetries go to zero in the limit of zero gauge coupling.