It seems like the simplest solution is to add X and Y fields to the polygon feature class and calculate the X and & values as you mention.9.3 code would look something like:# Process: Add Field
gp.AddField_management(parcels_shp, "X_COORD", "DOUBLE", "", "", "", "", "NON_NULLABLE", "NON_REQUIRED", "")
# Process: Add Field
gp.AddField_management(parcels_shp, "Y_COORD", "DOUBLE", "", "", "", "", "NON_NULLABLE", "NON_REQUIRED", "")
# Process: Calculate Field
gp.CalculateField_management(parcels_shp, "X_COORD", "!shape.centroid.x!", "PYTHON_9.3", "")
# Process: Calculate Field
gp.CalculateField_management(parcels_shp, "Y_COORD", "!shape.centroid.y!", "PYTHON_9.3", "")
Once you have those added to the polys as attributes, why not simply use the identify tool?You can copy and paste the X and Y values from the tool.If you don't want to see all of the other attributes when you use the identify tool you can "hide" the other attributes from the identify tool by simply unchecking the other attributes in the Field tab of the layer properties.