print "maximum %s: %g" % ("xyz"[dim], atoms[-1].coord().data()[dim] - atoms[-1].radius)
Put the above in a file ending in ".py" and open it in Chimera and it will print the min/max xyz values accounting for radius in the Reply Log.
--Eric
Eric Pettersen
UCSF Computer Graphics Lab