I have one point shapefile (POINTS) and one polygon shapefile (POLYS). i would like to create a routine that goes through each record in the polygon shapefile, selects all the points contained by that polygon, and the populates an attribute of each point (ATTR1), with an attribute value from the polygon dataset (ATTR2) from the polygon that contains those points.