correct interpretation of shape of dx keyword

Previously the documentation had 'shape len(phi)', which wasn't right since numpy.ndarray.__len__ returns the phi.shape[0] which is an integer.  The code used len(phi.shape), a shorter equivalent for which is phi.ndim.