div_particles and globals.py

[olzama]

There were some objections to introducing globals.py; also, there may be a way of computing the value for div_particles that is more consistent with what is done elsewhere in the code (e.g. a function in choices.py which gets called appropriately).