Is there a function in scipy for a hemisphere?
""" hemisphere over 2d X,Y """
h = np.sqrt(a-X**2-Y**2)
But with improved error checking, and perhaps returning an ma for the
NaN values, etc.??