Periodic smoothing spline on a 3D grid

Hello, I have a function (Vloc) defined on a 3D grid (rectangular regular grid). I would like to perform a smoothing spline on the function. I am currently using the following code :

x = {X,Y,Z};
[xx,yy,zz] = ndgrid(x{1},x{2},x{3}); 
y = Vloc;
noisy = y;
[smooth,p] = csaps(x,noisy,[],x);
mesh(x{1},x{2},smooth(:,:,5));

I would like to enforce periodic boundary conditions (which csaps does not automatically it seems). Does anyone know if a function can do that?