changed function to calculate diversity trajectory for mBDI
changed function to calculate diversity trajectory for mBDI

fixes a problem occurring with identical TS or TE values among species

dsilvestro/PyRate,dsilvestro/PyRate,dsilvestro/PyRate,dsilvestro/PyRate