my case is only for 2d, but should apply to Nd as well.
 
It would be convienent if np.max would return a tuple of the max value and its Nd location indices.
 
Is there an easier way than just using the 1d flattened array max index (np.argmax) and calculating its corresponding Nd location?
 
Chris