Assign 3D data value based on a 2D index profile

Question

I have a 3D numpy array:

data0 = np.random.rand(30, 50, 50)

I have a 2D surface:

surf = np.random.rand(50, 50) * 30
surf = surf.astype(int)

Now I want to assign '0' to data0 along the surface profile. Which I know for loop can achieve this:

for xx in range(50):
    for yy in range(50):
        data0[0:surf[xx, yy], xx, yy] = 0

Data0 is a 3D volume with size of 30 * 50 * 50. surf is a 2D surface profile with size of 50 * 50. What I am trying to do is filling '0' from top to the surface (axis=0) in the volume

Here, 'for' loop is very slow, and it is inefficient when data0 is very huge. Could someone advise how to efficiently assign the values based on the surf profile?