You are close.
In the meteorite_landings data step, you need to:
1. Convert reclong and reclat to radians in X and Y.
2. Added some more annotate code.
x= reclong*atan(1)/45;
y= reclat*atan(1)/45;
function="pie"; style="psolid"; rotate=360; size=0.2; color='red';
You need to change MAPS.WORLD into WORLD and x=-long; y=lat; to be using unprojected values.
You are doing this since you are not using GPROJECT.
In gmap, you need to use WORLD instead of MAPS.WORLD. And you need to use ID id; choro id;
... View more