interpolation in 3D

I want to set the lame parameters for solving the lame equation to values I have generated in numpy. In 2D it is no problem at all, using interpolateTable. Is there anything comparable in 3D?

Right now i use the Brick domain only, so I can slice the domain and interpolate for each 2D slice (see below), but i'd prefer to use other domains as well, any ideas?

Thanks for your time in advance,
Martin

mydomain = Brick(l0=l0, l1=l1, l2=l2, n0=nx, n1=ny, n2=nz)
sedmu = numpy.array(...)
x = mydomain.getX()
for i in np.arange(nz+1):
    m = x[0].interpolateTable(sedmu[:,:,i], 0., l0/nx, x[1], 0., l1/ny, 1e14)
    dom = whereZero(x[2] - l2/nz*i)
    mu = mu + m * dom