try fitting individual fluxes too with a forced-photometry pass
tried to forced-photometer multiple bands at once -- doesn't work
monkey around with stages, globals, exec, ...
more monkeying with stages()
Plot models for borked optimization attempts
tidy