Moving to a multi-species implementation where a list of densities is maintained. Also improved performance by allocating large numpy Arrays at DTRW object creation
removed stray pdb call