Dear MCusumano,
I nearly executed the program to perfection , but at the final stage, the Python Geoprocessing code wont work.
Possible Reason 1 -the centroid x and centroid y field names are POINT_X and POINT_Y, different from the code (point.x amd point.y ).
Solution : tried to change the code to match the field names.
Possible reason 2 - The building shape files have no cordinate system defined. Shall I try to define projection to the same CRS of the administrative boundaries/centroids ???
Thanks for the timely help again.