import numpy as np
A = np.zeros((32,32,32))
print (A)

import numpy as np
from copy import deepcopy

size, radius = 7, 2
A = np.zeros((size,size, size))

AA = deepcopy(A)
x0, y0, z0 = int(np.floor(A.shape[0]/2)),
int(np.floor(A.shape[1]/2)), int(np.floor(A.shape[2]/2))

deb = 0

for x in range(x0-radius, x0+radius+1):
    for y in range(y0-radius, y0+radius+1):
        for z in range(z0-radius, z0+radius+1):
            deb += radius - abs(x0-x) - abs(y0-y) - abs(z0-z)
            if (deb)>=0: AA[x,y,z] = 1
            deb = 0