The solver recognizes the constraints I set up in the log, but it apparently won't recognize it, and is therefore not giving me a bounded solution.

I tried this:

con jack: u2009[4]-d2009[4] = u_tot[4] - d_tot[4];

So when I print:
print u2009[4];
print d2009[4];
print u_tot[4];
print d_tot[4];


It gives me: u2009[4]:18391


d2009[4]: 0


u_tot[4]: 0


d_tot-4]: 0

So it obviously didn't obey the constraint, although it says it does in the log.