Programs for modeling binary star observables compute emergent intensity for given composition as it varies with local effective temperature, local gravity, and direction. With arrival of huge data sets from Gaia and other surveys, the benefits of fast, compact, and accurate computation of atmospheric radiation is likely to remain critical for the foreseeable future. Experience has shown that accurate radiative modeling is important for good parameter estimation. Here we augment the radiative treatment by Van Hamme & Wilson (2003) with a procedure by which individuals can generate the needed Legendre coefficients for arbitrary photometric bands. Resulting files can be inserted directly into the Wilson-Devinney (W-D) program without sacrificing portability or program unity, and should easily be adaptable to other binary star programs. We expect the new bandpass options to become part of the public W-D program. Limb-darkening tables will be placed at http://www.fiu.edu/~vanhamme/limdark.htm.