[SciPy-User] hemisphere function

[John]

Is there a function in scipy for a hemisphere?

Something like:

def hemisphere(X,Y,a):
    """ 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.??

--john

--