Add keyword for unique number of trajectories

It would be nice to have a keyword for the unique number of position-momentum pairs. It would keep selecting new unique position-momentum pairs until it reaches the number and then it would keep selecting new initial conditions only from that until `npsamples`. Users usually want just a unique number of trajectories. 

I guess I will get to this feature once I need to use it again. 