Axis Systems on points normal to a surface

In Catia V5 r32, I would like to create a lot of axis systems on a lot of points that are on a surface. I would like the axis systems to be normal to the surface relative to each point.  It is strange that you can do something similar to this using points repetition but you need to use a curve in points repetition. I am have a predefined surface and predefined set of points. I tried creating a macro without success.